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Aerodynamic Stability and Automatic 
Control 


The Fourteenth Wright Brothers Lecture* 


WILLIAM BOLLAY?T 
University of California at Los Angeles 


SUMMARY 


This paper summarizes the more recent developments in the 
field of aerodynamic stability and automatic control with par- 
ticular reference to the analytical and experimental methods that 
are used in the design of aircraft control systems. It traces the 
historical development of automatic control systems from the 
Watt Steam Engine Governor (1784) to the Sperry Gyroscopic 
Stabilizer (1910) and the automatic interceptor. It 
indicates that the analytical methods for studying automatic 


modern 
control systems are essentially the same as those presented by 
Routh in his Adams Prize Essay of 1877 and that the principal 
new step has been the reduction of these analytical methods to 
graphical techniques. The most useful of these are the ‘‘root 
locus” method of Evans, the amplitude and phase plots of Bode, 
and the vector response plots of Nyquist. These methods are 
illustrated by simple examples. The general methods of design, 
including the use of analog computers (simulators) and the 
analysis of flight-test data, are discussed. Some of the future 
potentialities of automatic controls in the application to aircraft 
are discussed: (1) automatic cruise control of propulsion sys- 
tems, (2) stabilization of an aircraft in the stall, (3) automatic 
stabilization of helicopters, (4) gust alleviation, and (5) automatic 
stabilization of seaplanes and hydroplanes. It is concluded that, 
while the past 25 years have brought revolutionary developments 
in aerodynamics and aircraft propulsion, the next 25 years will 
see comparably important advances in aeronautics by the appli- 
cation of automatic controls culminating not only in guided 
missiles, but also in transport aircraft that will be on automatic 


control from initial take-off to their landing at the destination 


(1) INTRODUCTION 


(A) Wright Brothers’ Experiments on Stability and 
Control in 1900 


p | {HE PRACTICAL ACHIEVEMENT of satisfactory sta- 
bility and control is probably the greatest contribu- 


Received May 8, 1951 

* Presented before the Institute of the Aeronautical Sciences 
in the U.S. Chamber of Commerce Auditorium, Washington, 
D.C., December 16, 1950 

t Visiting Professor of Engineering; 
tor, Aerophysics Laboratory, North American Aviation, Inc 


formerly Technical Direc 


tion of the Wright Brothers in the development of the 
airplane.* Their predecessors, Lilienthal, Mouillard, 
and Chanute, had carried out short soaring and gliding 
flights with primitive control systems; however, the 
Wright Brothers were the first to introduce effective 
aerodynamic control surfaces. Their first gliders and 
aircraft incorporated a tail-first configuration with a 
controllable elevator for longitudinal control, and they 
employed differential torsion of the wing tips to achieve 
lateral control (see Fig. I-1). Their first tests with these 
control systems were made during the summer of 1900 
on an unpowered aircraft, which they tested both as a 
glider and asa kite. These early flights in 1900 were the 
basis of their first airplane design. The topic of sta- 
bility and control is therefore a particularly appro- 
priate subject at the 1950 Wright Brothers Lecture 

a half century since the Wrights first verified their 
theories on stability and control. 


* Wilbur Wright in 1901 outlined the problems of the airplane 
very well when he stated,” ‘‘The difficulties which obstruct the 
pathway to success in flying-machine construction are of three 
general classes: (1) those which relate to the construction of the 
sustaining wings; (2) those which relate to the generation and 
application of the power required to drive the machine through 
(3) those relating to the balancing and steering of the 
Of these, two are already 


the air; 
machine after it is actually in flight 
to a certain extent solved. Men already know how to construct 
wings of aeroplanes, which, when driven through the air at suffi 
cient speed, will not only sustain the weight of the wings them 
selves, but also that of the engine and the engineer as well. Men 
also know how to build engines and screws of sufficient lightness 
and power to drive these planes at sustaining speed. Inability to 
balance and steer still confronts students of the flying problem, 
although nearly ten years (since Lilienthal’s success) have passed 
When this one feature has been worked out, the age of flying ma- 
chines will have arrived, for all other difficulties are of minor 


importance.” 
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Fic. I-l1. Wright’s gliding apparatus. 


(B) The Present—A Period of Transition—Subsonic to 
Supersonic Speed—Manual to Automatic Control 


At present, the airplane is going through another 
period of transition similar to that of a half century ago. 
The. power plants now available are powerful enough 
and light enough to permit airplanes to fly at speeds 
equal to, and greater than, the velocity of sound. The 
major problems that are encountered in flying in this 
regime are again stability and control problems. The 
transonic flow regime behaves in many respects like a 
stall regime. Below the so-called critical Mach Num- 
ber, as long as the angle of attack is below the stall 
angle, the airflow follows the contour of the wings. In 
the transonic regime, even though the free-stream 
velocity is subsonic, the local velocities of flow over por- 
tions of the wing are supersonic, and thus it is possible 
for shock waves to form. These shock waves have two 
important effects: (1) They change the pressure distri- 
bution over the wing and thus change the center of 
pressure of the wing. (2) Where they intersect the wing 
surface, they give rise to a steep pressure gradient that 
often leads to flow separation and thus buffeting. When 
the airplane speed is well up into the supersonic regime, 
the airflow again follows the wing contours; the principal 
shock waves are ahead and behind the wing, and the 
flow is smooth again. These phenomena are illustrated 
in Fig. I-2. The variation of the aerodynamic charac- 
teristics in the subsonic, transonic, and supersonic 
regime is illustrated in Fig. I-3 for the German V-2 
rocket. It will be noted that these coefficients vary 
smoothly in the subsonic and supersonic regime, al- 
though at different levels. In the transonic regime, the 
coefficients vary rapidly and, in some cases, relatively 
erratically as a function of Mach Number /, angle of 
attack a, and control surface deflection 6. The center 
of pressure on aerodynamic surfaces changes from the 
vicinity of the 25 per cent chord point at subsonic speeds 
toward the 50 per cent chord point at supersonic speeds. 
This is the principal reason for the large shift in the 
aerodynamic center. The hinge moments vary corre- 
spondingly because the center of pressure moves aft at 
supersonic speeds. An elevator that is balanced at sub- 
sonic speeds requires a large hinge moment to move it at 
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supersonic speeds. In sizing up the magnitudes of these 
moments, it must be remembered that hinge moment = 
Cy.(1/2)pV?.c.S; thus, the absolute magnitude of the 
hinge moment is proportional to the product of Cy and 
(1/2)pV?. 
that if the stick force per degree elevator deflection jg 
+1 Ib. at WM = 0.2, it is —94 Ibs. at WV = OS and 
+210 Ibs. at W = 1.2.* Such stick forces are too high 
for manual control of an airplane. 


Using Fig. I-3, we compute, for example, 


Thus, it is essential 
to use either power boost or, preferably, full power con. 
trol. The latter is the present trend. 

The increase in aircraft speeds brings with it another 


requirement—the necessity for decreasing the time of 
response of the pilot and the necessity for increasing the 


precision of his maneuvers. Consider, for example, the 
problem of shooting down a 600-m.p.h. bomber with an 
interceptor plane capable of a speed of 1,000 m.p.h 
The time available for interception is likely to be so 
short that the interceptor must be vectored by ground 
If the 
attack proceeds from the front, the relative speed be- 
If the intercep- 


control into the approximate attack position. 


tween the two aircraft is 1,600 m.p.h. 


* Such changes in trim could, for example, be produced by 
opening flaps or dive brakes. 
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- M= 0.90 
REGIONS \ 
M =0.95 
a.) FLOW FIELDS AND SHOCK WAVE PATTERNS 


IN THE LOWER TRANSONIC RANGE (M<!.0) 





M= 1.05 
SUBSONIC << 
REGIONS 
M=1.30 
b.) FLOW FIELDS AND SHOCK WAVE PATTERNS 


IN THE UPPER TRANSONIC RANGE (M>1.0) 


Fic. I-2._ Flow separation on airfoil surface.®** 
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AERODYNAMIC STABILITY 


tor pilot first sees the bomber 5 miles away, he has 11 
sec. to identify his target, to line up his airplane, to fire 
his guns, and to get out of the way. The likelihood of 
scoring a hit under conditions of such a high closing rate 
isextremely low. If he decides to go around to make an 
approach from the rear, even if he turns at 4g normal 
acceleration, it will require a 6-mile diameter circle to 
turn around and a time of 72 sec. During this time the 
bomber has traveled 12 miles and has therefore a good 
possibility of being under cloud cover. 

It is because of such considerations that the ‘‘auto- 
matic fighter’’ and, eventually, the guided missile will be 
required to intercept high-performance bombers. In 
reference 7 a typical interceptor mission, illustrated in 
Fig. I-4, 

“The automatic the 
navigate an interceptor into proximity of the 
track and attack the target, and navigate the intercep- 
tor back to its base. . .the target will originally be 
located by ground-based, The in- 
formation obtained as to bearing of the target will be 
manually set into the system by the pilot, and the air- 
craft will then take off, climb to cruise altitude, and 
navigate to the vicinity of the target At this 
point, airborne radar in the interceptor will take over, 
search for the target, identify if, and begin to track it. 
While tracking, the system will automatically maneuver 
the airplane into the attack position. At the proper in- 
stant, calculated to yield the best ‘kill’ probability, 
guided missiles or rockets will be automatically re- 
After sufficient time has elapsed for the system 


is described as follows: 
control system of future will 


target, 


long-range radar. 


area. 


leased. 
to determine the results of the attack on the target, 
The system then 
either to another 


automatic evasive action will begin. 
will return to automatic navigation, 
target area or to the home base. At the home base, an 
instrument landing pattern will be automatically 
entered, and the plane brought to touchdown. The 
action in any phase of this sequence will be completely 
automatic; the pilot's job will be to determine when to 
switch from one phase to the next.” 


The development of such automatic control systems 
involves two important types of technical problems: 
The combination of the many different elements 
into a single, well-integrated system. The airplane with 
its autopilot are only two of the elements of this system. 
These two elements must be capable of responding to 
the directions of (a) the ground radar, which vectors 
the airplane into the vicinity of the bomber; (b) the 
airplane's target seeker, which steers the airplane into 
firing position and fires the guns or rockets; and (c) 
the automatic-ground-approach system (Fig. I-4+). The 
overall system must be stable over the wide range of 
involved. The system must reject 
input signals that behave like tur- 
The system 


speeds and altitudes 
the spurious ‘‘noise”’ 
bulent fluctuations about a mean velocity. 
must be precise enough to provide accurate fire control 


during the attack. The landing approach has to be pre- 


AND AUTOMATIC CONTROL 571 
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Fic. I-38. Variation of aerodynamic characteristics of the 
V-2 with Mach Number. 
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Fic. 1-4 


cise enough so that the airplane does not exceed its de- 
sign landing-load factors upon touch-down. 

The requirement for rapid automatic control systems 
with reasonable control forces will require much closer 
teamwork between aerodynamicists and autopilot de- 
signers than ever existed in the past. There are good 
possibilities for new types of aerodynamic controls which 
were not necessary previously. Spoilers may be used 
for rapid control, leaving the trimming to the more con- 
New types of aerodynamic bal- 


ventional surfaces.” 
ances must still be invented which will be effective over 
a wide range of Mach Numbers, angles of attack, and 
angles of yaw. The dynamic analysis of these systems 
must be carried out as a joint enterprise by the aerody- 
namicist and the autopilot designer. 

(2) The achievement of a high degree of reliability in 


this complicated system.* This can be achieved par- 


* An excellent review of the problems involved in achieving a 
reliable powered control system for aircraft is given by Lyons.** 
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tially by careful initial design, such as the use of dupli- 
‘ate components for critical elements and the use of 
miniaturized “‘plug-in’’ boxes that can be easily checked 
and replaced. The high degree of reliability will be 
finally achieved only as the result of an intensive ground 
and flight-test program in much the same way as the 
current aircraft engines have evolved. 

After this has been achieved, the automatically con- 
trolled airplane will also have important peacetime 
potentialities. The transport aircraft of the future may 
well be on automatic control from the moment of take- 
off to the automatic landing at its destination. Such a 
system would free the air transport of its most serious 
handicap—the weather. The human pilot would serve 
mainly to oversee the flight instruments and to take 
over in case of an emergency. 


(II) STABILITY AND CONTROL PROBLEMS OF THE PAST 
DECADE 


During the past decade, airplanes have moved into 


the transonic regime. They have thus pierced what 
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Fic. IIl-1. Typical effect of flow separation on downwash angle. 
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used to be called the sonic wall. This was a wall not 
only from the standpoint of requiring tremendoys 
powers to penetrate, but even more, it was a region of 
great ignorance regarding the stability and control 
problems that it involved. 


The fighter airplanes developed during the early 
1940's were the first airplanes operating at speeds high 
enough to encounter these ‘‘compressibility effects,” 
During steep dives these airplanes reached Mach Num- 
bers of the order of 0.70 and experienced many severe 
stability and control difficulties. Buffeting occurred at 
these high speeds. This was associated with a flow 
separation from the wing when a shock wave was 
formed on the upper surface (see Fig. I-2). In addition, 
the airplane experienced a severe change in trim, which 
was associated with the unsymmetrical flow separation 
from the wing and a reduction in downwash over the 
tail (Fig. II-1). 
tribution gave 


The resulting change in pressure dis- 
a variation in elevator stick force, as 
shown in Fig. II-2. It is noted that, at 1J = 0.72, a re- 
versal in stick force occurs, followed by a severe diving 
The 
pilot, realizing that the pull forces were exceeding his 


or ‘‘tucking’’ moment as Mach Number increased. 


strength capabilities, normally used the elevator trim 
tab to trim the airplane nose up. This displaced the 
stick force curve as shown in the dotted portion of Fig. 
II-2. 


tudes, the speed of sound increased because of the higher 


As the airplane continued its dive to lower alti- 
ambient temperatures. Consequently, the airplane 
rather rapidly decreased Mach Number, and, if the pilot 
had used too much nose-up trim tab, the airplane 
rapidly pulled out of the dive, often with disastrous re- 
This was alleviated in the P-38 and P-47 air- 
These 


sults. 
planes by the addition of dive-recovery flaps. 
flaps tended to increase the downwash over the tail and 
to induce separation on the lower surface to balance 
the separation on the upper surface of the wing. The 
effect of these flaps is shown in Fig. II-3. In the design 
of the newer jet fighters, these troubles are first of all 
minimized by design changes —the use of sweepback of 
wings, the use of extremely thin wings, or the use of low 
aspect ratio. Fig. Il-4 shows the effect of reducing 
wing thickness. Fig. [1-5 shows the effect of sweepback. 
It will be noted that the effect of these design changes 1s 
to postpone the “‘transonic bump.” 


At supersonic speeds the effectiveness of trailing-edge 
elevators is much reduced (Fig. II-6), and thus it has 


become desirable to use an all-movable stabilizer.” 


The control moments on such a movable stabilizer be- 


come so high at high speeds that either power boost or 
power controls are essential. With power controls the 
surface is actuated directly by a hydraulic actuator or 
motor-driven screw jack. This has the additional ad- 
vantage that the undesirable reversals in hinge moments 
are not transmitted back to the pilot. The pilot's stick ! 
merely controls the hydraulic actuator. In order to 


provide a simulated feel in the stick, it is connected to a 
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AERODYNAMIC STABILITY 


spring for centering. In addition, the aerodynamic 
loads are simulated by a ‘‘q”’ 
sensitive to free-stream dynamic pressure; it increases 


the effective spring stiffness by shifting the spring 


actuator, a device that is 


bungee connected to the stick. 

In order to limit the accelerations that the pilot will 
apply to the aircraft, bob-weights may be connected to 
the control system. When the airplane experiences a 
normal acceleration, these bob-weights exert a moment 
on the stick in a direction tending to reduce the control- 
The effect of such bob-weights is 
This airplane 


surface deflection. 
shown in Fig. II-7 on a P-51 airplane. 
was overloaded by the addition of a fuel tank behind the 
pilot to increase its range. This made the airplane un- 
stable longitudinally, so that when the pilot applied a 
normal acceleration to the airplane by pulling back on 
the stick, an increase in acceleration required less stick 
force or even a push force. The addition of the bob- 
weights gave the airplane a more nearly normal feel. 
Another family of stability and control problems 
which showed up during the past 10 years was associ- 
i.e., the deformation of the 
Some of these 


ated with aeroelasticity 
aircraft structure due to the air forces. 
can be cured only by a stiffening of the structure. 
example, during high-speed dives, the elevator contour 
would tend to bulge outward and, in some cases, would 
The bulged surfaces resulted in a 
Vented sur- 


For 


rupture completely. 
decreased stability or stick-force gradient. 
faces and strengthened structures usually eliminated 
this problem. 

A typical problem of control-force reversal occurs 
when an aileron deflection, say, up, tends to twist the 
nose of the wing in the same direction—i.e., up. As the 
dynamic pressure is increased, the rolling effectiveness 
of the aileron is continuously decreased until, at a cer 
tain speed, the aileron effectiveness becomes zero. 
This is shown in Fig. II-8. A solution to this problem is 
the use of spoiler ailerons that produce a given rolling 
moment with less twisting moment on the wing. Thus, 
they are effective in increasing the speed at which 
aileron effectiveness is reduced to zero. 

Another aileron-control problem encountered at 
high Mach Numbers under high normal accelerations on 
a sweptback wing is ‘“‘upfloating’’—i.e., a deflection of 
the control cables permitting both ailerons to deflect 
upward. If the ailerons are far behind the c.g. (as for 
swept-wing aircraft), this acts like an additional eleva- 
tor deflection and therefore results in a still greater 
acceleration on the airplane. One cure for this problem 
is a mechanical restraint of the aileron. 

A dynamic instability that has been encountered on 
high-speed aircraft is ‘‘snaking.’’ This is a rapid direc- 
tional oscillation at the short period lateral frequency 
apparently caused by a wide variety of elements—flow 
separation in the vicinity of the vertical tail, fuel slosh 
ing, and low lateral damping associated with aeroelas- 
ticity or high-altitude flight. This builds up a vibration 
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Two methods ol 
The first 


involving bending of the fuselage. 
fixing this vibration have been employed. 
adds trailing-edge strips to the rudder to increase its 
However, this 
Another 


effectiveness at low deflection angles. 
‘fix’’ also increases rudder forces rapidly. 
means of eliminating this vibration has been by an 
automatic control system—a yaw damper that is 
actuated by a rate-sensitive gyro and which damps the 
vibration by deflecting the rudder in a direction oppos- 
ing the motion. This method, described in reference S2 


and shown in Fig. II-9, has been found effective. 
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(111) DEVELOPMENT OF AUTOMATIC CONTROL SYSTEMS 


The theory of dynamic stability of aircraft has until 
recently been considered largely as an academic sub- 
ject. The basic reason for the low practical importance 
of the earlier dynamic stability calculation was the 
artificial boundary conditions that were ordinarily as- 
sumed for the control system—-i.e., either that the con- 
trols were locked in a trimmed position or that the con- 
trols were completely free. Both of these assumptions 
are unrealistic, since actually the pilot sits at the con- 
trols and, if a disturbance increases the airplane attitude 
and thus reduces its speed, he counteracts this distur- 
bance through his controls. The difficulty with the 
mathematical analysis was that it was impractical to 
include the human pilot as a part of the dynamical 
system. 

With the development of automatic control systems 
this situation changed. The autopilot system and its 
input information (no matter whether it comes from a 
radar signal or a gyroscope system) is susceptible to a 
dynamic analysis. This dynamic analysis is ordinarily 
based partially on theory and partially on performance 
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SCIENCES 
tests of the actual equipment. The overall system, 
consisting of the aerodynamic, inertial, and elastic ele- 
ments of the airplane combined with the electrical, 
mechanical, and hydraulic elements of the autopilot, 
forms a closed-loop network. The analysis of such 
feedback networks forms the next subject of our dis- 
cussion [Section (IV)]. This entire system must be 
carefully engineered as a unit in order to combine 





these many elements into an overall precision device. 
This process is called the synthesis of the airplane-auto- 
pilot system. 

The development of automatic control systems has 
progressed far, particularly during the past 10 years, 
Automatic controls have, however, been known for a 
One of the first was the flyball governor for 
Routh (1884) 


long time. 
the Watt Steam Engine (about 1754). 
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AERODYNAMIC STABILITY 


analyzed the stability of this system’ and noted that, 
unless damping was introduced in the spindle on which 
the flyballs move up and down, the system would hunt. 
A great improvement in the achievement of uniform 
rotation was the invention (1843) by Siemens”! of a 
governor that controlled the steam valve of an engine 
by the error signal between a uniform angular velocity 
and the actual angular velocity. This instrument, the 
chronometric governor (Fig. III-1), was applied by Sir 
G. B. Airy, the Astronomer Royal, for regulating 
astronomical and chronographical instruments. This 
governor operated on the principle of what is now called 
a simple feedback servo network (Fig. III-3). 

The invention and description of the hydraulic servo- 
motor occurred almost simultaneously in England by 
Brown! in 1871 and in France by Farcot" in 1873. 
Lincke'® in Germany (1879) presented an excellent sur- 
vey of many mechanical feedback servosystemis. He 
showed, among other things, the analogy of these sys- 
tems to the human nervous and muscular system and 
indicated the future potentialities of servosystems for 
the automatic control of machines and factories. He 
was thus an early exponent of what Wiener®! now calls 
“eybernetics."’ The application of servos to the regula 
tion of steam turbines was studied experimentally and 
1893 with particular 
Later summa- 
ries of regulators were given by Tolle’ in 1921, by 
Winsch” in 1930, by Engel* in 1944, and by Hanny" in 
1946. 

A new element involved in servo-control was the 


theoretically by Stodola” in 


reference to the conditions of stability. 


combination of the mechanical feedback and power 
amplification. The application of such servomecha- 
nisms to the directional control and stabilization of ve- 
hicles started with the control of torpedoes by White- 
head (1862). It is interesting to note that Sir Hiram 
Maxim* considered the use of a gyro-controlled servo- 
system for the stabilization of his first airplane in 1891. 
The practical application of a gyro for the stabilization 
of an airplane was started in April, 1910, when Elmer A. 
Sperry? and H. C. Ford installed a gyro in a monoplane 
belonging to Mr. Stanley Y. Beach. This device was 
known as the ‘“‘Gyroscopic Stabilizer.’’ These experi- 
ments were continued by Lawrence Sperry in coopera- 
tion with the U.S. Navy on a confidential basis. In 
1914 a public exhibition was given in Paris of a Curtis 
flying boat under the control of an automatic pilot (see 
Fig. III-2). For this exhibition, L. B. Sperry received 

* According to Lanchester! (Aerodonetics, p. 424, 1908): “In 
1891 Sir Hiram Maxim proposed to secure the longitudinal 
stability of his flying machine by means of a gyroscope arranged 
to actuate horizontal bow and stern rudder, through the medium 
of a steam or pneumatic relay, in a manner similar to the main 
tenance of course in the Whitehead torpedo. As this machine 
never got beyond the captive stage, it may be presumed that the 
arrangement in question was not carried into effect.”’ 

+ This information was kindly made available by Mr. P. R 
Bassett, President of the Sperry Gyroscope Company, from his 


company’s private files 
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ic. III-1 Chronometric governor (Siemens,*' 1843) 


the Grand Prize; he had demonstrated to the world, for 
the first time, the practicability of automatic flight. 
The analysis of servo-control systems for the steering of 
ships was made by Minorsky" in 1922. This study led 
to the application of autopilots for the control of ships. 
Feedback control systems were also applied to elec 
trical networks, and the theory of feedback systems was 
expanded particularly in connection with the design ot 


electrical feedback amplifiers. Nyquist’? (1982) im a 
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Fic. L11-2. Prize-winning flight with Sperry gyro stabilizer in 
Curtiss F-boat, Paris, 1914. L. B. Sperry, pilot, standing in 
cockpit with hands off controls while mechanic walks out on 
wing. 
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Fic. I1I-3. 


Simple feedback loop. 


noteworthy paper applied the techniques developed by 
Cauchy and Routh to the feedback amplifier. Bode* ex- 
tended these techniques and developed a simple graphi- 
Hazen'' presented a summary of the 
The most wide- 


cal analysis. 
theory of servomechanisms (1934). 
spread application of these techniques came during the 
early 1940's in connection with the development of anti- 
aircraft gun directors. This work is well presented in a 
number of excellent books '* ° by the M.I.T. staff. A 
simple, but effective, extension of these techniques was 
developed by Evans.’ The application of these tech- 
niques for the study of the dynamic performance of an 
airplane-autopilot system is illustrated in Section (IV). 

Probably the major development in aircraft stability 
during the past 10 years has been the evolution of 
analytical and experimental techniques that permit an 
eugineering calculation of the motion of an airplane 
when under the control of an automatic pilot. None of 
these techniques appears particularly new or novel when 
viewed by itself. Routh in 1878 used the fundamentals 
of all of these theoretical techniques in his stability cal- 
culations. Thus, to the mathematician, there 1s little 
new represented here. The new element is largely that 
during the past 10 years the engineers have appro- 
priated some of these powerful tools of the mathematical 
physicist, that they have learned to think in terms of 
them, and that they have learned thereby to conceive 
an overall control system that is able to perform not only 
stably, but accurately according to the input com- 
mands. 

The principal techniques that have made this develop 
ment possible are the following: 

(1) The Laplace transform method of solving differen- 
lial equations..This technique, which is also called 
operational calculus, was popularized by Heaviside** 
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(1893), who introduced it to electrical engineering, 
The mathematical fundamentals of this technique are 
more than 150 years old. Its great usefulness is not only 
its function as a shorthand method of solving differential 
equations, but even more, its aid in helping the engineer 
in synthesizing a suitable control system. The elements 
of this method, particularly as they apply to the design 
of automatic control systems, are summarized in Ap- 
pendix (A). 

(2) The graphical methods of analyzing the dynamic 
performance of a system, particularly the diagrams often 
associated with the names of Nyquist, Bode, and Evans. 
These methods are summarized in Section (IV). 

(3) The experimental techniques of simulating the 
physical dynamical system by an equivalent electrical 
network. 

Last but not least in importance in the evolution of 
these techniques has been a series of excellent sum- 
mary reports and textbooks® !* 1% that have intro- 
duced these new techniques to a wide circle of engi- 
neers. 

ANALYSIS OF AUTOMATIC CONTROL 
SYSTEMS 


(IV) DyNAmIc 


(A) Introduction 

The general approach in the development of auto- 
matic control systems is similar to that of the aerody- 
namicist in the development of an airplane configura- 
tion. It involves the following steps: 

(1) Use of a simplified theory with empirical correc- 
tions in order to determine the general circuit design. 
This step is commonly referred to as the synthesis of the 
system. 

(2) Use of models and application of the theory of 
similarity in order to obtain a better approximation of 
the performance of the system. Ordinarily, the aero- 
dynamic and inertial elements of the airplane are re- 


placed by a “‘simulator,”’ an electrical network whose 
constants are adjusted so that it obeys the same dif- 


ferential equations as the airplane on a one-to-one time 


In the early phases of the design, the servo- 


basis. 
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AERODYNAMIC STABILITY 
system elements may also be replaced by an equivalent 
electrical network. In the later phases of design, the 
actual servosystem elements are included in the cir 
cuit. Thus, the nonlinear characteristics of the hy- 
draulic systems or the effects of mechanical backlash 
may be correctly introduced.* 

(3) Flight tests to determine the actual performance, 
including the nonlinearities of the elastic and aerody- 
namic elements, interactions between longitudinal and 
lateral modes of motion, and the effects of environment 
(vibrations, pressure, and temperature). 

Many of the research methods of aerodynamics also 
carry over into the analysis of automatic control sys 
tems. Thus, the theory of turbulence laid the founda- 
tion for the theory of “‘noise”’ in electrical circuits. 

We shall discuss first of all the simplified theory in 
which the airplane and its control system are approxi- 
mated by linear differential equations with constant co- 
The actual airplane and its control elements 


It is ordinarily, however, a good 


efficients. 
are nonlinear systems. 
approximation to linearize the equations for small os- 
cillations about an equilibrium position. When the 
nonlinearities become pronounced—for instance, when 
an airplane goes through a maneuver involving large 
changes in angle of attack —then the simple analytical 
methods no longer yield quantitatively correct answers, 
and they will only be used to obtain a qualitative idea 
of how to synthesize the system. Ordinarily a combina- 
tion of analytical and experimental methods is then 
most conveniently employed in order to obtain the 
quantitative estimates of the stability and the response 
of the system. 

The linear differential equations of the system are 
solved by the method of Laplace transforms {see Appen- 
dix (A)]. This technique converts the differential equa- 
tions in terms of the time variable ¢ into algebraic equa- 
tions in terms of a complex frequency variable p = 
\ + w. The solutions of the differential equations are 
of the type Ce” = Ce .e'“’. Thus, it is seen that X, the 
real part of p, corresponds to the damping coefficient of 
the solution and that w, the imaginary part of /, corre- 
sponds to the angular frequency. 

In order to maintain a good physical insight into the 
behavior of the system, it is ordinarily desirable not to 
proceed with the solution in a purely formal manner. 


* The most commonly used simulators are of the analog com- 
puter type, in which amplification and integration are accom 
plished electronically. Such a computer is shown in Fig. IV-1. 
This type is particularly well adapted for tests at the true time 
scale in connection with autopilot elements for frequencies be- 
tween 0.005 and 20 cycles per sec. The straightforward electrical 
analog computer, consisting of resistances, capacitances, and 
inductances, is suitable for calculation of autopilot responses but 
not for connection with the actual autopilot hardware. The 
reason is that inductances and capacitances must be constant 
and have low losses; thus, the magnitude of LC is limited to 
approximately 10~‘, and correspondingly, the frequency must be 
20 cycles per sec. or greater. This type of computer operates 
therefore, at a much shorter time scale than the actual autopilot. 
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Fic. IV-8. Block diagram of a second-order system. 


Ordinarily it is helpful to break down the system by 
means of a block diagram into its major components. 
Thus, a simple airplane roll control system consists of 
the major components shown in Fig. IV-2. 

(1) An input signal directs the airplane to roll through 
an angle gj. 

(2) The measured roll angle ¢, is compared with the 
input by a potentiometer and the error angle ¢ = 
¢; — ¢- 1s determined. 

(3) This error angle ¢, is the input to an amplifier and 
computer. Its output is gc. 

(4) This signal ¢g,, is introduced into a hydraulic 
aileron servo, resulting in an output aileron angle 6,. 

(5) This aileron angle 6, is applied to the air frame, 
resulting in an actual roll angle gp. 

(6) This roll angle gp is measured by a roll gyro, re 
sulting in an output signal ¢, (which may be slightly in 
error due to instrument imperfections). 

(7) This gyro signal ¢, is compared with the input ¢, 
as in step (2). 

Each of these block components may then have a 
further breakdown into its block elements. 

The characteristics of any block are determined by 
the input y,(¢), and the output yo,(t). A more con 
venient description is obtained in terms of their Laplace 


transforms (Fig. IV-3). Thus, we shall define 
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transform of a linear velocity, Vi(p) = p.Yi(p), and 
CANTILEVER SPRING the output may be the transform of a displacement, 
A ae I Yo(p). Then the ratio 
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First-order system: 


Vi(p) = Laplace transform of 


input signal y;(f) 


output signal yo, (f) 


The output yo(¢) and, correspondingly, its transform 
Y'y,(p) would in general consist of two parts: 


vont) - yvo(t) + Vo(t) 


(a) The output yo(/) and its transform 


p / e yo(t) dt 


which are due to the input signal y;(f). 
(6b) The output y,(¢) and its transform Y,(p), which 


Yo(p) = 


are due to the boundary conditions at ¢ = 0 — e« —that 
is, to y,(O), yo’(O), vn” (0), ete. 
We define 

Vu(p)/ Vip) = K.G(p) = ratio of output to input 
(IV-1) 

where 

Y,(0)/V,(0) = lim Yo(p)/¥i(p) = A = “gain” 

a (IV-2) 

and 

G(p) = (1/ AK) Vo(p), Vi(p) = “transfer function” 
(IV-3) 

Note that the input and output do not necessarily have 


the same dimensions. Thus, the input may be the 


Yo( ) Yo 
P = iP) = K.G(p) 
Vi(p) pb. Yi(p) 
and 
Yo 
- p) = KG,(p) 
Yi(p) 


where Gi(p) = p.G(p). 

Fig. IV-3 illustrates the definitions of the gain and the 
transfer function for a second-order system. It is seen 
from Fig. IV-19 that the gain of a series circuit is the 
product of the gains of the elements and that the overall 
transfer function is the product of the transfer functions 
of the elements. 

The basic elements into which these blocks are re- 
solved are first-order or second-order systems. In our 
simple analysis, we shall therefore study particularly 
the characteristics of the first-order and second-order 
systems and the methods for combining the solutions of 
these simple systems. The purpose of the analysis will 
be two-fold: (a) to determine the range of parameters 
for which the overall system is stable and (b) to deter- 
mine: the accuracy of control—i.e., to compute the over- 
all output response for a given input. 


(B) Analysis of First-Order Systems 


Consider a cantilever spring (Fig. IV-4) which has a 
The spring stiffness is k, the 
The end A of the 
cantilever spring is moved upward; its position being 
given by y,(¢). End B lags behind the input due to the 
damping; it experiences the displacement yo,(t). (This 
simple case approximates, for example, the motion of a 


dashpot at its end B. 
damping constant of the dashpot is c. 


recording pen having viscous friction.) Neglecting the 
mass-acceleration forces, the forces acting on the dash- 


pot are 
= R(y; — Von) 


Spring force—upward 


Damping force—upward = —c.dyo,/dt 


Thus, the equation of motion is 


C(dyon/dt) + Rym = Ry; (IV-4) 
or 

T,.(dyon/dt) + You = Yi 
where 7, = c/k. The boundary condition is yo = 


y(O) at t = 0 — e where y,(0) is a constant. The 


Laplace transform of this equation is 
(T\.p + b?. Voo(p) = Y;(p) + T\.p Vp (0) 


and thus 
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| Vi(p) + ie (0) (IV-5 
= -},(p =< * VY) a) 
Yoo(P) Tip 1 i Tip l - 
Yo(p) Y,(p) 


The term }(p) represents the Laplace transform of the 
output motion which is due to the input signal y,(¢). 
The term Y,(p) represents the Laplace transform of the 
output motion which is due to the initial boundary con- 


dition y,(Q). 
The transfer function in this case is found from 


KG(p) = Yo(p)/Vilp) = V/ (Mp + 1) (IV-6) 
and thus 
K = |] 
G(p) = 1/(Mp + 1) 


The poles and zeros* completely characterize the trans- 
fer function G(p). We shall represent G(p) in the p- 
plane by indicating each zero with a cross and each pole 
with a small circle. The transfer function G(p) = 
1/(T\p + 1) has no zeros, but it has a simple pole at 
p = — 1/7), as shown in Fig. IV-5c. 

Special case: y(t) = h(t) 

For the special case of a unit step input in displace- 
from Eq. (B-9) | 





ment, we have 
yi(t) = h(t) 
Yi(p) = 1 


The corresponding solution for the first-order Eq. 
(IV-6) is [from Eq. (B-8) | 


Yo(p) 1/(Tip+ - 
-t/7 (IV-7) 
vo(t) = (| l =~ \f 
and from Eq. (B-1) 
op) = [T »(O Tip - ) 
Yo(p [7\.p y0(0) (Tip + L)I (IV-8) 
y(t) = w(O).e f 
and thus the complete solution is 
yo(t) = (1 -—e' Mt) + yr(O).e és 
yo(t) V(t) (IV-9) 
The term y,(t) = y,(0).e~""' represents a damped 


subsidence due to the initial boundary condition y,(0). 
This term enters into the analysis of automatic control 
systems only when it is necessary to compute the re- 
sponse of a system that is not initially at rest and in 
equilibrium. 

The term y(t) = 1 — e~’"' represents the output 
motion due to the step input /(¢). 
in Fig. IV-4b. 
proach to the input. 


This term is shown 
It is represented by an exponential ap- 
It indicates that the output lags 


* We define a simple pole of a function as a point at which the 
We define a simple 
sero of a function as a point at which the function becomes zero, 
Thus, the function Y¥(p) = p/(p - 
simple pole at p = a and a simple zero at the origin p = 0. 


function becomes infinite, like 1/x as x ~ 0. 


like x as x —> 0. a) has a 
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c) DIAGRAM OF END-POINTS OF VECTOR G (iw) AND sis) 


Simple lag circuit: G(iw) = Me? 1/(1 + tw7 
response to sinusoidal input 


Fic. 1V-6 


behind the input, that at a time / = 7), the output has 
reached 63 per cent of its asymptotic deflection. The 


error signal is 


y(t) = wilt) — w(t) = e 


We note that the error y,(t) ~ 0 ast ~ o and that 


the initial slope of the output is 
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dyo(t | l t ‘| l 
eS : ware 
dt —: 7 t=0 7 


Special case: y;(¢) = Vim .e 


THE 


The special case of a sinusoidal input with frequency 
w is most conveniently represented by writing the input 
yi(t) = Vim .e" Thus, the 
real part of the input is y.,,. cos wf, the imaginary part is 
Yim. sin wt. The Laplace transform of the input is 
Yi(p) = Yin. p/(p — tw) from Eq. (B-2). 
sponding solution for the first-order equation [Eq. 
(IV-6)] is 


t wt - 
where e“” = cos wt + 7sin wf. 


The corre- 


Y,(p) = [Tip .¥o(0) | (Tip 4. 1)) ; 
e -1¢ 
y(t) = w(0).e07" ( (1V-10) 


p | 


Yo(p) = Vim° = ae 
: p—-wT\ptl 


(IV-11) 


— Vi m 


> aa 
1+ wT] 


V . 
t/1 7 tm toot 
ie eh 
l + tw] 
—— 


| vo(t) rs 


Vo(t) 





[vol t) = 


The terms },(p) and y,(¢) are the same as obtained in 
the previous special case. The solution yo(t) may be ob- 
tained from Eq. (A-15) or by expressing Vo(p) in partial 
fractions of the form 


p fie Pp 


Yo = Y,(0) A ” y 
(p) + Ay p+il T,) + 


p — tw 


We then have from Eq. (A-9): 


Y,(0) = 0 
A, = —Vin (1 4. lw ]’;) 
As = yin/(1 + iwT)) 


We note that the solution yo(‘) has as its first term 
pO ai damped subsidence. This term represents a 
transient disturbance associated with the initial applica- 
tion of the force. 

The second term of yo(t) is vim/(1 + tw7\)e. This 
term represents a steady sinusoidal vibration at the 
same frequency w as the input. 

The ratio 

eee. (IV-12) 
yilb) 1+ wl 
It will be noted that this ratio is the same as would be 
that is, 


obtained by writing ‘w instead of p in G(p) 
G(tw) = 1/(1 + tw7)) 
and thus for a sinusoidal input 
[yo(t) |steaay = G(tw) . yi(t) 
This result follows as a special case of Eq. (A-15). It is 


convenient to express the complex vector G(iw) in 


terms of its magnitude .\/ and its phase angle 6 —that is, 


G(iw) = Me” 
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In the above example of a first-order system we then 


find 
l L ~ dot 
Me = —— = at 
+ lw l + wT? 
that is, 
M = 1/V1+?T7;7 (IV-13 
tan @é = —wl| (1V-14) 


The above relations indicate how well the output is 
able to follow a sinusoidal variation in the input: 
corresponding experimental records are shown in Fig, 
IV-4e. 
output y(t) follows the input y;(f) in magnitude and 


We note that at very low frequencies the steady 


phase, at higher frequencies the ratio .\J becomes less 
than one and the output lags the input, and at very 
high frequencies ./ — 0 and the phase angle lags by 
90°. The variation of .V and @ with the dimensionless 
It will be 
It is con- 


frequency “« = wT) is shown in Fig. IV-6a. 
noted that for large values of u, \J =~ 1/u. 
venient to plot the curves of .V/ vs. u on log-log paper 
(Fig. IV-6b). 
Then we note that for large values of u, since 


This corresponds to plotting logy MJ vs. 
logio u. 
logw.W = logw l/u = — logy u, the slope becomes —1; 
for small values of u, logy./ starts with a slope of 0. 
Thus, the first-order system can be approximated by the 
two dotted straight lines on this diagram. 

In the acoustic and electrical literature it is customary 
to plot 20 logy) ./ in order to convert the amplitude units 
into decibels. A doubling in frequency is called an 
octave, and thus the region where our curve of logy M 
vs. logw u has a slope of —1 is denoted by a slope of 
—20 logy 2 = —6.03 db. per octave. In order to read 
the electrical literature, it is necessary to know this 
definition; there appears to be little reason, however, to 
introduce this complication into the aerodynamics 
literature. 

It will be noted that on this logarithmic plot the 
approximate curve logy J/ goes through 0 at u = | 
that is, at w = 1/7. 
perimental response of a first-order system and plot it 


If we therefore measure the ex- 


as indicated above, the delay time 7) is easily esti- 
mated by noting the frequency at which a straight-line 
approximation for large values of w crosses the axis. 
The above described methods of plotting the response 
curves to a logarithmic scale and the corresponding 
approximations were pointed out by Bode.’ 

Another convenient method of representing G(iw) = 
Me" is directly as a complex vector of magnitude W and 
argument 6. Since both of these functions have al- 
ready been plotted in Fig. IV-6b, it is simple to con- 
This diagram for a 
We note that 
this diagram consists of a half circle, starting at G(iw) = 


struct the vector diagram of G(iw). 
first-order system is shown in Fig. 1V-6c. 
1 at w = O, going through G(ww) = 1I/(1 + 72) = 
(1/Y2)(1 — 7) at w7| = 1, and ending at G(iw) = Oat 
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yervoutace “| S(P)F'D $e) poe $;(t)= ANGULAR | Vo(t)= VOLTAGE | 
INPUT ELECTRIC MOTOR a RATE OUTPUT | 
ROTOR GYRO 
p- PLANE p- PLANE 
(a) (b) (a) (b) 
w 
1 2 4 6810 +90° 
—— TIM 
ei | til 
M | | M 
A | Went A 8 
-90°! 
| 
0! Lt iti 
! 2 4 6 810 
w 
tc) (c) 
; - ; ire oe ie _ k di: Fic. [V-8. Differentiator circuit: G(iw) = twA. (a) Block 
Fic. 1V-7. Integrator circuit: G(iw) = A/iw. (a) Block dia- diagram of G(p) = Ap. (b) Representation in p-plane (simple 
gram of G(p) = A/p. (b) Representation in p plane (simple zero at p = 0). (c) Response to sinusoidal input: G(iw) = 
pole at bs = 2 c) Response to sinusoidal input: G(iw) = Me?®. M = Aw = magnification ratio (logarithmic scale), 6 = 
Me, J = A w = magnification ratio (logarithmic scale), 0= +90° = phase angle, and w = exciting frequency (logarithmic 
—90° = phase angle, and w» = exciting frequency (logarithmic scale) ; 
scale 
wT; = ©. This vector representation of the re- \n electric motor whose speed d¢o/dt is proportional 


For many 
It will 
be noted that, in the case of a first-order system, this re- 


+ lw] (Fig. 


sponse curves was introduced by Nyquist. 
purposes it is more convenient to plot 1/G(iw). 


sults in the straight line 1/G(iw) = 1 

IV-6c). 
Summary. 

for a simple first-order system we note the following: 


I (Tip + 7. 


Summarizing the important conclusions 


(1) Its transfer function is G(p) = 
which has a simple pole at p = —1/7). 

(2) For a step input the output approaches the input 
exponentially with the time constant 7}. 

(3) For a sinusoidal input there is a transient and a 
steady-state output. The ratio of steady-state output 
to sinusoidal input, which is most important for the 
stability analysis, may be written down directly from 
the transfer function 

ee, Ps Cis) = — 


1+ iw] 


Vi (0) 


(4) The transient disturbance due to the initial 
boundary conditions at ¢ = 0 — ¢ and the transient dis- 
turbance associated with the initial application of the 
force are damped out like e~'’"'. These terms are im- 
portant for calculation of the response curve but not for 
the stability analysis. 

Examples of Other First-Order Systems.—The follow- 
ing physical systems may be approximated by first- 
order differential equations, and they frequently occur 
in the synthesis of aircraft control systems: 


(1) Integrator: 


G(p) = A/p (IV-15) 


to an input voltage v; follows the equation dgp»/dt = Av; 


where A is a scale factor. Thus, the angular rotation 
“e 


¢o of the motor is proportional to v;.dt. This rela- 
0 
tion is represented by the block diagram of Fig. IV-7a in 
terms of the Laplace transforms %o(p) and V;(p). This 
transfer function represents a limiting case of G(p) = 
1/(Tip + 1)as7,\— ©. Thus, the transfer function is 
represented by a pole at the origin. The output re 
sponse to a sinusoidal input is obtained from 
tr/2 


Me” = G (iw) = A/iw = (A/w)e 


Thus 


M = A/w (IV-16) 


i.e., the output is attenuated at high frequencies propor- 


tional to 1/w and 


(IV-17) 


0 = —x/2 = —90° 


i.e., the output /ags the input by 90°. These relations 
are shown in Fig. IV-7c. 

(2) Differentiator: A rate gyro gives a voltage output 
vp proportional to the angular velocity d¢g,/dt of the pre- 
cession axis—i.e., % = A(dg;/dt) where A is the scale 
This case is, therefore, the inverse of the pre- 
ceding one. The block diagram is shown in Fig. IV-Sa. 
The transfer function G(p) = Ap has a zero at the origin 
(Fig. IV-Sb). The output response to a sinusoidal input 


is obtained from 


factor. 


Me” = G(iw) = Aiw = Awe'*’” (IV-18) 


Thus, 


(IV-19) 
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i.e., the output is proportional to the input frequency w 


and 


6 = r/2 = +90° (IV-20) 


i.e., the output /eads the input by 90°. These relations 
are shown in Fig. IV-Sc. 


(3) Simple lag network: 


G(p) = 1/(1 + Tip) (IV-21) 


The resistor and capacitor network shown in Fig. 
IV-5b has the transfer function 


Vo(p)/Vilp) = G(p) = 1/(1 + RCP) 


which is the same as that of the spring-damper system 
of Fig. 1V-4a with 7; = RC. This circuit is frequently 
used in order to introduce a phase lag into a system. 
If the input signal is sinusoidal with frequency w, the 
phase angle 6 may be read conveniently from Fig. IV-5c, 
since tan 6 = —w/'. The output response to a sinu- 
soidal input is given in Fig. IV-6. 
convenient to deal with the reciprocal of the delay time, 
and thus we define w, = 1/7) and correspondingly the 
transfer function for a simple lag network 


It is generally more 


G(p) = 1/[1 + (p/or)] = wi/(o1 + Pp) 
(4) Lead network: 
1 ew, + p . 
Gy) = -: » Ke IV-22 
( . p € (IV ) 











1?) a 
‘Ol 0.1 ae 1.0 10.0 
c) aw, 
Fic. IV-9. Lead network AG(iw) = [1 + t(w/ew,)]/{1 + 


(iw/w,)]. (a) Block diagram, V2(p)/Vi(p) = KG(p) = (ew: + p)+ 


(w, + p). (b) Representation in p-plane. G(p) has a simple 
zero at p = —ew, and a simple pole at p = —a;. (c) Response to 
sinusoidal input G(iw) = Mei®. M = magnification ratio (loga 


rithmic scale), 0 = 0, — 0. = phase angle, w = exciting frequency 
(logarithmic scale), and e = A = gain. 
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The resistor and capacitor network shown in Fig, 
IV-9 obeys the following relations: 
V2(p) : R, + RiR2Cp 
- en KG(p) = = — (V-33 
| i(P) (R + R:) + R,R2C Pp . 
or 
F 1 + (p/€o) 
G(p) = then 
] + (p @)) (IV-24 
K = Ro (Ri + Ro) 
where 
oO = (R, + R:) R,R2C IV me 
e = R/(Ri + R:) sae. 
The transfer function G(p) has a zero at p = — ew, anda 
pole at p = —w,. The value of € varies ordinarily be- 


tween 0.1 and 1. The characteristics of this system are 


as follows for a sinusoidal input: 
1 ew; + lw 


Me” = G(iw) = 
€ @ + lw 








where 
M = V [1 + (w?/e*%w,”)]/[1 + (w?/w,?)] (1V-26 
6 = tan"! w/ew, -ytan~! w/o, (IV-27 
—— 6. 





Gz = 

















(Cc) 

Fic. IV-10. 
[1 + 2(w/ew )]. 
E(w, + D)/ (ew, + p). 


+ i(w/o,)] + 
KG(p) = 
(b) Representation in p-plane; G(p) has a 


Restricted lag network G(iw) = [1 
(a) Block diagram, V.(p)/Vi(p) = 


simple pole at p = —ew, and a simple zero at p = —w);. (c) Re 
sponse to sinusoidal input G(iw) = Met? AJ = magnification 
factor (logarithmic scale), 6 = 6, — 0 phase angle, w = exciting 
frequency (logarithmic scale), and K = 1 = gain. 
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The angles 6; and @, are indicated in Fig. 1V-9. We note 
that for verv low frequencies (w > 0) 6 = 0, ./ = 1, and 
KM ~ e; for very high frequencies (wo > ©) @ = 0, 
Vx 1. and A.V = 1; and for intermediate frequen- 
cies 0 << 6 < 90° and e < AA <1. Thus, this circuit 
introduces a Phase lead into the system over a limited 
band of frequencies. The values of ./ and 6 are shown 
in Fig. 1V-9e. 
(5) Restricted lag network: 


G(p) = ear + p)/(ew + Pp), K = 1 (IV-28) 


The resistor and capacitor network shown in Fig. 


IV-10 has the transfer function 


I'(p 1 + RiCp 
Li = K.G(p) ~ — (IV-29) 
V(p 1+ (Ri + R)Cp 
G(p) = [1 + (p/w) ]/[1 + (p/ ear) 
K = 1 
where 
wo = 1/RC (1V-30) 
e = Ri/(Ri + R) (1V-31) 
The transfer function G(p) has a zero at p = —o, anda 
pole at p = —ew;. The value of ¢ varies ordinarily be- 


tween 0.1 and 1. 
The characteristics of this system are as follows for a 
sinusoidal input: 
1 + (w/w) 


Me” = G(iw) = ; 
1 + 1(w/ €w) 


where 
Mo= VEL + (w?/w7)]/[1 + (w?/e?a,?)] = (1V-32) 
6 = tan! w/w, — tan! w/ ew, (IV-33) 


6; A. 


The angles 6; and 6 are indicated in Fig. IV-10b. We 
note that this circuit introduces a phase lag over a 
limited band of frequencies. The values of V/ and @ are 
shown in Fig. 1V-10ce. 

(6) Simplified aerodynamic stability equations: The 
simplified equations for the rolling motion of an airplane 
are 

deg x 1, = bs 


dt? ” a (iV-s4) 


where / is the moment of inertia about the roll axis, Lg is 
the aerodynamic damping in roll, and ké is the torque 
applied due to the aileron deflection 6. Writing » = 


dg/dt, then 


I(dn/dt) + Le.n = ké 
In terms of Laplace transforms we find, if 7(0) = 0 


(Ip + Le).N(p) = k.D(p) 


STABILITY AND 


AUTOMATIC CONTROL 583 


Thus 

KG(p) = N(p)/D(p) = k/ Up + Le) 
or 
1/(U/Le).p + 1) (IV-35) 
K = &/Ls 


G(p) = 


The analysis of these simplified equations follows that of 
the simple lag network, Figs. IV-5 and IV-6. Some- 
times the aerodynamic damping Ly, is so small as to be 
negligible; then the analysis follows that of the simple 
integrator circuit, Fig. IV-7. 
C) Analysis of Second-Order Systems 

Let us consider the cantilever spring of Fig. IV-4 
which has a dashpot at its end B and include the finite 


mass m of this dashpot. Then, the additional force act- 


ing upward on the dashpot is —m.(d*yo, df*?). The 
equation of motion then beco:nes 
don Von 
m +c + kv = ky (I1V-36) 


dt* dt ! 
with boundary conditions 


yo(t) = w(O) att =O-—e 
dyop(t)/dt = y‘(0) att = 0 — « 


We introduce the following ratios 


wy = VR/m = natural frequency of undamped 
mass-spring system 
¢ = (c/m)/2w) = ratio of actual damping coefficient 


to critical damping coefficient 
Then the equation of motion becomes 


d Von : Yon ‘ je 
2Fa° Wo”. Von = Wo”. Vi 


dt® dt ‘ 
and its Laplace transform is 
(p? + 2a». p + wo”). Yoo(p) = wo”. Vi(p) 
Thus 


: wo" 
) oo(P) = < 


-. ¥i(p) + 
p° i 2fwop + a" 


Yo(p) 


p yy’ (0) + (p? + 2Fanp) yn(Q) 
p? + 2fwop + wy” 


V,(p) 


The term }(p) represents the Laplace transform of the 
output motion which is due to the input signal y,(/). 
The term Y,(p) represents the Laplace transform of the 
output motion which is due to the initial boundary con- 
ditions y,(0) and y,’(0).. The transfer function in this 


case 1S 
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b) RESPONSE TOA STEP-INPUT: 9j(t)=h (t) 
Fic. 1V-11. Second-order system characteristics 
Yo 9” 
K .G(p) = a) saa dae = 
Vi(p) p* + Woop + a? 
and thus 
K =] 
G(p) = wo” (p? + 2fwop + wo" (IV-38) 


This transfer function G(p) has no zeros, but it has two 
poles at p, and ps, respectively, where /; and p» are the 
roots of the denominator D(p) = p? + 2fwop + wo? = 0. 
These are 

pi = w(—-¢ + VE? — 1)/ 


; raxc> |] 
pe = wo( — F — Vv o — 1)\ 


hemnt + ert - = (IV-39) 
A + iv = ane’ 
ao <I 
pe = wo =< = iv | — c*) = 
A — iv = we 
We shall find that \ = —apf is the damping factor and 
v= V1 — §&°. wis the natural frequency of this second- 


The roots p; and pz: are complex conju- 
For a fixed value of a, the locus of these 


order system. 
gates if¢< 1. 
roots as the damping coefficient c is varied is on a circle 
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SCIENCES 


of radius w) about the origin, Fig. [V-11. This is easily 


seen, since \? + vy? = wy’, which is the equation of 4 
circle in the p-plane. For ¢ = 0, the roots are p, = ij 
and p2 = —iw. For ¢ = 1, the roots are p, py = 


—wy. This point corresponds to critical damping of the 


system. For intermediate values 0 < ¢ 1, the roots 
lie at intermediate points on the upper and lower half oj 
the circle, respectively. As ¢ increases above ¢ L, 


the root /; moves along the negative real axis toward the 
origin, and the root p2 moves along the negative real 


axis toward \ = —o., 

The solution y,(/) due to the initial boundary condi 
tions is found as follows for the case ¢ < 1. We intro- 
duce AX = —fey and vy = V1 — £&?.winto V,(p). Since 

D(p) = p® + 2aop + wo? = (p — A)? +0 
we have 
P. yn’ (0) p.(p — 2d). ¥(0) 
Y',(p) a Ae ae 
(Dp = A)? + ¥ (p — A)* + 
p.v yp’ (0) P(p — Xd) 
— = = + 3 
OHA a y (p — rX)* + 


p.v —X 
y,(O) + s te v»(O) 
(p — X)* + Pv’ v ; 


and thus from Eqs. (B-6) and (B-7), 


Vo "¢ 0) 
Vol {) = € 
Vv 


M 4 Al 
sin vt + y,(0).e”. cos vt + 


am ‘ 
yr(O).e”. sin vt (1V-40) 
y 


Since \ is a negative number, it is seen that the initial 
motion due to the boundary conditions is a damped sine 
wave with damping factor \ = —fwo and with the 
natural frequency vy = V1 — (a. The solution for 
¢ > 1 is obtained similarly. 

Let us examine more carefully the forced motion due 
to the input y,(f). 

(1) Special case: y(t) = h(t). 
of a unit step input in displacement, y;(¢) = A(t) and 


For the special case 


V(p) = 1. The corresponding solution for the second- 
order equation is then 

wy” wo" e 
) ol P) = = (IV-41 


Pp? + 2Canp + wo" (p — ))° + y* 


and thus from Eqs. (B-6) and (B-7) we find 


y(t) = 1 — [eos vt — (A/v). sin vt].e” (IV-42) 


This solution is plotted in Fig. IV-11b. We note that 
the output y(t) approaches the input y,(/) for large 
values of ¢. After the initial application of the force, the 
output overshoots the equilibrium position if ¢ < 1; 
it approaches the equilibrium exponentially if ¢ > 1. 
The region 0.4 < ¢ < 1 is ordinarily considered a good 
compromise between quick approach to the input and 


excessive overshooting. 





(2 
sinu 
ven 


and 


and 


and 


TT 
tran 
catic 
ane 
sent: 
was 


has « 


wher 


and 


At 
valu 
note 
phas 
shov 


scale 


of the 
’ TOOts 


half of 


mm 


| 


rd the 


e real 


‘ondi- 
intro- 


Since 


V,(0) 


-40)) 


itial 
sine 
the 


for 
due 
‘ase 


ind 
nd- 














AERODYNAMIC STABILITY AND AUTOMATIC CONTROL 585 
(2) Special case: y;(t) = Vim.e'.-—The special case of N(p) = a»? 
sinusoidal input y,(t) = ym.e is solved most con- : 
veniently by applying Eq. (A-15). We have Db) = 2p + un = 2(p — d) 
Vi(p) = Yin. [p/(p — iw)] pf = AX +4, D'(pi:) = 2iv 
and : 
‘ o = — wy ) o) = —Vryp 
N(p) wo Leas = 
G(p) = = — - Ss a , 
D(p) pp? + wy P+ wo?) (TV-43) Thus, 
G(p) = wo” [(p =Ppuy — p2)] | Yo(p) = Vim {[p/(p ie w)|.G(p)} (IV-44) 
and 
a “ N(pi) N(po) 
yo(t) = Yim | G(iw) .e" + = : a _ ae 
(pi — tw).D’(p1) (p2 — tw).D" (peo) 
, |G ) mT rt | wo I | (A+iv)t i‘ wy” A—iv)t 
= Vim | G(tw).e Yim- —— ie ii ee 
A +2(v—w) 2 A—t(v+tw) —2y 
A B 
and simplifying, we find 
vo(t) = y | i —— et] +y Ae“e'™ + Be™.e~™ (IV-45) 
° (wo? —- w*) oo 21f wow nia , _ 


[yo(t) lematy 


The second part of this solution, [yo(t) |transient» 18 the 
transient disturbance associated with the initial appli- 
It is damped out, since \ = —f{wp is 
a negative number. The first term, [-o(t) ]steaay, Tepre- 
sents the steady-state vibration, at the same frequency 
was the input, which results after the starting transient 
We again have the relation 


cation of the force. 


has damped out. 


[vo(t) lsteaay/Vi(t) = Me” = G(iw) 


where 


w”) + 21Fwyw 


Me™ = w?/[(wo” — 


Thus, for the second-order system 


(w? a”) |? + [2¢(w/wo) |? 
(IV-46) 


M = |G(iw) 


Vii — 


and 


(w/ay)?]  (IV-47) 


tan 6 = —2¢(w/ay)/[1 
At w = w, we find M = 1/[2¢.(w/wo)). The maximum 
value of M occurs for a slightly lower value of w, as 
noted in Fig. 1V-12. Also we notice that again there is a 
phase lag. The magnification ratio and phase angle are 
shown in Figs. IV-12 and IV-13 to linear and logarithmic 


scales, respectively. We note that at w = w) the phase 


[vo(t) aw 


angle is —90°, and, as w/w) > ©, the phase angle ap- 
proaches — 180°. 

The logarithmic plot of / indicates that for large 
values of w/w, we have VM ~ 1/(w/w)*, and thus 
logi M =~ — 2 logy (w/w); i.e., the curves approach a 
slope of —2. The acoustic engineer would call this a 
slope of —12 db. per octave. 

The transfer function G(iw) = Me" is easily plotted 
on the basis of the previous calculations. The resulting 
curves are shown in Fig. IV-14 for various values of ¢. 
The vector G(iw) starts from M +1, 0 0 when 
w = 0; the vector increases to a maximum at resonance 
and then decreases in magnitude again as \J — 0 and 
6 = —180° asw> >. 

(3) Other Second-Order Systems.—The 
physical systems that occur in the automatic control of 
aircraft may be approximated by second-order systems. 


following 


(a) Mechanical spring-mass-damping systems 


(1) Hydraulic servo system.—-A measured re 
sponse curve for such a system is shown in Fig. 
IV-15. 

(2) Rate gyro. 
previously used equation AG(/p) 


A better approximation than our 
Ap is 


A.p 


(p a)? + 2¢(p wo) + 


= (IV-48 


KG(p) 
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Fic. 1V-12. Second-order system: Response to sinusoidal input 
G(iw) = Me? (linear scales). 
(3) Accelerometer 
Ap* 


KG = 
(?) (p/w)? + 2¢(p wo) + 1 


(IV-49) 


(4) Electric motor 


KG(p) = A/{p.[p + 1/7) 


: (1V-50) 


(b) Electrical networks 
Inductance, capacitance, and resistance networks 
(c) Aerodynamic systems 
(1) The short pitching oscillation of an airplane 
(2) The phugoid pitching oscillation 
(3) The yawing oscillation 


The above-mentioned aerodynamic systems can only 
approximately be represented by second-order systems; 
the actual motions are described by the equations of 
The solution of these equations is well 
A typical solution of 


Appendix (C). 
known to the aerodynamicist. 
these equations for a modern plane is given in Fig. IV-16 
which shows the dynamic longitudinal response for an 
F-80 airplane as computed by Walkowicz.*! This cal- 
culation, incidentally, also illustrates the great im 
portance of aeroelastic deformations. 

(4) Analogy Between Transfer Functions and Hydro- 
dynamic Potential Functions.—There is a relation be- 
tween the elementary transfer functions G(p) and 
corresponding complex potential functions w = ¢ + 1p 
expressible byw = log G(p). This analogy is summarized 
in Figs. IV-17 and IV-18 for some simple examples. 


SCIENCES-—SEPTEMBER, 1951 


We write p = re’ and G(p) = R.e’’. We note that a 
hydrodynamic source corresponds to a zero of G(p 
(Fig. IV-i7a), a hydrodynamic sink to a pole of G(p 
(Fig. IV-17b), and that the stream lines Y = constant 
correspond to lines of equal phase angle 0 for the vector 
G(p). A simple lead or lag network corresponds to 4 
source and a sink of equal strength on the x-axis (Fig. 
IV-18a). A simple second-order system corresponds to a 
pair of sources at conjugate complex points (Fig. IV. 
ISb). 


(D) Composition of Elements 


(1) Series Circuit.—Consider a series arrangement of 
nelements. The rth element has a gain A, and a trans. 
fer function GAp) = M,.e”. 
equivalent to a single block having a transfer function 


K .G(p) where 


This series array is 


a? on Co Coe ee 
G(p) = (Me™). (Mee™).. . (Me)... .(M,e*") 
(M,.M....M,...M,).e%T*" ee 


II 


(1V-51 


Such an arrangement is illustrated in Fig. IV-19a for 
three elements in series. The Laplace transform of the 
input to this system is },(p), the Laplace transform of 


the output due to the input signal is 14(p), and thus, 
Vs(p)/Vil(p) = K.G(p) 


where 


K = ky : Ko. K; 
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Second-order system: Response to sinusoidal input 


Fic. IV-138 
Ve'? (.\/ and w to logarithmic scales) 
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. - ° ’ Ta 

The overall transfer function G(p) Me™ may be ob 
tained graphically by adding the logarithmic response 
curves, since 


log Me" = logy [( Mie") . (Mee™) . (Msc) 


and thus, 


logio .V/ logio M, + logyw Me + logy M3 (IV-52) 


and 


6=0+0+4+ 0; (1V-53) 


This is illustrated in Fig. IV-19b for the series arrange 
It will be noted that it is 
possible to approximate the product by adding the 


ment of linear networks. 


straight-line approximations for each of the elements. 
(2) Parallel Circuit. 


in a similar manner in Fig. IV-20a. 


A parallel connection is shown 
The outputs add 
following the same law as electric currents. 

(3) Simple Feedback Servo. 
is shown schematically in Fig. IV-20b. 


A simple feedback servo 
The equations 


that apply are 


E(p) = Y,(p) — Yol(p) (1V-54) 
Vo(p) = K.G(p).E(p) (IV-55) 
and thus, 
E(p)/Yi(p) = 1/[1 + AG(p) (IV-56) 
and 
Vo(p)/Vi(p) = KG(p)/[1 + AG(p)] (1V-57) 


This overall ratio ¥o(p)/ Yi(p) could be represented by a 
system transfer function A,.G,(p). 
The feedback servo is the most important element for 


automatic control systems. The input }j(p) is com- 
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Vi(p) 


plane as a control signal. 


ference E(p) = },(p) is passed on to the air 
Chis trick permits the output 
te he controlled by an error signal rather than by direct 
input, and, as a result, the accuracy and speed of re- 
sponse are much increased. It will be remembered that 
this technique is the same as that used by Siemens in 
1843 in his ‘“‘chronometric governor” (Fig. III-1). 

(4) General Feedback Servo.—The general feedback 
servo is shown in Fig. IV-20. Its arrangement is similar 
to the simple servo, except that the signal that is fed 
back is modified in amplitude and phase by a corrective 
network with the transfer function K»y.G2(p). We have 


the relations 


E(p) = Yilp) — Yy(p) (IV-58) 
Vi(p) = KxG2(p). Yo(p) (IV-59) 
Yo(p) = KiGi(p).E(p) (IV-60) 
and thus, 
Yo(p) K,G,(p) 
. - ire pees (IV-61) 
Y;(p) 1 + [KyGi(p)|.[KeG2(p) 


(E) Stability of Feedback Networks 


(1) Routh Criterion for Coefficients of the Frequency 


Equation.—We have seen that for the simple feedback 
network the overall system transfer function is defined 
by 

Yo(p) ; N,(p) KG(p) 

: I = K,.G,(p) = f = a\P} (IV-62) 

Yi(p) D.(p) 1+ KG(p) 


The function G(p) will ordinarily be expressible as the 
ratio of two polynomials, each of which is given in 
factored form. Thus, for example, the roll-control sys- 


tem shown in Fig. IV-2 has the gain and transfer func- 


pared with the actual output }o(p), and only the dif- — tion: 
ple, N(p) j ki/Le ) sl + elope + Tsp « + spl 5 pay” } 
KG(p) = ee ast me oe ape ts Sree = (IV-63) 
Dip) lp. + pT, 6 414+ Top 14+ Tsp 14+ Tip Up? + 2eop + wo?f 


due to airplane 


and computer 


The numerator N(p) of this function is of the third de- 
gree in p, the denominator D(/p) is of the seventh degree 
in p. 


the three zeros and seven poles of G(p). 


Each is given in factored form so that we know 
We note that 
N,(p) 


K,G,(p) = D,(p) = 


N(p) 


16 
D(p) + N(p) (IV-64) 


Thus, the zeros of G,(p) are the same as those of G(p). 
The poles of G,(p) are, however, changed. The degree 


of D.(p) is the same as that of D(p), since ordinarily 


due to amplifier 


due to aileron 
control servo 


D(p) is of equal or higher degree than N(p). Let us 


write, in general, 

(p — 2).(p — 2). . (Pp — Sm) 

- pi).(p — po)...(p — pnd 
(IV -65) 


N.(p) 


Kitt) = = 
s(P) D,(p) (p 


where 2}, 22, . . . 3m are the zeros of G,(p), pi, po, ... Pn are 
the poles of G,(p), and A is a constant. 
The principal mathematical problem of the servo 


analysis is ordinarily to go from the function G(p) to 
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Second-order system: Vector plot of response to 
sinusoidal input G(ijw) = Ae'?®. 


Fic. IV-14. 


the function G,(p) in factored form. For after the latter 
has been written in factored form, the solution of the 


problem is immediately available. If, for example, the 


; ee “ N3(pr) ; 
yolt) _ Faas | KiGl i + > . f , ef" 
r=1 (pr — tw).Ds'(pr) 


steady-state vibration 


where the p,’s are the roots of the so-called frequency 
equation D,(p) = 0. 
order for the terms e’” to vanish as t > ©, it is neces- 
sary that the real part of each of the p,’s be negative. 
This means then that all of the roots p, must lie in the 
negative half of the p-plane. The determination of the 
roots p, is ordinarily not simple if the equation D,(p) isa 
In the above example of the 


It is obvious, therefore, that in 


high-order polynomial. 
roll-control system, it would be necessary to solve a 
seventh-order equation. This problem is generally 
analogous to the problem of determining the dynamic 
stability characteristics of an airplane. In the past, the 
aerodynamicist studied this problem, say, for a fourth- 
degree equation, by expanding D,(p) in the form 


Dp) = p'+ Bp? + Cp? + Dp+E 
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Fic. IV-15. Frequency response of a hydraulic servo: G(iw) = 
Me'? (experimentally determined ). 
input is a sinusoidal force F; = Fra, .e'’, then 
Vi(p) = Frnax.-P/(p — tw) 
and the solution from Eq. (A-15) is 
(1V-66) 


initial transient 


and then applied Routh’s criterion*! that for stability 
the coefficients B, C, D, and - and X should be all 
positive where Y = BCD — D*? — EB?. This technique 
has the disadvantage that one loses all insight into the 


physical behavior as the various parameters are 
changed. The following techniques are superior in this 
respect. 


(2) Root Method Method).—-Let us 


examine the problem of finding the poles of G,(p) accord 


Locus (Evans's 


ing to a method recently introduced by Evans. 
He considers the function 


KG(p) 
1 + KG(p) 


_ Ns(p) 


K,G.(p) = 
t D.(p) 





and studies the variation of the poles of G,(p)—1.e., the 
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Fic. 1V-16. Dynamic longitudinal response of an F-80 air 


plane.*!' M = 0.6, h = 10,000 ft. for (b), (c), and (d). = 
rigid airplane, - = elastic airplane. 
roots of the frequency equation D,(p) = O as the 
parameter A is varied from AK = 0 to A = ©. Let us 
assume that G(p) is given in factored form 
; (oe — SoD — Se) 2.0.5 {0 = Za) — 
G(p) = A (1V-67) 
LP — Gye — Oe) ss (OP — @) 
where 


A= (q1-q2 Qn) / (S81. S2- Sun} 


Ordinarily the denominator of G(p) is of equal or 


higher order than the numerator—i.e., n > m. Let us 
express each of the factors in vector form. 

p aa = ke Ie gag 

p —- Ss = Se.e™ 

pb — Sm = Se 

p—-a = Qe" (IV-68) 

b-@ =Q ef” 

P— dr = Qe” J 
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Fig. IV-2la shows that the vector p — s, = S,e’* 
starts from the end point of the vector s, and goes to the 
end point of the vector p. The length of the vector 
pb — s,is S,; its angle with respect to the positive axis is 
g,. Similarly, the vector Q,.e"” goes from the end point 
of the vector g, to the end point of the vector p. The 
points s, and g, are, of course, the zeros and poles of 


G(p). The point p is a variable point in the p-plane. 
Thus, 
; (SPO ASE) «ilo 
G(p) - A 10 10e 10, = 
(Oy.€°°) (Cet)... He") 
S1S9 — be [ —— ee o (0, 4-004 +6n)] yr 
fortet...tend—m 1 (TV-69) 
Q1.Qo... Qn 


The function G(p) may therefore be interpreted as 
another vector of magnitude R and phase angle 0 
where 

G(p) = Re 
and 


R = A(S,Se eee be (QiQ2 a Qn) 

















6 = (¢g1 t go t.--t Om) — (A + 2 + (IV-70) 
ee of 
- * © 
ZERO AT ORIGIN (a) SOURCE AT ORIGIN 
Gr)= 2 w (2) = Qt iy = LOG 
R o-<-r $=LOGr p 
& +8 w =106 Ge) i 
® =LOGR 
¥ = 
Peat 
POLE AT ORIGIN — (6) SINK AT ORIGIN 
| ows 
G(p)= wi(p)= Qtiy= -LOG p 
R -« di 9 =-LOGr 
r y =-@ 
® =-8 w=LOG Gip) 
9 =LOG R 
y:e 
Fic. IV-17. Analogy between transfer functions G(p) and hydro 


dynamic potential functions. 
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/ TABLE IV-1 
\ n—m Asymptote 
l 180° 
SINK 2 +90°, —90 
3 +60°, —60°, 180° 
\ 4 +45%, —45°, , +135°, —135 
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Analogy between transfer functions G(p) and hydro 


dynamic potential functions. 


Fic. 1V-18 


The poles of G,(p) are found from the condition 
D.(p) = 0-—that is, 

1+ KG(p) = 0 
or 


1+ KRe® = 0 
which may be written 
KRe? = —1 


or 


, 0 
KRe™ = 1.e" 


This vector equation may be resolved into two scalar 
conditions: 
(1) Magnitude condition: 


KR=1lorK =1/R (IV-7 la) 


Thus, the vector AG(p) will have a magnitude of one. 
(2) Angle condition: 


0 = ror 8 = 180° (IV-71b) 


Thus, the vector AG(p) will have an angle of 180°. 
Both of these conditions must be satisfied if the point p 
is a pole of G,(p). (The angle condition is also satis- 
fied by 8 = mw + integer X 27.) 

Evans studies this problem in two steps. He deter- 
mines first of all the locus of all points which satisfy the 
angle condition (2). Then he computes the values of K 
for various points along this locus by measuring the 
lengths Si, So, ...S, and Qi, Qo, ...Q, and computing 
the ratio R. Fora given value of K, the points ~i, po, ... 
pn along these loci are then the roots of D,(p) = 0—i.e., 
the poles of G,(p). 


Evans has developed a number of simple rules that 
determine the locus of points satisfying the angle con- 
dition (2). These may be verified by simple examples: 

Rule 1. 0, the poles of G,(p) coincide 
with the poles of G(p). 
the poles of G(p). This can be easily verified. 

Consider, for example, 


For a gain K = 
Thus, the root locus starts from 


G(p) = [(p — 51)/(b — Qi) )-qi’ss 
where s; and g; are on the real axis. Then, 
D(p) = 1+ Kl{(p — s)/(p — @))-G/ 1 = 0 


This is only possible if [(p/g.1 — 1)/(p/s1 — 1)] =—K. 


As K — 0 itis necessary that p — qi: ~ O—1-e., that 
p — qi, which is the pole of G(p). 
Rule 2.—For a gain K = ©, the poles of G,(p) coin- 


cide with the zeros of G(p) or they are at p = ~. If 
the number of poles 1 exceeds the number of zeros m, 
the direction of the infinite asymptotes of the loci are 
shown in Table IV-1. Thus, the root locus ends on the 
zeros of G(p) or on the above-mentioned infinite asymp- 
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Fic. IV-19. Series circuit. 








ol 


th 


po 


alt 
loc 
th 


po 
pa 
spt 
thi 
all 
all 
if t 
of 


the 
fro 
fro 
Zer 


effe 


or 


S that 
eC con- 
nples: 
incide 
‘ from 


0 


> —K, 
_ that 


coin- 

li 
Os mM, 
Ci are 
n the 
ymp- 


NTS 


IN OF 

















AERODYNAMIC 


This can be verified, for example, if G(p) = 
Then 


totes. 


((p — )/(P — 1) )-(Q1/ 51). 


Dip) = 1+ K\(p — 51)/(b — q))-(q/s) = 0 
only if 
(q/si (Pp — 1)/(P — m))] = —(U/K) 


or |/K — 0, it is necessary that p — s;—> 
If the function 


Thus, as A > 
0, or Pp > Si, which is the zero of G(p). 


G(p) = (1/p) [(p — 31) /(p — q) )-(Qi/s1) 


then 


Dip) = 1+ K(p — 51)/[P-(b — q)I-@/s1) = 0 


(a sU/p)l(p — s)/(p — qa)] = — (1/K) 


As K > o, we find that either p > s; or p > © are 
possible solutions. 

Rule 3.--The root locus along the real axis is along 
alternate segments connecting the zeros and _ poles 
located on the real axis, starting with the one farthest to 
the right. 

This rule may be easily verified by considering any 
point on the real axis. The angle to this point from a 
pair of conjugate complex poles is +@ and —8@, re- 
spectively, and thus their sum is zero. The angle to 
this point from the poles or zeros on the axis is @ = 0 for 
all poles or zeros to the left of p, and it is @ = 180° for 
all poles or zeros to the right of p. Thus, the sum is 180° 
if there is an odd number of poles and zeros to the right 
of p. 

Rule +.—If there is a breakaway of the root locus from 
the real axis, the point of breakaway may be estimated 
from the conditions that, for a small displacement Aw 
from the axis, the increase in angle due to the poles and 
zeros on the axis to the left must be just balanced by the 
effect of those to the right. 

Example: 


G(p) = 1/[p.(p — qd(p — @)] 
We note from Fig. IV-24 that the point of breakaway 
I’, is given by 
Aw/Xp = (Aw/x1) + (Aw/Xe) 
or 


1/X%o = (1/1) + (1/2) 


Rule 5. 
the imaginary axis into the positive half plane can often 
be estimated by taking advantage of the properties of 


The point U at which the root locus crosses 


the right angle. 
Example: 


G(p) = 1/|p.(p — u)(p — @)) 


At the point u the angle 6, from p = 0 is 90°, the angles 


6 and @; must therefore also add to 90°. Take a trial 
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(©) SIMPLE FEEDBACK SERVO 
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(c) GENERAL FEEDBACK SERVO 


Fic. 1V-20. Block diagrams of various circuits, 


point U’ and draw a 45° line. Then for u to be on the 
root locus, it would be necessary that a = £6 in Fig. 
IV-24. 

Rule 6._-The direction of locus departure from a pole 
(or locus approach to a zero) may be easily estimated by 
computing the angle at the pole under consideration 
from all of the other poles and zeros in the field. 

For example, in Fig. IV-26 we wish to determine the 
angle at which the locus departs from the pole gy. At 


the pole qs, we have 


6, = 150° ¢1 = 120° 
6 = 90° g = 45° 
6, = 90° 


6, + 0 + 0; = 330° g1 + go = 165° 


Now let us consider a point P on the locus a small dis- 
tance away from the point g;. This point P must be at 
an angle 6; from qs such that 

(gy + ¢2) — | (A, + Ao + 6) 4 A, ) = 180° 


ee 


165° 330° 


or 
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VECTOR p = OP VECTOR p = OP 
VECTOR S,= OA VECTOR G = OC 
VECTOR # ~5y* Spa" AP’ VECTOR -F Q,< CP 


LENGTH OF VECTOR P —Sr=Sr= AP LENGTH OF VECTORS-J, =Qy =CP 
ANGLE OF VECTOR ® —Sr+@r= ANGLE DAP ANGLE OF VECTOR? -Jr =@, = ANGLE ECP 


(a) VECTOR REPRESENTATION OF 2-5 (6) VECTOR REPRESENTATION OF # -gr 


Ke Q @ 2,0, 
S, S, 
8 = (Qt O,)— (6,+ 0, -6, +6,) 


TRIAL POINT 





rf 


4 
(c) ILLUSTRATION OF ROOT LOCUS CONSTRUCTION FOR Gip)= 


(-s, (ep - Se) 
(pq) (P-gep-gsXP-ga) 
Elements of root-locus construction 


Fic. IV-21. 








K=@ K=2.0 K=O 


» P q x 


= + 








PS, Ha, 
Pp 





Fic. IV-22. Root locus on axis: G(p) = [(p — s1)/(p — qi)]> 
gi/S,;. Locus moves from a pole (for A = 0) toa zero (for K = 
o), Example: g, = —1,5 = —4,p = —2; |p—-—H] =Q = 
1; |p —5;| = S, = 2; therefore, K si/qi: Cp — qi)/(p — 51) 
== ©) 


i.e., the point p is at an angle 6; = +15° with respect to 
qs. 

For intermediate locations, the root locus is found by 
taking a number of trial points, measuring the angles 
6, and ¢,, and taking the algebraic sum a = 
This addition can be conveniently performed by a 
special plastic device called the “‘spirule’ (Fig. IV-27). 
A spiral curve on the device gives the logarithm of a 


2¢, — 26,. 


length in the form of an angle so that addition of angles 
gives the logarithm of the product 


K = AWOy. . .Qe)/(Si.Sr.. . S) 
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Fig. IV-28 shows the root locus for the roll stabiliza- 
tion system of Fig. IV-2. A gain A = 670 would give 
satisfactory stabilization in this case with the roots 
indicated in the figure. The values of the gain corre- 
sponding to various points along the locus are shown on 
Fig. [V-28. 

The above described technique is instructive in that it 
suggests methods for correcting adverse trends in the 
directions of the loci. For example, Fig. IV-29a shows 
the root locus for a simple system. 
the effect of a lead network in stabilizing this system up 


Fig. IV-29b shows 


to larger values of the gain. 

The above technique also has a great advantage in 
visualizing clearly the trend of the roots in conditionally 
stable systems. For example, Fig. I[V-30a illustrates 
the problem of an autopilot stabilizing an airplane that 
is statically unstable (without the autopilot). The 
locus starts at low values of gain in the right half plane. 
It crosses the axis at a gain Ay, and the system becomes 
stable. At a still higher value of gain A», another 
branch of the locus moves into the right half plane, and 
the system becomes unstable again. Fig. [V-30b shows 
the root locus for an autopilot controlling a statically 
stable air frame. It will be noted that both of these 
cases lead to a stable airplane-plus-autopilot combina- 
tion. This example, taken from a systematic investiga- 
tion of a large number of systems by Martin, illustrates 
his conclusion that the autopilot designer would most 
prefer to have an airplane possess small inherent “‘static 
stability’ or even a small degree of instability. 

Another interesting proposal suggested by Evans is 
to use an electric analogy to determine the root locus. 
This method is still under development; thus far, the 
graphical technique has proved to be the most con- 
venient and rapid procedure. 

High-Order 


criterion for 


(3) Routh’s Criterion for Stability of 
Systems.—Routh deduced 
stability of linear vibrating systems by using a theorem 


his famous 
from the theory of complex variables (due to Cauchy). 
This theorem states that if f(p) is the vector represent 
ing the function of a complex variable p, this function 


* within a closed contour 


{(p) having s poles and r zeros 
C of the p-plane, and as the point p moves around the 
closed contour C encircling it once in the clockwise direc- 
tion, the vector /(~) carries out r-s clockwise revolutions 
about the origin. Routh applied this theorem to find 
the number of roots of the frequency equation D,(p) = 0 
in the right half of the p-plane. He took as boundaries 
for the contour C: (a) the imaginary axis, p = /w |along 
this axis D,(p) = D,(iw)|; and (b) the infinite half circle 
to the right, p = Re”, where R— ©. Nyquist carried 
over the application of these criteria to the simple feed- 
back network for which 
K,G;(p) = N.(0) = — 
Dp) 1+ KG(p) 


* A pole or zero of order m is counted m-fold. 
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Fic. 1V-25. Infinite asymptotes: 


(1)-(3 


») (7 





G(b) = tn = — 
Wp (p + 1).Cp + 3). FH 9).(19 F 
= : © = 
K=@ At K = 0, locus starts from poles at —1, —3, —45, and 7 
Fic. [1V-23. Root locus moves from poles to zeros, or to infinity Infinite asymptotes at + 15° and +135", since for large p, G(p) 
1/p*. Note breakaway from \-axis at 4, = —1.75, where 
, o + 5 (13/5) A A A A 
G(p) ; - - = ) w v ) ' 
Nb p+1 [p+ (34+ 2) lp + (8 — 2i —= + - + —- 4 1, =0 
; Mi +7 i+ 5 1 + 3 i + 1 
Locus moves from poles (for A = 0) to zeros or to p (for 
" = oo), Infinite asymptotes at +90°, and —90°, since for 





large p, G(p) 1/p? natant a ‘iceland 





,iw ASYMPTOTE 
v7 














K= 8680 
- x 
- 420 
9, 
- / 02290°\ 6, =150° 
00| $2=45° $, =120° . 
hte Umar! a 
S2 G2 3; 91 
\K= 8680 @*+90° 
\ ie @3 
Fic. [V-24 Breakaway from axis and crossover to unstable q3 “a 
region 
; (0.001) (2) (6) 
G(p) = a 
(p + 0.001) (p + 2) (p + 6) 
At K 0, locus starts from —6, —2, and —0.001 Infinite 
asymptotes at +60° and 180°, since for large p, G(p 1/p’. 
Note breakaway from \-axis (point V) at 4; = —0.9, where 
Aw Aw Aw 
+ 4 0 
i + 6 A + 2 A, + 0.001 
Note crossing of w-axis (point U) where a = 8B. Fic, 1V-26. Direction of locus departure from pole. 
or The point [D;.(p)]/[Ns(p)] = 0 corresponds to | 
[K.G(p)] = —1 or 1/[G(p)| = —K. Thus, if the 
; D.(p) | vector |D,(p)|/|N.(p) | encircles the origin as the point p 


= = | i ; 
K,.G.(p)  Ns(p) KG(p) moves clockwise around the contour C of the right halt 
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p-plane, then the vector 1/[AG(p)] encircles the point 
—1, or the vector 1/{G(p)] encircles the point —A. 

Let us illustrate the application of this method for the 
simple transfer function 


G(p) = 1/[p.(1 + pT) + p72) ] 
and thus 
1/G(p) = p.(1 + pl).C1 + pT) 
and 
_ Db) _ p.A + pT). + pT») 7 
K,.G,(p) 7 N;(p) i K 
(a) Along the imaginary axis: Pp = Iw 
—— = tw.(1 + ww7)). (1 +1072) 
G(1) —_—_—— —_— — 
M,e™' Moe 
\t ao => 0 
1/(|G(iw)|] = +0.7 
Atw= +a 
| 


= lim (7Z).(1 + 1L7)).(1 + 1L72) > —i@ 


L—>« 


G( lw ) 


The vector Mye and Mee which multiply iw turn 
this vector progressively through +1S0°—1.e., from 
+90° through + 180° to +270°. The end points of the 
vector |/|G(iw)| are represented by the curve OAB. As 
p traverses the portion of the imaginary axis from 
w = 0tow = — ~, theend points of the vector 1 /(G(iw) | 
follow the curve OA’B’, which is the mirror image of 
OAB in the real axis. 


i ? ~ 
2 " io FG 8 7 6 
il > i 
A2 ~~ 
SPIRULE e 
“As Universal Equipment Company ~~ 
AIRCRAFT DIVISION 
7 
af CULVER CITY. CALIFORNIA 
“a COPYRIGHT 1950 NORTH AMERICAN AVIATION. ine 


Fic. IV-27. 


AERONAUTICAL 
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(b) Along the infinite contour: p = Re’ as R > « 


= (Re). (1 + Re'*7T,).(1 + Re’®T2) > 


7(p) 
R?.e"* Ty7, 


As the vector p rotates 180° clockwise along the infinite 
circle from B to C to B’, the vector 1/[G(p)] rotates 
3 X 180° = 540° clockwise along the infinite circle. 

Let us examine now whether the vector 1/[G(~) | en- 
circles the point —A as p describes the contour C. 
Imagine a string connected to a pin through the point 
—kK,. As the end of the string traverses the contour 
C’ in the plane of 1/[G(p)], we note that the string has 
i.e., the 
positive and negative rotations are canceled out. We 
thus deduce that the system is stable for the gain A, 
i.e., it has no roots with a positive real part in the right 
half of the p-plane. 

Let us now examine similarly whether the system is 


carried out no net rotation during this process 


stable for the gain A»,—1.e., assume the pin anchored at 
the point —A, and trace the number of revolutions 
We find that the string 
carries out two clockwise revolutions about the point 
—K». We thus deduce that the gain A, makes the sys 
tem unstable and that there are two roots with a positive 
real part in the right half of the p-plane. 

It is obvious by inspection that the system becomes 


carried out by the string. 


unstable as the point —A passes through the point A. 
This criterion for stability is called Nyquist’s criterion; 
the diagram (Fig. IV-31d) is called Nyquist’s diagram 
for the direct transfer function G(p); and the diagram 
(Fig. IV-31c) is called Nyquist’s diagram for the inverse 
transfer function 1/[G(p)]. We do not, however, obtain 
any quantitative estimate from this diagram to indicate 
how unstable the system is. The root locus diagram 
(Fig. IV-29a) gives the solution of this same problem 
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Root locus 


Fic. IV-28 


including a quantitative measure of the negative damp- 
ing coefficient as a function of the gain A. 

In order to obtain a quantitative measure of the 
damping coefficient, positive or negative, it is in princi- 
ple possible to carry out the conformal mapping not 
only of the line p = 7w, but also of neighboring lines of 
\ + iw, keeping A = constant, 
keeping w = con- 


constant damping p = 
and of lines of constant frequency, 
Such a map reproduced from Hanny 
This process is ordinarily too tedious for 
Osborne” showed that for typical 


stant. 'S is shown 
in Fig. [V-32. 
practical application. 
control systems an approximate measure of the damp- 


ing coefficient ¢ for the most critical root is given by 
(IV-72) 


¢ = (1/60)(a/m) 


‘phase margin”’ in degrees at the cutoff fre- 
1 and m = slope of the curve 


where a = 
quency where AG(iw) = 
of logy AG(iw), vs. logy w at the point A.G(iw) = 1. 
Thus if a = 30° and m = 1.7, then ¢ = (1/60) X 
(30/1.7) = 0.3 
Another word of caution is necessary when the func- 


tion G(p) has r zeros in the right half of the p-plane. Then 


the vector 1/[G(p)] will carry out r counterclockwise 


rotations as p carries out one clockwise rotation around 


(c) 


plots of roll-control system 


the contour C due to these r zeros. If the string whose 
end is anchored at the point —A carries out r counter 
clockwise rotations in traversing the contour C, the 
system is stable. If it carries out r—s counterclockwise 
rotations, then the system is unstable and has s roots 
within the right half of the p-plane. 

It is seen from Fig. IV-31b that the function K 
[AK,G(p)| = {1/[G(p)|} + A and thus the system trans- 
fer function 1/{G,(p)] is obtained simply from 1/[G(p) | 
by a shift of the origin. This is the principal advantage 
of plotting 1/{[G(p)] instead of the direct function G(/). 

Bode presents the function KG(iw) = KMe” from 
w = 0tow = by plotting its magnitude AM and 
phase @ separately vs. w(K .M and w on a logarithmic 
scale). This is generally a more convenient method of 
presenting the data than Nyquist’s vector diagram. 
The transition from stability to instability is again in 
dicated when G(iw) = —1/K—1.e., when KM = | 
(i.e., logy A.W = 0) and @ = 180°. Fig. IV-33 presents 
the Bode curves for the same examples as Fig. IV-29a. 
In this example @ = 180° atw = 3.4. At this frequency 
M = 0.074. Thus, the transition from stability to in 
stability occurs at 1/AK = 0.074 or K = 13.5. A de- 
sirable design condition in terms of the Bode diagram 1s 








596 JOURNAL OF THE 








bz % 


(2) Si) RRP 





a a vi] b 
%2 %e Se | 


I _P- Se 
P(P-4.)(P-4.) P-4¢ 





(b) G(p) = 





(a) Root locus of simple system. (b) Root locus of 
system with lead network. 


Fic. IV-29. 


frequently stated to be a “‘phase margin’’ of about 50 
e., 0 = 150” to 130" at 
In the example of Fig. 


to 50° at “gain cross over” 
the frequency where AM = 1. 
IV-33, @ = 150° at w = 1.35 where M = 0.45. In 
order to have AK. = 1 at this frequency, the gain 
would be K = 2.22. 

In comparing (a), the Bode or Nyquist presentation, 
with (b), the Evans root-locus method, we note that the 
Bode or Nyquist diagrams permit the direct use of ex- 
perimental data; thus they do not require the ap- 
proximation of the frequency-response curves of Fig. 
IV-15 by those of a second-order system. 

On the other hand, the Evans root-locus method pre- 
sents directly a complete picture of the stability and 
transient response characteristics that are most 1m- 
portant for aircraft control; for airplanes are rarely sub- 
ject to purely sinusoidal inputs. Its merits are as 


follows: 


AERONAUTICAL 
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(1) The root locus gives the roots of the closed-loop 
system directly and, by a simple calculation, the tran- 
sient response. 

(2) The degree of stability can be read from the root 
locus directly. 

(3) Complicated systems can be set up in such a 
fashion that the effect on transient response and sta- 
bility of changing any parameter in the problem can 
easily be visualized. 

(4) There is no ambiguity or confusion concerning 
interpretation of the plots even for complex systems 
having any number of roots and poles in either the right 
or left half of the complex plane. 

For comparison, Fig. IV-34 presents the Bode dia 
gram for the roll-control system of Figs. IV-2 and 
IV-28. 


(F) Illustrations of Simple Feedback Systems 


(1) Roll Stabilization.—The roll control of an air- 
craft or missile may be approximated by the differential 


equation 
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Fic. IV-30. (a) Autopilot control of a statically unstable air 
frame. (b) Autopilot control of a statically stable air frame 
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AERODYNAMIC 


I,.(d*¢/dt?) + Le.(de/dt) = ké (IV-73) 


where 
g(t) = roll angle 
(p) = Laplace transform of ¢(¢) 
6(t) = aileron angle 
D(p) = Laplace transform of 6(¢) 
@ = moment of inertia about the longitudinal 
axis 
Es = coefficient of damping in roll (usually de- 


noted by L,; we shall not use the common 
aerodynamic definition of p = dg/dt to 
avoid confusion with the Laplace operator 
p) 


In Laplace transforms, Eq. (IV-75) reduces to 


(,.p? + Le.p).®(p) = k.D(p) — (IV-74) 


assuming a system initially at rest. Thus, the airplane 


transfer function will be defined by 


P( p) k k Le 


"1.G,(p) = = = 
eid D(p) I,.p?>+Le.p p.l+ pT) 
(IV-75) 


In many practical cases, the aerodynamic damping Lg 
is so small that it may be neglected; then 


K,G\(p) = k £*) 


The simplest servo-control system would be one that 
measures the angle « = ¢; — g by means of a position 
gyro and gives aileron deflection 6 proportional to e. 
Then 6 = u(¢g; — go), and the transfer function for the 
servo would be A»G2(p) = uw. The transfer function for 
the two systems in series would be 


KG(p) = (uk/Le).1/[pUl + p7)] 


while that for the feedback system would be 


K Gp) = KG(p)/[1 + AG(p)] 


The function G(p) has a simple pole at p = 0 and p = 
=f l T | } 


cates that, as the gain A is increased, the new roots move 


The root-locus diagram (Fig. IV-35a) indi 


toward each other and then along the vertical line 
p = —(1/27)) to +7 © and — 1 W. 


therefore always stable as long as Lg is positive. 


This system is 
The 
damping is, however, extremely low; the damping ratio 
is 

2a T 7™ i 


f=] 2V uk! 


A better approximation to a servo would allow for a 
small time delay 7), between the input signal ¢« and the 
resultant aileron deflection 6. This relation could be 


represented in terms of a first-order system by the 


transfer function A,.G,(p) = uw/(l1 + T)p). This 
would give the open-loop transfer function 
KG(p = i(k Lo |p (1 + pT)) } [u (1 + Tp) | 


and for the feedback system 
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The root locus in this case starts from the poles al p 0, 
p = —(1/7)), and p -1/T),. 


this cubic system in Fig. IV-24, and we note that at a 


We have already seen 
gain A > K, the system becomes unstable. In order to 
extend the regime of stable operation over a wider fre 
quency regime, a lead and lag network is introduced, 
represented by the transfer function 


me 1+ elop 6, + 13d es + 14p 
K.G.(p) = a = = 
1 + T»p 1+ 736 1 + 7yp 
Lag ea Lead 


The experimentally determined transfer function for the 
hydraulic system 1s approximated by 
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Fic. IV-33. Phase-amplitude diagram (Bode plot) for . § : a iets : 
: functions are shown in Fig. IV-36. The corresponding 
G(iw) = io + iw)I( + tw)/10)" root-locus plot is shown for the conditions of high dy- 
namic pressure (Fig. IV-37) and low dynamic pressure 
(Fig. IV-38), respectively. The servo gain selected for 
K;.Gi() Me the high dynamic pressure condition is A, = 0.40; for 
3. 3(P = . : +): - ~ wh 
die the low dynamic pressure condition, A; = 5.3. The 


(p?/wo”) + 26(p/wo) + | 


and thus, the resultant transfer function AG(p) shown 
in Eq. (IV-63) is obtained. The root locus for this sys- 
tem is shown in Fig. IV-28. It is seen that this system 
is stable above a gain of A.,i¢ (;) and below a gain Ait (2). 
The design range varies from a gain of A = 9.75 at low 


dynamic pressures up to a gain of A = 670 at high dy- 
namic pressures. The roots corresponding to the latter 
case are shown as ~, po, .. . pi. 

(2) Pitch Stabilization.—The previous example 


showed how an airplane could be roll stabilized with a 
position gyro and a servo feedback network. Let us now 
study how a rate gyro could be used for pitch control. 
Let us assume that the speed of the airplane is kept con- 


corresponding roots for the frequency equation are 
shown in Table IV-2. 

(3) Computation of System Response. 
a simple feedback system whose open-loop transfer 


Let us consider 


function is 


TABLE I\V-2 

High Dynamic Low Dynamic 

Pressure, Pressure, 

Condition K = 0.40 K = 53 
pi —28 + 50: —16 + 437 

p» —28 — 50i —16 — 43: 

p —27 + 22.51 —27 + 221 

ps —27 — 22.51 —27 — 221 

ps — 8 — 2 
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Fic. 1V-36. Pitch stabilization system with rate gyro 
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Fic. IV-37. Root locus for longitudinal dynamic stability 
analysis-high dynamic pressure. (Note: zero of origin due to 
rate gyro cancels pole due to airplane. ) 


i Sa) ae SP ss 
Le — a) AP — Ge)... AP — Ge 
(IV-76) 


G(p) = 


and whose closed-loop transfer function is determined 
from the relation 


K.G.(p) E a 
| + KG(p) 
G,(p) = N;(p) ss (p — S1).(p — Se) (Pp — Sm) K 
D.(p) (p — fr).(p — po) (p — pn) Ks 


(IV-77) 
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Fic. IV-38. Root locus for longitudinal dynamic stability 
analysis-low dynamic pressure. (Note: zero of origin due to 
rate gyro cancels pole due to airplane. ) 


It will be noted that the zeros of G,(p) are the same as 
the zeros of G(p). The response due to a step-input 


yi(t) = h(t) is then 
" N,(p,) 
(Tt) = | G,(0) pPri 
mi) = KG. + 9 prep) © | 
= K, [Ag + > A,e’"| 
; l 
(IV-7S) 


These terms may be interpreted graphically from the 
positions of the poles and zeros in the p-plane. Con 
sider, for example, a third-order system having the 
closed-loop transfer function 


Gs(p) = 1/[(p — pi)(p — p2)(p — ps)] 


where p) = —1 +1, po = —1 — 1, ps = —10 and thus 
(pi — po) 2 condi pi — ps = 9 O5e° agit pe — Pps 
9.05e"°*. We find 
mi | 
A, = G,(0) = - = 
pi. po. Ps 20 
| N(p1) | 
7 pi. Ds" (pr) pi.(pPi — p2)(pi — ps) 7 
I ~tel2o 
é 
25.6 
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(a) 


P4-P,* V2e7iiss° 


P, 
.— > 


P, Bel s— , 


Py- Pp * J2e 138° 





(b) 


iw-p,*100.5¢ ‘°*? 
iw-p, = 990ei18* 


iw- Pp, = 101 ei 89.4" 


c 
ic. 1V-39. Vectors Pe terms of the response 
_ N,(p2) a l oe 
py. Dy'(p2) pa. (D2 — pr(p2 — ps) 
a. i129 
25.6 
l 


_ N (ps) i” l — 
PsD,' (pa) Ps-(Ps — pi)(bs — pe) $20 


The vector lengths p,, p; — po, pi — ps, ete., are easily 
measured in the p-plane (Fig. IV-39a). It is obvious by 
inspection that A; will be much smaller than A; and A» 
because the vector lengths to the point ; are much 
larger than the others. Thus, the term A;e”" will be of 
small importance in the response compared to the first 
two terms. 

If there is a zero in the p-plane, the effect is to reduce 


the effectiveness of a pole close to it. Thus, if to the 


above three poles we add a fourth py; = —2.00 and a 
zero S; = —2.01, then we would find 
AY —1 
ie ag te = 1.005A,y 
bs Pipops 
4 , <a Pr I 
Ay = . = 
Ps — Pr Pr (Pi — pe)(pi — ps) F 
1.007e°"- A, 
BP ee ee “i 
= ps — po pe: (p2 — pr)(p2 — ps) ; 
1.007e-°"""A, 
§ —_— l 
in. = 0.999A; 


Ps — Pps ps.( bs — pr) (ps — pe) 
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tt Si — pa ; l _ l 

Ps (pbs — pi)(ps — pro) (ps — ps) 3200 
The vector lengths are shown in Fig. IV-39b. This 
example illustrates that the coefficient A, is extremely 
small, and, if it were possible to have the zero s; and the 
pole ps coincide, they would completely cancel each 
other, and the coefficient A, would be zero and the 
motion e?” would be completely eliminated. 


twl 


If the exciting force is sinusoidal (y; = e’“’), then the 


response is 


N;(pr) “os 


wit) = K,. | Gut ») et 4+ 
: ” a (p, — Iw) D,'(p,)° 


r= 1 


steady-state vibration transient vibration 


Xn 

= K,(A,.e"% + > A,e*”) 

r=1 

It is clear from the preceding example that the co- 
efficients A,, A;, A, ... A, can again be determined 
graphically from the measured vector lengths and direc- 
tions in the p-plane. If, for example, w = 100 rad. per 
sec. in the previous example, then 


A, = G,(w) = = = ; 
(tw — pi).(tw — p2).(tw — ps) 
A, = —_ a 7 
(py — iw)D,'(p1) (pi — tw). (pi — pe). (pi — ps) 

- N.(p2) 7 

' (po — iw)D,'(p2) (pe — tw)(pe — pi) (pe — py 
ee 7 
(ps — iw)D,'(ps) (bs — iw)(ps — pr)(Ps — pr) 
where 


iw — p» = 99e'°* 
iw — p2 = 101e"** 
iw — ps = 100.5e"*” 


The vector lengths are shown in Fig. IV-39c. If there 
were a zero of the system near the exciting frequency 
w = 100 rad. per sec., then the response would be 
small; if the zero coincided exactly with the exciting 
frequency, then the response would be completely elimi- 
This case is of great practical engineering im- 
It corresponds to the so-called “‘tuned vibra- 


nated. 
portance. 
tion damper’’ discussed, for example, by Den Hartog 
in reference 6. It corresponds physically to adjusting a 
system so that the exciting force acts at a node of the 
vibratory system. 

The system response to any other input can be com- 
puted analogously. Thus, if, for example, it were de- 
sired to compute the output y(t) for an exponential 
input y,(t) = e * then, since Y,(p) = p/(p + 8B), we 
would have 


Vo(p) = [p/(p + B)].A.G,(p) (IV-79 


and 
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an” 
yo(t) = Bs.e~” + > B,.e?" 


r=1 


(IV-SO) 


The coeflicients Bg, and B, could be written down by the 
partial fraction rule [Eq. (A-10)] and evaluated graphi- 
cally in terms of the vector lengths in the p-plane as 
before. 

The error of the response for large values of time ¢ 
may be determined for a single feedback network for 
various simple inputs by applying the rule [Eq. (A-16)] 
that lim e«(f) = lim E(p). This rule applies if the 


— pod 


(a) For a sample step input: 
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closed-loop transfer function G,(p) has no poles on the 
imaginary axis or in the right half of the p-plane. If 
the transfer function is of the type 

G(p) = 1/[p"(p — pi)(p — pe) ] 
we find from Eq. (IV-56) 


E(p) = Vi(p)/{1 + 1/[p".(p — pid(p — pe)]i 


Consider the following examples. 





yi(f) = h(t) Vi(p) = 1 
l si 
; : I : fn = 0 
lim ¢(f) = lim x = {1+ (K/pips) (IV-81) 
oe ae “ 0 ifn > | 
P \P — pi)(P — pe) 
(b) For a simple velocity input: 
yi(t) = Vo.t Vi(p) = Vo p 
oo ifn <1 
Vo Vo “¢ r 
lim ¢(f) = lim re. = ifn = 1 (IV-82) 
t—>o p—>o 4+ K l + (K pipe) 
pb" (p — pi)(b — po) 0 ifn>| 
(c) For a simple acceleration input: 
yi(t) = (1/2)at? Vi(p) = (1/2)(a/p?) 
oo iin <2 
: ; (1/2)(a/p?) (1/2)a : ote 
lim ¢(¢) = lim : = = ifn = 2 (IV-83) 
Ross po, K 1 + (K/pip») | 
p"-(p — pi)(P — pr) 0 ifn >2 


(G) Multiple-Loop Networks 


The autopilot system for an airplane consists of 
multiple-loop networks, since the six equations of mo- 
tion of an airplane involve interactions between the 
pitch, up and down, and fore and aft motion, and also 
between the roll, yaw, and side motions. If the airplane 
is in any maneuver except straight flight, the first three 
longitudinal motions interact with the latter three 
lateral motions. In addition, the sensing instruments 
may measure mixed combinations of pitch, yaw, and 
roll angles when the airplane is undergoing maneuvers. 

In order to design such multiple-loop networks, it is 
usually desirable to consider a single loop at a time 


using simplified assumptions. After the system param- 


eters are tentatively fixed, it is then instructive to 
expand the analysis one loop at a time using the root- 
locus method for each loop. The principles involved 
are illustrated in Appendix (D) for the case of longi- 
tudinal stabilization of an aircraft. 

Another method that is recommended by Hanny'® 
is as follows: Write down the Laplace transforms of 
all of the equation and arrange them in determinental 
form. Then insert a cut in the network, and introduce 
an input signal ¢,(¢) at one side of the cut. The condi- 
tion for closure is that the output signal ¢,(¢) at the 
other side of the cut be equal to the input signal. 
Fig. IV-40 illustrates schematically such a control sys- 


tem for a steam turbine. When an analog computer is 
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Fic. 1V-40. Multiple-loop network representing a steam-turbine 


regulator.!* 


available, it is ordinarily most convenient to utilize it 
to study such multiple-loop networks experimentally. 


(V) FLIGHT-TESTING TECHNIQUES 


Flight tests are an essential step in the development 
of automatic control systems, both from the standpoint 
of verifying the performance characteristics of the 
system as well as the achievement of high reliability 
under the actual conditions of pressure, temperature, 
In the present section, we propose to 
namely, the deter- 


vibration, etc. 
consider another type of flight test 
mination of the aerodynamic coefficients in order to 
compare them with the coefficients estimated on the 
basis of wind-tunnel tests with stationary models. 
The basic problem of this flight technique is to 

determine the aircraft transfer function 

KG(p) = [Oo(p)]/10:(P) | 
or 

KG(iw) = [Oo(tw) |/[0;(iw) ] 
where 0,(p) represents the Laplace transform of a 
control-surface deflection and 0)(p) represents the re- 
sulting motion (linear or angular). This technique 


Milliken of the Cornell 
Laboratory in an extensive series of tests on a North 


was used by Aeronautical 
American Aviation B-25J airplane to determine the 
transfer functions for the longitudinal motion of the 
airplane.*® The elevator was moved +1'/2° sinusoi- 
dally at various frequencies from | to about 5.5 rad. per 


sec., and the resulting normal acceleration » and 
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attitude angle 6 were measured. These results are 
shown in Figs. V-la and V-1b, including a comparison 
with the theoretically estimated values. This technique 
of introducing a sinusoidal input deflection 6 and 
measuring the output motion 6 is the most straight. 
forward one in principle. Practically, this method has 
been rather tedious, in that it required about three 
hours of flight testing to cover the frequency range for 
one flight condition and required many hours of evalua- 
tion to reduce the data. Even so, the phase angles 
between the motions could not be deduced accurately 
enough from the separate measurements. There appear 
to be good potentialities, however, for the development 
of electrical systems that measure directly the magni- 
[Oo(1w) |/[0;(tw) ] 
With such equipment, it 


fication ratio 
between output and input. 


and the phase angle 


would then be simple to determine the transfer func- 
tion for various flight conditions from subsonic through 
transonic and up to supersonic speeds. 
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(b) 
Fic. V-1. Airplane response for a sinusoidal elevator deflection 
56 = dye! = elevator deflection. @ = Het(o! + 4) = pitch angle 
n = noet(wt + «) = normal acceleration. XX X X = experimental; 


= theoretical. (According to W. F. Milliken, Jr., Cornell 
Aero Laboratory, B-25J airplane tests. ) 
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In order to proceed from the transfer function to the 
determination of the principal aerodynamic coefficients, 
it is necessary to make assumptions about the magni- 
Thus, the flight- 
test technique provides a check only on the principal 
aerodynamic coefficients. This is sufficient for verifying 
Further 


tude of the secondary coefficients. 


the assumptions of the autopilot analysis. 
development of the flight-test techniques is desirable 
to separate more clearly the effects of aeroelastic de- 
formations and the effects of aerodynamic time lags 
on the motion of the aircraft. 

Other methods for determining the aircraft transfer 
function for longitudinal motion have been to apply an 
elevator displacement (a) as a step function and (b) 
as an impulse (Fig. A-4). Milliken experimented with 
the step-function input. 
test data is much simpler; however, the interpretation 
of the data is more difficult since the control charac- 
teristics (particularly the hinge moments) depend not 
only upon the elevator angle but also upon the rate 
and acceleration at which the angle is applied. Another 
difficulty with a step-function input, particularly in 
roll and vaw, is that the mean flight path may be altered 
radically, and thus the validity of the linear theory 
may be exceeded. Furthermore, the step input excites 
the low-frequency modes of the system much more than 
the high-frequency modes, and thus important high- 
frequency modes may be overlooked or located inac- 


He found that securing the 


curately. 

These last two difficulties are appreciably reduced by 
the use of a pulse-type input. In this case, a surface 
deflection is introduced and quickly removed, and the 
airplane remains near its equilibrium position. In 
addition, it may be seen that the various modes of 
oscillation of the aircraft are excited more uniformly. 
In the limit, as a perfect pulse (infinitely large ampli- 
tude and infinitesimal time duration) is approached, all 
modes are excited equally. Practically, the speed with 
which the surfaces may be moved and the resulting 
fidelity of the input pulse determine the range of fre- 
quencies over which the pulse spectrum is flat. 

Several pulse shapes have been considered, the first 
resulting from a rapid surface deflection, followed by a 
period in which the deflection is held constant, and 
finally, by a rapid restoration of the surface (Fig. 
V-2a). 


0,(p) = a(l — e °’) 


Oo( Pp) 


Input: 
Output: 


The transfer function is given by 


KG(p) = Oo(p) - B0(b) 
O(p) a(l—e ®’) 
In the low-frequency range, the process of reducing the 
data to obtain AG(p) may be simplified by considering 
the output to be made up of a number of pulses similar 
in form to the input, as illustrated in Fig. V-2b. With 
this approximation 
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A(t) A ho, O<t< ¢ 
A(t) al hy, rx F< 2r 
A(t) A he, 2r <t < Sr 
and 
0,(p) = cm bl —e "2 
n= i 
Thus, 
| m 
G(p) =- > me 
a n 1 
? 
ot t 
a) STEP INPUT: b) OUTPUT: _ 
“ = O(t)=h, o<te 
a[h (t) h(t T)] eh pcteey 
=he 2T<t<3T 
etc. 
DEFLECTION 








c) GRAPHICAL COMPUTATION d) ACTUAL DEFLECTION 
OF TRANSFER FUNCTION 
Fic. V-2. Determination of transfer functions from rectangular 
pulse inputs. 
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Determination of transfer function from triangular 
pulse inputs 


Fic. V-3. 
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and 


m 


te h ew net 
Rev 


1 
an = ] 


KG(iw) = 


Since the factor e~'"*” represents only a phase shift, 
the transfer function is obtained by adding the h, with 
the proper relative phase as shown in Fig. V-2c. The 
required vector addition may be accomplished simply 
with either punch-card or elementary analog equip- 
ment. 

Practically, it is not possible for the surfaces to be 
deflected instantaneously as indicated; a finite deflec- 
tion time is required, as shown in Fig. V-2d. Since the 
frequency spectrum of the input is constant over the 
largest frequency range when the time interval of the 
pulse is a minimum, it appears desirable, from the stand- 
point of uniform mode excitation, to employ a triangular 
pulse in which the slope of the sides corresponds to the 
maximum rate of surface deflection. This technique 
was first applied to the experimental determination of 
the aircraft transfer function by Seamans, Blasingame, 
and Clementson of the M.I.T. Instrumentation Labora- 
tory (Dr. C. S. Draper).” In this case the input signal 
of Fig. V-3a has the transform 


a | , a | 
-—o(] — Qe"? 4+. @ 2) = -.-. 
tT p tT p 


(l—e ")? 


0,;(p) = 


Again, for low frequencies, the transfer function may be 
obtained most simply by considering the output to 
consist of a number of pulses similar in form to the 
input. Thus 

ho l hy l 


7 
Oo(p) = eisai ee 4 
tT p 7 p 


e *(1—e *)? + 


In this case, the output is approximated by a series of 
straight-line segments connecting the points at time 
equal to 0, 7, 27, ete., as shown in Fig. V-3b. As in the 
case already considered, 


m 


KG(iw) = pe esa 

eos 

A transfer function obtained in this manner is satis- 
factory for low frequencies. For the higher frequencies, 
it is necessary to consider the fine structure of both 
the input and the output and to use the exact expression 
for G(iw). This is particularly important when it is 
desired to determine the influence of factors, such as 
aerodynamic time lags, which primarily affect the high- 
frequency response of an aircraft. Before this effort is 
considered, however, it appears desirable to study the 
problem of noise in the measuring equipment and its 
possible effect on the results. The amplitude of the 
high-frequency response is normally small, and delicate 
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In this 
regard, it appears that some improvement might be 


measuring equipment is required to detect it. 


obtained by averaging the results of a number of tests, 
thereby reducing the effect of random disturbances, 

In general, the pulse method is superior to the 
steady-state techniques: Simpler equipment and con- 
siderably less flight time are required to obtain the 
results. From the standpoint of accuracy, however, 
neither method appears to offer substantial advantage. 
Considerable improvement is required before either 
method will be suitable for extending our knowledge of 
aircraft dynamics. 


(VI) FUTURE PROBLEMS AND POTENTIALITIES 


(A) Future Problems 


(1) Nonlinear Systems.—The theory and practice of 
automatic control systems are still in the early stage of 
development. Thus far, only linear systems have been 
thoroughly analyzed. Nonlinear systems have been 
studied only for especially simple cases. 
of the theoretical work in this field are available.*’: 
The advent of the new computing machines may make 
it possible to carry out systematic studies of nonlinear 


Summaries 


83, 85, 86 


elements and thus provide background information 
that will permit the intuitive synthesis of nonlinear 
controllers of desired characteristics. An energetic 
attack of this problem is under way by the Office of 
Air Research, U.S.A.F., under the direction of Dr. 
Frank W. Bubb. 

(2) Treatment of ‘‘ Noise.” 
in electromechanical networks is similar to that of 
turbulence in a fluid flow. Taylor,** in his early studies 
on the statistical theory of turbulence, introduced the 
Wiener, Kolmogoroff, 


The problem of ‘‘noise”’ 


idea of correlation coefficients. 
and others have carried over the use of these concepts 
to the treatment of noise in electrical filter networks. 
With the aid of this theory, the communication engi- 
neers have learned to design optimum linear filter net- 
works that reduce to a minimum the effect of the noise 
Many of the 
“smoothing networks” 


84, 87, 89, 90 


relative to the signal.'® 
principles for the design of 
which have been developed by the communication 
engineers are being taken over in the design of automatic 
control equipment. The application of these tech- 
niques to controller design promises to be one of the fer 
tile fields for improvement in the accuracy of response 
of control systems. 


(3) Optimum Controller Design.—The design of a 


corrective electrical network for a control circuit which 
will have prescribed characteristics is one of the means 
of improving the response of control systems. The 
general techniques for the design of such networks 
have been worked out by the communication engi- 
neers.*:'! These corrective networks are variously 
called ‘‘controllers” or ‘‘equalizers.’’ These more com- 
plicated networks may have application to aircraft 
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control systems, although probably only for specialized 
applications. The reason is that the transfer function 
of the airplane is not a fixed quantity but varies over 
wide limits, and thus it is not generally worth while 
to go to extreme pains to design an optimum controller 
circuit for any one design point. 

There are, however, a number of simple tricks that 
may well be incorporated in the synthesis of the con- 
troller. Moore describes two methods for designing a 
controller that converts the equations of the unequal- 
ized system into the desired form. The first alters 
the system equations by the addition of certain terms; 
the second method inverts certain terms. An example 
of the former occurs in the case of aircraft control 
where undesirable rolling moments occur as the result 
of rudder deflection. Here, the controller anticipates 
the rolling moment associated with rudder deflection 
and automatically calls for a compensating aileron de- 
Another example occurs in the case of an 
resulting from rolling 


flection. 
undesirable yawing moment 
angular velocity. Here, the rudder is deflected by an 
amount necessary to compensate for the calculated 
yawing moment based on rate gyro measurements of 
roll rate. The use of the inverting controller is illus- 
trated when there is a necessity for compensating auto- 
pilot gain for variations in dynamic pressure. Here, 
the system gain is increased as the dynamic pressure 
decreases so that the overall system characteristics are 
as independent of dynamic pressure variations as is 
economical. 

The methods of direct synthesis involve measurement 
of functions of dependent variables of the system to 
compensate for effects of these functions on the outputs. 
However, since such controllers are inside the system 
feedback loops, they must be used with caution, be- 
cause any uncertainty in the form of the system equa- 
tions may cause overcompensation and actual insta- 
bility. 

Another method for improving the performance of 
control systems is to use an open-cycle controller to 
schedule the approximate value of an output and to 
use a closed-cycle system to supply the vernier control. 
Such a system, proposed by Moore, is shown in Fig. 
VI-1. 
given in reference 1S. 
tion due to rolling angular velocity can be compensated 


A more extensive discussion of this same topic is 
As an example, the yawing mo- 


in open-cycle control by anticipating the rolling rate 
that will result from application of a given aileron de- 
flection and calling for a proper amount of rudder de- 
flection to eliminate the yawing moment resulting from 
such a roll rate. 


(B) Future Potentialities for Aircraft 


There are many applications of automatic control 
systems to aircraft beyond the automatic flight and 
landing systems discussed so far. Among the more 


important ones are the following. 
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OPEN CYCLE 
CONTROLLER 
(INVERSE OF 
UNALTERABLE 
ELEMENTS 
AND LOAD) 
LOAD TORQUE 
(IF OPEN CYCLE CONTROLLER IS 
PERFECT THERE WILL BE NO STEADY 
STATE ERROR) | 
+ 
8; RROR| SERIES LOOP UNALTERABLE 8. 
= * CONTROLLER [-*> ELEMENTS Sueur 
CONTROL LOOP 
PARALLEL s 
CONTROLLER 
Fic. VI-1. Combination open-cycle closed-cycle system (J. R 


Moore ). 


(1) Automatic Cruise Control of Propulsion System. 

Most propulsion systems have two or more parameters 
that can be varied to produce a given thrust. During 
cruise, the important problem is to operate the system 
at the optimum specific fuel consumption expressed in 
pounds per mile. An open-cycle, closed-cycle system 
appears well adapted to schedule the operation at 
approximately the correct speed and altitude and to 
depend thereafter upon the direct measurement of fuel 
flow and distance in order to adjust the parameters for 
optimum performance. The adjustment of the controls 
must be carried out consistent with the limitations of 
the engine and the aircraft. The various adjustable 
parameters and the limiting factors are shown in 
Table VI-1. All of these power plants suffer under 
certain conditions from instabilities. The first four may 
experience a pulsation associated with the compressor; 


TABLE VI-1 


Principal 


Adjustable 
Operating Limits 


Parameters 


Propulsion 
System 
Piston engine r.p.m.—by propeller (1) Maximum r.p.in 
and propeller pitch (2) Maximum mani 
Manifold pressure—by fold pressure 
fuel flow 


Turbojet r.p.m.—by fuel flow (1) Maximum r.p.m 
Turbine inlet tempera- (2) Maximum turbine 
ture—by exhaust inlet temperature 
nozzle area (3) Pulsation limit 
(4) Combustion blow 


out limit 
(1) Maximum r.p.m. 
(2) Maximum turbine 
inlet temperature 
(3) Pulsation limit 
(4) Combustion blow 
out limit 


Propeller tur-  r.p.m.—by _ propeller 
bine pitch 
Turbine inlet tempera- 
ture—by fuel flow 


(1) Maximum tem 
perature 

) Pulsation limit 

) Blowout limit 


Maximum tempera 
ture—by fuel flow 
Pressure recovery —by 
exhaust nozzle area 


Ram-jet 


9 
2 
o 


( 
( 


Rocket Fuel flow (1) Maximum pressure 
Fuel ratios (2) Combustion in 
stability 
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the last four may experience instabilities associated 
with combustion. There is a good possibility that auto- 
matic control systems may either prevent operation in 
these unstable regimes or, under certain conditions, 
perhaps permit the damping of these oscillations by 
introducing a feedback link into the system. 

(2) Automatic Stabilization of Aircraft Under Stall 
Conditions.—An aircraft is ordinarily unstable in roll 
as it approaches the stall. It tends to fall off on one 
wing or the other. It should be possible to eliminate 
this instability by an automatic control system that 
applies a rudder deflection to counteract the unstable 
roll. 

(3) Automatic Stabilization of Helicopters. 
pointed out by von Karman” in 1925 that helicopters 


It was 


are dynamically unstable. This prediction was verified 
when helicopters were successfully flown in the 1930's 
in that they had to be continuously controlled. There 
are aerodynamic methods for stabilizing the helicopter. 
However, the automatic stabilization by an autopilot 
may well prove to be the generally most flexible system. 

(4) Gust Alleviation. Many 
carried out on the possibility of reducing the effect of 
gusts on aircraft both from the standpoint of permitting 


studies have been 


reduced design load factors and to improve passenger 
comfort. Most of the previous studies were carried 
out on a semiempirical basis. With the advent of 
rapidly responsive servosystems, it may well be possible 
now to achieve a better measure of success on this 
problem. 

(5) Automatic Stabilization of Seaplanes and Hydro- 
planes Riding on the Water.—Seaplanes experience a 
longitudinal dynamic instability called porpoising. It 
should be possible to devise an automatic control sys- 
tem that will damp out these oscillations. For example, 
a hydrofoil surface arranged either ahead or behind the 
main float could be actuated by an automatic control 
system. 

Hydrofoils have also been repeatedly proposed as a 
more efficient method of cruising either ships or air- 
craft in the water at high speeds (by Tietjens and 
One of the problems of hydrofoil systems is 
Auto- 
matic control systems should be able to provide sta- 


others). 
again the achievement of longitudinal stability. 


bility either with the use of a control surface in front or 
behind the main lifting hydrofoils. 


CONCLUSIONS 


We have seen that the current high-performance 
aircraft are flying in the transonic speed regime that is 
the most difficult from the standpoint of stability and 
control. In order to overcome the requirements of 
high stick forces, power controls had to be introduced. 
In order to increase the precision of maneuvers, auto- 
matic control systems were required to supplement the 
human pilot. The fighter airplane with automatic and 
power controls therefore represents an intermediate 
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version between the conventional aircraft of Worid War 
II and the guided missiles of the near future. 

The application of automatic control systems to air- 
craft promises to bring about the most important new 
advances in aeronautics in the future. The past 25 
years have seen a tremendous development in the fields 
of propulsion, aerodynamics, and structures. The next 
25 years will probably see comparable major develop- 
ments in the application of automatic control systems 
to stabilize and control accurately and reliably the 
aircraft of the future from its initial take-off to its 
landing. The student of aeronautics would, therefore, 
be well advised to include in his program a study of the 
fundamentals of dynamics and of servomechanisms, 
The practicing aeronautical engineer will find it neces- 
sary to gain an understanding of these principles in 
order to design the optimum aircraft systems. 


Appendix (A) 
METHOD OF LAPLACE TRANSFORMS 


A) Solution of Linear Differential Equations 


The theory of Laplace transforms represents one of 
the most useful tools of dynamic analysis. A complete 
discussion of this theory is given in a number of ex 
cellent texts.**~“’ In order to follow the analysis of the 
present paper, the following condensed summary may 
be of help. The more important rules are illustrated by 
some problems whose solution ts well known to the aero- 
nautical engineer. 

The Laplace transform* of a function y(t) is denoted 
by ¥(p) and is defined as follows: 


Y(p) = » f e ” y(t).dt (A-1) 
0 


It is assumed that y(t) = O for t < 0. 
ordinarily denotes the time and is a real number. 


The variable / 


The variable p may, however, be a complex number. 
The above integral defines Y(p) as long as the integral 
converges. For example, if y(t) = e~™, then 


V(p) = p ee" dt = ag re Se 
0 = 2 


*In the mathematical literature it is conventional to define 
the Laplace transform 


F(s) = | wt) .e . 
0 


7 


1 a-wtrt 
yt) = = / F(s).e% . ds 
ot Ja a 


This form has advantages of greater symmetry in its trans- 


and its inverse 


formations compared to our }(p)/p. Our form |[Eq. (A-1)] has 
the advantage, however, that the dimensions of Y(p) are the 
same as those of y(t) since p has the dimensions of (1/time). 
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At the upper limit, the function e p—avt + Qast— « 


if (p — @) is positive—i.e., if p > +a. At the lower 
‘e~9' _, | ast— 0, and thus the integral defines 


limit, é 
V(p) = p/(p — a) for p > +a. As part of the defini- 
tion of the Laplace transform of y(t) = e”, we define 


V(p) = p/(p — a) forp<a 


This part of the definition follows by analytic continua- 
tion (i.e., by using the same analytic expression of the 
function) rather than as a consequence of the integral 
definition, Eq. (A-1). 

In a similar manner, it is possible to substitute other 
simple functions y(/) under the integral sign and to 
compute the corresponding values of Y(p). A short 
table of such transforms is given in Appendix (B) (in 
both the engineer’s p and the mathematician’s 5s 
notation). It is understood that y(¢) has the values 
shown only for ¢ > 0 and that y(t) = O for ¢ < 0. 

If the function y(¢) or its derivatives have a discon- 
tinuity at the origin, it is necesary to define more care- 
fully what we mean by y(0)—1.e., the value of y at a 
time ¢ = O, and y’(0)—i.e., the value of dy/dt at time 
t= 0,ete. If, for example, an elevator rising at a uni- 
form velocity Vo picks up a load at position fp», the 
position of this load at time ¢ is shown in Fig. A-lc. 
We shall, as a matter of definition, agree that 


y(0) = y(t) att =O 
y’'(0) = dy/dt att = 0 —e 


where € is a very small quantity that approaches zero; 
ie., we take the values of these functions to the /eft 
Thus, in the present case y(0) = hy 
Let us now study the 


of the origin.* 
and y’(O0) = 0 (not y’(O) = V»). 
Laplace transform of a derivative dy/dt by carrying 


out an integration by parts: 


,, dy : 
» fi e — -dt = p.[e™.y(t)]° — 
0 ora C ne an 


"eal u v 


dv 
p / y(t).(—p).e ” dt 
Aa =~ —_——_ 


v du 


The value of e”. y(t) > O ast — ©& if y(t) either is 
finite at f = © or diverges less strongly than e”. The 
value of e ”. y(t) > y(0) as t > 0, and, as agreed above, 


we shall take y(0) = y(O — e€). 


p / e ny = 
J0 dt 


and similarly, 


Thus, we get 


—p.y(0) + p.V(p) (A-2) 


* It is also possible to build up a logically correct system with 
other definitions of y(0) and y(0). For example, references 5 
and 32 use y(0) = y(O + €) and y'(0) = y(O + e) and define the 


Laplace transform of dy/dt as s.F(s) — y(O + e). 
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ia mo UY , ; ' 
p e "-—~—-dt = —p.y'(0) — p?.y(O) + p?. Y(p) 
0 dt? 5 
(A-3) 


With these two rules, we can convert a linear differential 
equation with constant coefficients into an algebraic 
equation. For example, the differential equation of a 
mass-spring system with velocity damping and an input 


force [;(t) is 
m(d*y/dt?) + c(dy/dt) + ky = fi(t) (A-4) 

with the boundary conditions 
y = y(0) 


‘ ati =O0- 
dy/dt = y'(0)§ ’ , 


We multiply this equation by pe” and integrate from 


t=O0Otot= oo 


! - a t+ / oe dy it 4 
mp Mss. “¢ cp e -—-@. +- 
0 dt? f /J0 dt 
kp [ e ” y(t).dt = p / e ” fi(t).dt 
/0 J/0 


We define as before 


Y(p) = rf e 
0 
Fi(p) = p J e ” fi(t).dt 
J 0 


and evaluate the integrals according to our definitions 


pt y(t) dt (A-5) 


and 


(A-6) 


m(|—p.y’(0) — p*.y(0) + p?.¥(p)] + 
cl[—p.y(0) + pY(p)] + k.V(p) = Fi(p) 


or rearranging terms 


(mp? + ch + k).Y(p) = Fi(p) + mpy’(0) + 


(mp? + cp).v(0) (A-7) 
and solving for 
. Fi(p) 
V(p) = - 
mp?>+cp+k 
m.p.y'(O) + (mp? + cp).y(0) 
2 = (A-S8) 


mp? +cp+k 


In order to find y(t), the inverse transform of }(/), 
we first resolve ¥(p) into terms of the type A,.|p 
(p — p,)], for which we know the inverse transform is 
A,.e’". We use the following “partial fraction rule’ 

If G(p) = N(p)/D(p) where V(p) and D(p) are poly 


nomials in p, we have 


N(p) = ado + a1.p + Q2.p? + + aap" 
D(p) = (p — pi) (p — fe») (p — pn) 
— p, in D(p) 1s 


and m < n and none of the factors p 


repeated; then G(p) may be separated into fractions as 


follows: 
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ae TTT HHI Fil?) | KG(p) [Yo(P) 
No AHHH 
} Ui, BLOCK DIAGRAM 
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(c) 





oT t 

h (t-7)= UNIT STEP AT t=T } a. UNIT IMPULSE AT ts 7 
(e) (f) 

Fic, A-1. 





Definition of singularities. 


P P 


G(p) =Ay+ Ai: + A»: + ...+ 
= p — pr 
A . p 
p — Pn 
where* 
Ay = N(0)/D(0) = G(0) 
and 
A; = N(p1)/[pi.D'(pr)] 
A>» = N( pe) ‘[pe.D'( po) ] 
A, = N(pu)/lPn-D' (pn) | 
or alternately 
as N(pi) 
et py. 1. (Pi — pa)(Pr — Pa)... (Pr — Po) 
* To find A,, note that 
N(p) _G(p)_ At 
p.Dip) bp pb 
Ao A» Aa A. ® )- ( _ D 
( . +4 )- (p) — (bp — pide(P) 
b Pp — pe Pp — Pn (p — pi).(p) 
o(p)/P(p) 
and that 
N(p1) = A;.P(p1) 
1.e. 
A, = N( Pi) ( p;) 
where 
P(pi1) = pi(pi — f2)...( pi — bn) 


which could also be written 


d 
(pi) = pr. D = p,.D"(p,) 
(Pi) = pr E wo ae pi. 


The other coefficients are obtained analogously. 
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— N(p2) 
p2.(p2 — pi).1. (pe — ps)...(po2 — pn) 
ae N (Pn) 


Pn-(Pn ie Pi) (Pn a Po)... (Pn ca Pn 1) l 


where D’(p) = dD/dp. 
form of g(t), we find 


g(t) = Ay +- A,.e?" + As” + : 


If G(p) is the Laplace trans. 


. + A,eP” 


This expansion theorem may therefore be stated as 
follows: 
If G(p) = N(p)/D(p), this may be expanded into 


. ~  N(p,) p 
G(p) = GO) + >& ; (A-9 
,=1),-D'(p,) p — P, 


and its inverse transform is 


n N(p,) 
g(t) = G(O) + —— @F  (A-%5 
2 p,.D (p,) 


where i, po, ... Pp; ... Pp» are the roots of the frequency 
equation D(p) = 0. 
Let us consider the application of this rule to Eq. 


(A-8). We note that Y(p) consists of two parts: 


‘ Fi(p) 
Yo(p) _ ms P - 
mp? + cp+k 
Laplace transform of output mo- 
tion due to the input force /;(¢) 
; m.p.y'(O) + (mp? + cp).y¥(0) 
Yo(p) = b p i 7 


mp? + cp +k 
Laplace transform of output mo- 
tion due to the boundary condi- 
tions y(0) and y’(0) 


In the analysis of automatic control systems we shall 
be primarily concerned with Yo(p). We denote the 


ratio 
output/input = ¥o(p)/Fi(p) = KG(p) (A-11 
where 
K = Y,(0)/F;(0) (A-12 
is a constant called the “‘gain’”’ 
G(p) = 4 ' =A?) = transfer function (A-13 
K_ F;(p) 
E the above example AK = 1/k and G(p) = 


l 
(m/k)p? + (c/k)p +1 
by a block diagram as shown in Fig. A-1d. 
The solution of Eq. (A-10) for two important special 
cases of input force may be written down immediately: 


4 
| We denote this relationship 
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(a) If f(t) = F,,.h(t) where F,, is a constant and 
h(t) is the unit step function, then F;(p) = F,,. From 
Eq. (A-11) we then have 

Yo(p) = F,.K.G(p) 
and assuming G(p) 1s given in the form of Eq. (A-9) 
; “ N(p;) et 

y(t) = F,,.K | Geo + 2. f ais (A-14) 
; , 1 P, D (p,) 


is the motion due to the input force f;(f). 
(b) f(t) = F,.e”, then 
Fi(p) = Fry. p/(p — tw) 
and from Eq. (A-11), 
Yo(p) = Fu. |p/(p — tw)|.KG(p) 


This expression is expanded similarly to Eq. (A-9) 
considering Gi(p) = [p/(p — tw)|G(p). The new 
numerator is Ni(p) = p.N(p). The new denominator 
is D,\(p) = (p — tw)D(p), and thus we find 


yo(t) = Fn .K. Glia) et 4 





steady-state vibration 


n N( : 
~ ae eortl (A-15) 
ro 1 (pp, — tw).D'(p,) J 


transient vibration 


Similarly, the motion y,(¢) due to the boundary condi- 
tions may be found by expanding Y;,(p) in partial frac- 
tions. The motion y,(f) is independent of the exciting 
force /;(t) and must be added to yo(t) if the resulting 
motion is desired. In many automatic control problems 
the initial boundary conditions may be neglected; in 














these cases y,(t) = 0. 
4 
hit) 
F) t 


h(t)* STEP INPUT y,(t) = OUTPUT 











t(t) y(t) 
t t 
f(t)" GENERAL INPUT GENERAL OUTPUT 
e Y(t)=¥ (o-€)+ 
e-T d¢(T) “ 
| y,( — dt 
Fic. A-2. Superposition of elementary solutions. 
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The similarity between the coefficients of Eqs. (A-14) 
and (A-15) is apparent. If the zeros and poles of G(p) 
are known, it is possible to construct graphically the 
function G(iw) by the methods outlined in Section (IV). 

The second term of Eq. (A-15) represents the tran- 
sient oscillation as a result of the initial disturbance, 
and the first term represents the remaining vibration 
after the transient has damped out. 


(B) Limiting Value of y(t) as t — 


Another special case that follows from the above ex- 
pansion theorem is the limiting value of the function 


y(t) ast—+ ©. We note that 


V(p) = p/(p — a) 


for y(t) = e”. 
When a = —0 + iw where bisa positive real number, 
then 
. at ° t t+ tw) 
lim e” = lim e ”.e*'“ +0 
t— « t— « 
and 


lim p/(p — a) ~0 


pO 


When the function Y(p) can be expanded by Eq. 


(A-9) in terms of 


: n p 
Y(p) = Ay + A,: 
f Py p oa P, 


we therefore find that 


lim Y(p) = lim (4. + > A,- P ) = Ay 
o> = £p, 


p> p70 r=1 


This is the same as 




















git) yp(t) 
~+t ~t 
8(t)=UNIT IMPULSE y,(t) OUTPUT 
INPUT 
i 4 Llp 4 ~~ + 


GENERAL OUTPUT 


t 
yitdey(one+fyelt-T)-ftt) dT 
Superposition of elementary solution. 


f (t)= GENERAL INPUT 


Fic. A-3. 
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the time r will be represented by h(t — r) (Fig. A-lc), 
TYPE at y(t) Y(p) Similarly, a unit impulse located at time 7 will be 
t " ; 
p v6 represented by 6(¢ — 7) (Fig. A-If). 
STEP INPUT | || | h (t) : 
AT ORIGIN =| gM, 
t 
(D) Superposition Theorems 
STEP INPUT a nee e-eT 
iT =7, ‘ ‘ ‘ * - 
ten oO 7 Methods for assembling solutions from these elemen- 
= —— tary functions are given by the following rules. 
| RECTANGULAR | h (t-7,)-h (tT) e-P Ti _6-P Te (a) A dynamic system has the response yi(f) as the 
PULSE: T, TO Ta} o = result of a unit step input A(t) at? = 0. If an input f(t) 
__- —_——_———. is applied starting at time ¢ = 0, then the resulting 








Fic. A-4. Construction of Laplace transforms of a rectangular 
pulse. 















































oe ee | 
TYPE SHAPE y(t) i! Y(p) 
oe! — a 
| -— 
} , y(t) 
LINEAR INPUT » 
STARTING AT A-+ a 
t=0 | P 
—_ rs | | 
| 
| y(t) FOR’ t+ <T,: 
LINEAR INPUT gitre A -py | 
STARTING AT ; | =er" | 
t+7, | 1/ P 
, My FOR t> 7: 
7, y(that-7,) | 
=a cece 1 
(t) | FOR t<7, 
4 | y(t) = 
FOR 7,<t¢7, Se eee 
| y(tealt?7) [pe “pe : 
7 T, FOR t > T; 
‘ y(t) «tT, - 7) 
ee — | 
- Tr 
y(t) 
SIMILAR TO 
ABOVE SPT A o-PTs| 
/ EXAMPLE P P 
Ty T: % | 
}— ——— - T = es | 
y(t) ha = A(T-T) = A(Ts-T) | | 
} | 
| AlePTi-2e-Pl2+eP Ts 
i | 
= ee 














Fic. A-5. Construction of Laplace transforms of a triangular 
pulse. 


Ay (A-16) 


nN 
lim y(t) = lim (4. +> A,e*) = 
i—0 t—0 r=1 

if and only if none of the roots p, is on the imaginary 
axis or in the right half plane. 


(C) Singularities 


The method of Laplace transforms presents in a 
simple manner the response of a system when the input 
is a unit step function /(¢) (Fig. A-la) or a unit impulse 
6(t) (Fig. A-lb). The latter is the equivalent of a rec- 
tangular impulse of amplitude A and duration from 
t = O0Otot = 0 + esuch that the area under the curve is 
A.e = 1. The width e of the pulse is made shorter, 
the amplitude larger, such that the product remains 
constant at 1. The limiting case as « ~ 0 represents 
6(t). It will be noted that a step function starting at 


motion at any time ¢ is 


> 


: df(r) 
y(t) = yO — €) + w(t — 7)°- dr 
0 dr 


—_—_—— 


Ay {A-17 


This theorem is obvious by inspection of Fig. A-2. 
The input force f(r) is broken up into a number of hori- 
zontal rectangles. [df(r) /dr]. Ar. 
The length of time since the initiation of the force 
The resulting motion at time ¢ is 


Their height is Af = 


increment is ¢ — rf. 
the sum of all of the elementary increments. 

(b) A dynamic system has the response yo(f) as 
the result of a unit impulse 6(/) at ¢ = 0. If an input 
f(t) is applied starting at time ¢ = 0, then the resulting 
motion at any time ¢ is 


bd 


y(t) = y(0 — e) + | yo(t — 7).f(r).dr (A-18) 
0 


This theorem is obvious by inspection of Fig. A-3. 
The input force is broken up into a number of vertical 
rectangles of height f(r) and width Ar. The effect of 
each rectangle at time / is (7) Ar times as much as that 
of a unit impulse placed at the same point 7. A unit 
impulse at the point 7 corresponds to an ‘‘influence 
function”’ The overall 
effect is the sum of these contributions. 


of unit force at the point r. 


(E) Construction of Simple Pulse Shapes 


It follows from the definition of the Laplace trans- 
form that if 


rf y(t).e "dt = 
ee 
rf y(t — r).e ”.dt=e 


Y(p) 
then (A-19) 


V(p) 


This relation is useful for finding the transforms of some 


of the pulse shapes used in the analysis of flight-test 
data. Figs. A-4 and A-5 illustrate the method of build- 
ing up response curves in terms of simple singulari- 


ties. 
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Appendix (B) 


TABLE OF LAPLACE TRANSFORMS 


Equation 
Reference 


v(t) \(p) F(s) Number 
ent p/(p —»X 1/(s d) (B-1) 
pb ) S lw 
elw : - P (p + Ww) = (B-2 
> — tw p* + w § w- 
cos wt p?/(p? + w?) s/(s w? (B-3) 
sin wl wp/(p* + w?) w/(s? + (B-4) 
> (p +A) + iw ( Ay +1 - 
¢ V+ tu b ~ = p b = =f cee 4 = (B-5) 
(p +A) — tw io Tap te (s + A}? + w* 
b(p +X) is + A . 
Mi-cos wl PAF . = : - (B-6) 
(p +r)? + (s +A)? +a 
pw w - 
e—-sin wl ; <s - ; ; (B-7) 
(Sf + AP Pw (s + A)* + w* 
(1 e—)/r 1/(p + X) 1/[s(s + A)] (B-8) 
Unit step: h(t) l l/s (B-9 
t 1/p 1 /s? (B-10) 
t"/n! 1/p’ 1/p"* (B-11) 
te-M P = (B-12) 
(p +X) (s +A)? 
Unit impulse: 4(¢) p l (B-13) 


Appendix (C) 
THE EQUATIONS OF MOTION OF THE AIRPLANE 
The equations of motion of the airplane are 
F = mi(dV/dt) + & X V] 
and 
CL = (dA/dt)+axHi 
where F is the resultant externally applied force due to 


gravity, the thrust, and the airforces; L is the corre- 


sponding externally applied torque; V is the resultant 
airplane velocity; @ is its angular velocity; and /7/ is 


its angular momentum. 


(a) Body Axes 


(1) Longitudinal equations: 


E £U dé ors we OX pz nq 
m dt 0a di 20 COS % = a a 
Ee =) uid | eS OZ» * 
= dt dt alas teal Oa 7 
d*0 OMn 
Y = Qa 


dt? Oa 


These equations that were first expressed in moving 
body axes by Euler (1707-1783) were linearized by 
Routh” with the assumption that the body carries out 
small perturbations about a mean equilibrium condi- 
Bryan and Williams (1903)*' first applied them 
Bryan 


tion. 
to study the dynamic stability of an airplane. 
expanded this theory in his book (1911) into what is 
still the classic treatise on the subject.’? Milliken®® 
recently presented an excellent summary of the history 
of dynamic stability of aircraft, and further review of 
this topic therefore appears unnecessary at this time. 

The resulting equations of the linearized theory are 
given below——(a) for body axes and (b) for so-called 
stability axes. The definitions of the terms are given in 
Figs. C-1 and C-2. 


OX B 4 OX B 7 ; + OX B 5 
ia” lCUr ~~ 
Zp OZ» OZ» 
T i 
oo a a 
2M, dd Mp WM, . . Mp 
, r bp 
O(d6/dt) dt t Ou oiled oT ud 05x ? 
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(2) Lateral equations: 
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1951 


(dB. dy ) | OV, OY, 
l o( = 0 — _ ‘0S Oy = B bp 
m| a se: — 0B ~ 
d* = OLp OLn dy OL pn do 4 OL . ' OL; 
7d? dB O(dy/dt) dt O(d@ dt) dt Obr ™ 06, me 
I d’y = ONe ONp dy ONez do ON» n ON: : 
* dt? rey} O(dy/dt) dt O(dgd dt) dt bx - 06, ™ 
(b) Stability Axes 
(1) Longitudinal equations: 
| eae ‘ ox, Pi OX, OX,,, OX,. 
m 1(9 — a) cos (A — a) | = : , 
a g a) cos a da a oo Uu oT + 5, Or 
g (= 7 + (6 in (8 | OZ, 4 OZ, 4 OZ, rT OZ, : 
m 0 = 4 = ) sin oo 0) = Be 4 
. “a: =" * ‘ x" ee T° uu 
d°6 OM, OM, 1 On, OM, . OM, dé 
/ = u — nee ; 
dt? da. | om rv’ em (de dt) dt 
(2) Lateral equations: 
fap *) | or, oY, 
Ud -f ey ‘OS (A — 0) = 3 
m| "as ?™ ws ag | *™ 
1? £* OL OL 1 OL, 1 Le s 
eS ee ae a n dy he ca 5, 
dt? dt? oye} O(dy/dt) dt O(dg@ dt) dt ObR 064 
ay d*p ON, ON, 1 ON, 1 ON, ON, 
r. i. = — A = + = + Or a O+ 
dt? dt? aye} O(dy/dt) dt O(d@ dt) dt Obp 06, 


Since modern aircraft and missiles may take off and 
climb almost vertically, the previously used assump- 
tions that lift = weight and drag = thrust are no longer 
generally valid. It is then most convenient to study the 


dynamic stability at a discrete number of points 
along the trajectory and to substitute numbers directly 
into the above equations rather than to transform them 
in terms of dimensionless coefficients. This method is 
illustrated in the example, Figs. I1V-37 and IV-3s. 

One of the questions at the present time is still the 
magnitude of the ratio (nonstationary air force quasi- 
steady air force). Studies of this question have been 


started by Miles*® and others.°! 


For airplanes equipped 
with autopilots, an exact solution of this problem is not 
critical, since most of the effective damping is intro- 
duced by the autopilot. 

A simple similarity analysis of dynamic stability 


yields some interesting results. Let us consider the 


simplified pitching equation, writing @ = y + a where 
y = flight-path angle. Then we have 
ee ae ee ee 
y° 9 ws q° - a — . a ( ) 
dt? dt . 


where 


VW - Cu, (1 Z)eV?.0 1 V 
M. = Cua: (1/2)pV?- i? 1 
and 
Pm = Mass density of aircraft 
p = mass density of air 


/ = a typical length dimension of the airplane 


Comparing this equation with the spring-mass-damping 
system, 


d°0 dé 
| 2¢ Oy Wy” oe = fif 

dt? d 

we notice that 
: I" 

wo = N (1 2)¢ Ma \ Ph Pm" 
and 

¢ = Cu, \ (1 2) Crzq Vv P/ Pm 


The dimensionless ratios that we find are the following, 
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TABLE C-1 

Altitude, Ft 0 10,000 20,000 30,000 40,000 50,000 60,000 

t/ts V/ p/p 1 0.86 0.73 0.61 0.50 0.39 0.31 

(a) Aerodynamic coefficients, Cy, and Cy,, which 
depend largely upon other dimensionless ratios—-angle 
of attack, Mach Number, Reynolds Number, etc. 

(b) Frequency ratio, F = w)//V. 

(c) Damping ratio, ¢ ~ V p/p, when the aero- 


dynamic coefficients Cy, and Cy, are constant. 

The conclusions that follow are that for dynamic 
similarity of maneuvers: 
(a) Fo = v= 
~ 1/land uw = V. 
Thus, a missile 6 in. in diameter flying at 2,000 
m.p.h. has 100 
times as high as a geometrically similar airplane of 60 
Small high-speed 


wl constant. This implies that 


Wo 


a natural frequency approximately 


in. in diameter flying at 200 m.p.h. 
missiles represent, therefore, the most difficult autopilot 
problems. In order to ease the autopilot problems, it is 
desirable to reduce the static stability coefficient Cy, 


to as near zero as possible. (See Fig. IV-30, a and b.) 


(b) ¢ = V —/pm. 

If we compare two identical airplanes, the first at 
sea level (s/) and the second one at altitude flying at 
the same dynamic pressure, then the natural frequency 
w is the for both for Mach Number 


effects). The damping coefficient ¢ is, however, much 


same (except 


less at altitude since €/f.) ~ V p/ps:. This reduction in 
the damping coefficient with altitude is shown in Table 
C-1. 


Appendix (D) 


THE “PInN-BALL METHOD” OF SOLVING HIGH-ORDER 


ALGEBRAIC EQUATIONS 


Evans’ root-locus method is essentially a graphical 
method for solving the frequency equation D(p) = 0. 
This is an algebraic equation, and thus the root-locus 
method is generally useful for solving high-order alge- 
braic equations. Evans calls the following method the 
“pin-ball technique’ because of its analogy to the well- 
known game: A new ball (representing a pole) entering 
upon the scene disturbs all of the other balls (poles) 
more or less. 

Let us illustrate the procedure for simple algebraic 
equations and follow the previous convention, indicating 
a root of D(p) = 0 by a small circle in the p-plane. A 
root of D(p) = 0 is a pole of the general response 


Y(p) = N(p) D(p). The first-order equation 
D\(p) = (pb + 3) = 0 
has as its root p = —3. This is shown in Fig. D-la. 


The second-order equation 


D2(p) = p* 4. 3p + ao = 0 
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may be written 


Dip) = (bp + 3)p + aw = 0 
or 
[p(p + 3)]/a2 = —1 mes 
The left-hand side had roots at p = 0 and p = —3. 


The root locus is shown in Fig. D-1b for any value of 
a). Suppose a2 = 2; then the magnitude condition is 
satisfied at p = —2 and —1 and 

Dp) = p? + 3p + 2 
The third-order equation 


D3(p) = (p? + 3p + 2)p + az 


may be written 
[((p + 2)(p + 1)pl/as = —-1 


The root locus for this case is shown in Fig. D-lec. 
The values of the parameter a; can be easily computed 
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Coordinate systems associated with aircraft stability. 
body axis and subscript S = stability axis 


_K FORCE 


Fic. C-1. 
Subscript B = 








L,M,N = MOMENTS 
m= MASS OF AIRFRAME 
I= MOMENT OF INERTIA 
THRUST 
PRODUCT OF INERTIA 
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CONDITIONS 


Axis system and notation 
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for various points along this locus. At a; = 2 the roots 
are p) = —0.25 + 0.827, po S —0.25 — 0.821, ps & 
—2.5. 

This process can be continued. The fourth-order 
equation 


Dip) = p* + 3p? + 2p? + 26 + mw = O 
may be written 
(bp — pr(p — po)(p — ps)pl/a = —-!1 


The root locus for this case is shown in Fig. D-Id. For 
a, = 0.9 the roots are 


pi’ = 0.81 

po’ = —0.87 
p;’ = —2.4 
pa’ = —0.6 


Every time the order of the equation is increased by 
one, a new pole enters the picture. The effect is to give 
a displacement to all of the existent poles, the amount 
of the displacement depending upon the value of the 
parameter a. 

It is convenient to use a similar technique for the 
study of the motions of coupled systems. Each time a 
new degree of freedom is introduced, the characteristic 
determinant has more terms that are combined step- 
by-step into a single factored term. This procedure is 
illustrated by the following steps. 

(a) The airplane is mounted in a wind tunnel with 
an axis of rotation permitting a rotation in pitch; the 
resulting equation of motion ts: 
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Dy = (p? + 0.25p + 0.8)6 = 0 


The corresponding oscillatory ‘‘short period’’ roots 
Pi and py» are shown in Fig. D-2a. 

(b) The airplane is permitted an additional degree 
of freedom—namely, an up-and-down motion. The 


resulting equations of motion are: 


(p? + 0.25p + 0.8)0 — 0.87 = 0 
0:2 6— (p+ 0.2)7y = 0 


The expanded form of the characteristic determinant 
is then: 


D, = —(p* + 0.25p + 0.8)(p + 0.2) + (0.2)(0.8) = 0 


The general procedure is to find the zeros of the sum 
of two terms by plotting the locus of values of p which 
make the ratio of the two terms equal to —1. Since 
the signs of the terms are opposite in this case, the ratio 
must equal |, as shown in Fig. D-2b. 

One of the roots is at the origin, as can be checked by 
noting that the constant terms cancel. Denoting the 
new roots as py and p», the determinant can now be 
expressed as: 


D, = —p(p — pu)(p — par) 


(c) An additional degree of freedom is introduced in 
that the airplane is permitted to move fore and aft. 
The new equations of motion are: 


(p? + 0.25p + 0.8)0 — 0.8 y + 0. u=0 
O.2 06— (p+ 0.2)7 + 0.01 u = 0 
0.4 6+ 0.6 y + (p + 0.03)u = 0 


This determinant should be expanded by minors 
about the third column in order that the previous 
determinant appears as one of the minors. 


Dz, = (p + 0.03)D, — 0.01[(p? + 0.25p + 0.8) X 
0.6 + 0.4 X* O.8] 


As 


The second minor A: is simply a quadratic having 
zeros Ss; and s, and the locus plot for the ratio of the two 
Note that the 
phugoid oscillation is introduced by the horizontal 


remaining terms is shown in Fig. D-2c. 


degree of freedom but that the short period roots are 
shifted slightly to pi, and pr». 

(d) If a rate gyro is introduced to control the 
elevator and a perfect servo is assumed, the equations 
of motion are: 


(p? + 0.25p + 0.8)6 — OSy+tO0. ut+ C;.6=0 
0.2 6— (p+ 0.2)7y + 0.0lu + 0.6= 0 
0.4 6+ 0.6 y + 
(p+ 0.03)u+ 0.6 = 0 
Ko.p 6+ 0. y+0. u- $= 9 


The expansion of the determinant by minors about 
the fourth row or column gives: 
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Ds = (-1)D2 — GsKop[—(p + 0.2)(p + 0.03) — 
(0.6) (0.01) } 


As 


Denoting the zeros of the new minor A; by s; and 
s,, the root-locus plot versus A,C; as the variable is 
shown in Fig. D-2d. Selecting A,C; to be 1, the roots 
are at Pis, Pos, Ps, and ps. 

(e) Damping of the phugoid terms can be achieved 
by introducing elevator control proportional to upward 
velocity. The resulting equations of motion are then 
as follows: 


(p? + 0.25p + 0.8)0 — OS y +0. ut+ C;.6 = 0 
0.2 6 — (p+0.2)y +0.0lu + 0.6 = 0 

0.4 6+ 0.6 y+ 
(p+0.03)u + 0.6 = 0 
Ky.p gé+ K,y+0. u- 6 = 0 


The complete determinant is then: 
D, = Ds + K,C;[0.2(p + 0.03) — (0.4)(0.01)] 


The locus versus A, is shown in Fig. D-2e. 

The above example, which was worked out for the 
author by W. R. Evans, shows nicely the effect of vari- 
ous terms upon the roots of the frequency equation and 
indirectly therefore upon the transient response. 
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Discussions of the Lecture 


Preston R. Bassett, President, Sperry Gyroscope 
Company, Division of The Sperry Corporation: 
I wish to compliment Dr. Bollay on his most excellent 
and comprehensive paper. When I say comprehen- 
sive, I mean something much greater than what we 
used to consider the scope of this field. The addition 
of the two great new zones in aerodynamics—namely, 
transonic and supersonic-——has tremendously increased 
the complications of stability and control problems. 
Dr. Bollay’s paper is interesting in showing the rela- 
tionships and the differences of the problems and 
methods of solution in the three speed fields. 

As I have personally been associated with the prog- 
ress in this field from the time in 1914 when Law- 
rence Sperry won his first prize for the gyro stabilizer 
in Paris, I would like to call attention to the changing 
approaches to this problem which have been made 
during the 40 years of its history. 

When Mr. Sperry was developing the stabilizer 
from the years 1910 to 1915, his objective was to sup- 


ply an accessory to the airplane which would make it a 
more practical device. I mean by this that the 
stability of the early designs of airplanes had so little 
margin that it was only the extreme and continuous 
alertness of the human pilot which kept them in the 
There 


air. It seemed that some aid was necessary. 


were two approaches—Sperry’s approach of a gyro 
accessory and the aerodynamic approach of better 
stability designed into the airplane itself. During the 
war years from 1915 to 1920, so much was learned 
about stability in airplane design that the gyro ac- 
cessory was found to be nonessential, and inherent 
stability reached a point where the human pilot 
could at least relax a little. 

But by the late 1920's, airplane performance had 
improved to the point where the duration of flight 
and the range of the airplanes were so great that the 
element of fatigue entered the problem and automatic 
control again came into the picture—not so much to 
give improved stability to the craft as to relieve the 
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pilot from long continuous hours of playing servo- 
motor. I think that Wiley Post’s flight around the 
world alone in 1933 with an automatic pilot first 
demonstrated the advantages of automatic flight 
from the fatigue point of view. 

Both military and air transport planes have used 
automatic flight to considerable extent from that 
time on. But up to World War II, automatic flight 
consisted of two parts: (1) any design of airplane 
built for human control and (2) an added accessory 
called an automatic pilot. During World War II, 
with the increased performance and speed of the 
military planes plus the increased demand for ac- 
curacy of the flight control, a new concept became 
necessary, and that was the treatment of the aero- 
dynamic plane, its control system, and the automatic 
flight controller as a single system. Great progress 
has been made along this line during the past ten 
years. I highly endorse the statements that Dr. 
Bollay has made throughout his paper of the impor- 
tance of the control designers’ and the airplane design- 
ers’ working more closely together and treating their 
problem as a single overall system. 

When it comes to problems in the transonic and 
supersonic and those of guided missiles, we have 
reached the point where automatic controls have be- 
come such an integral part of the unit that it would be 
quite foolish to design one portion without relation 
to the others, and, in guided missiles, the human 
element has already dropped out of the picture. 

It is interesting to note, therefore, that we have 
more or less completed a cycle in the field of automatic 
control and are back to the original objective of in- 
creasing the stability and improving the control of 
the entire aerodynamic unit. This time, however, 
the function has become vital—tt 1s here to stay. 

Dr. Bollay's paper is impressive in pointing out the 
numerous analytical methods of approach to the 
complicated problems. I wish to comment on a few 
specific points that Dr. Bollay mentions in his paper. 
He predicts that the future interceptor will carry out 
automatically all the necessary functions in inter- 
ception, with the pilot’s job being that of switching 
from one phase to the next. We agree with this and 
would call attention to the fact that in our automatic 
instrument landing approach this system of having 
the human pilot switch from phase to phase is al- 
ready in successful use. The so-called ‘‘push button” 
flight of the Air Force from the United States to Eng- 
land in 1947 was also a demonstration of these pos- 
sibilities. 

Dr. Bollay comments on the importance of aero- 
elasticity in automatic control. Mention is made of 
FSO. In our experience in controlling the FSO auto- 
matically, we have not run into aeroelasticity troubles, 
but, in larger high-speed planes, the effects may be 
serious. We might add that although aeroelasticity 
usually makes the control system problem more dif- 
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ficult, in some cases it does provide an advantage: 
for example, the twist of wings which varies with 
speed and altitude and is in a direction to reduce con- 
trol effectiveness, may be used to minimize the need 
for gain changes under varying flight conditions of 
speed and altitude. 

On this point it might be mentioned here that the 
big problem in automatic control at present is the 
design of a system that will work under all flight condi- 
tions with widely changeable parameters. This is 
analogous to the aerodynamic problem of designing a 
plane that flies at high speed and is capable of land- 
ing at low speed. 

Dr. Bollay comments upon the use since 1940 of the 
techniques devised by electrical engineers for the 
analysis of automatic control systems, such as the 
work of Nyquist and Bode. This is important, and it 
can be pointed out that, although these methods 
might be somewhat new to the majority of aero- 
nautical engineers, they were first extensively used 
in the design of fire-control equipment, then later 
applied to automatic airplane control. 

Dr. Bollay gives a description of the use of analog 
type simulators in the solution of automatic control 
Examples are given based upon the use 
This has been an out- 


problems. 
of the Reeves type computer. 
standing improvement in the technique of pretesting 
control systems. It is interesting to note, however, 
that prior to electronic analog computers, similar 
studies were made by a mechanical analog we called a 
synthetic airplane. 

It might be mentioned that a simulator solution is 
only as good as the data used and is only as informa- 
tive as the questions asked of it. In the past, the 
use of simulators has been limited by the knowledge 
of aerodynamic parameters. Even with the data fully 
known, it is necessary for the full use of a simulator 
for the engineer to have such imagination and experi- 
ence as to be able to anticipate all the conditions that 
might occur in flight and to be able properly to define 
them in mathematical form. 

Dr. Bollay mentions a method of open-cycle control 
which improves automatic control. As an example, 
he states that the yawing motion due to rolling angu- 
lar velocity can be compensated open cycle by antici- 
pating the rolling rate that will result from applica- 
tion of a given aileron deflection and calling for a 
proper amount of rudder deflection to eliminate the 
yawing moment resulting from such a roll rate. This 
has been used for some time in our own automatic 
pilots by a simultaneous introduction of aileron and 
rudder signals to initiate a coordinated turn. 

Dr. Bollay mentions among the future problems 
the automatic control of helicopters. Great progress 
has been made on this problem during the past year, 
and I can state that the solution of the problem has 
proven much easier than anticipated. 
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AERODYNAMIC STABILITY 


Philip Colman, Chief Aerodynamics Engineer, Lock- 
heed Aircraft Corporation: Dr. Bollay is to be com- 
plimented on an excellent summary of where we have 
been and where we are going in the research and ap- 
plication of the science of stability and control. 
There is little one can criticize or add to such a com- 
plete treatment of the subject. I do have a reaction, 
however, which results from the title of Dr. Bollay’s 
paper, which could be slightly altered tosay: ‘‘Weare 
gradually abandoning aerodynamic stability and be- 
coming dependent on automatic control.” 

In the past we have “pulled little tricks’ to im- 
prove what we could produce in the basic aerody- 
namic stability and control aspects of our airplanes. 
These ‘‘tricks’’ included such items as: 

(1) Aerodynamic balance and geared balance tabs 
to tailor aerodynamic hinge moment to a desired 
variation. 

(2) Flying tabs for pilot actuation which could 
float the major control surfaces to the desired position. 

(3) Spring tabs to incorporate the benefits of bal- 
ance and flying tabs by a mechanical spring relation- 
ship for airplanes with increased flight-speed range. 

(4) Variable mechanical advantage control sys- 
tems to utilize the favorable hinge moments in a dif- 
ferential aileron control system or to take advantage 
of the difference in required elevator or rudder surface 
angular positions required at low and high speeds. 

(5) Surface ‘“‘bungees’’ to create a constant me- 
chanical spring hinge moment to be balanced by a 
variable aerodynamic hinge moment to attain ap- 
parent level flight longitudinal stability or proper 
lateral ‘‘feel.”’ 

(6) Control system ‘‘bob weights’’ to simulate the 
favorable mechanical hinge moment of the ‘“‘bungee’’ 
with the added apparent stability in accelerated ma- 
neuvers. 

Our schemes have expanded with the incorporation 
of boosted control systems. These presently incor- 
porate proportional control ‘‘feel’’ in the majority of 
systems, but the trend toward irreversible systems 
with ‘‘synthetic feel’’ is clearly discernible. Thus, it is 
only one more step to the fully automatic control sys- 
tem in which the pilot only monitors the mechanical 
brain. 

The dynamic stability problems became manage- 
able several years ago when we perfected the aircraft 
designed to fly at half the speed of sound. These 
aircraft encouraged the establishment of dynamic 
stability criteria, which we are presently finding in- 
creasingly difficult to meet on the newer “improved” 
Our troubles stem from: (1) increased 
increased flight-altitude 
higher wing loading; (4) 


airplanes. 
flight-speed range; (2) 
range; (3) smaller wings 
revised wing shapes, changing the relationships of the 
fundamental aerodynamic parameters; (5) reduced 
effective aspect ratios and redistribution of airplane 
weight components increasing the importance of iner- 
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tia factors; (6) variations in aerodynamic parameters 
in the transonic range; and (7) markedly increased 
flexibility of airplane structural components. 

Thus, aerodynamic stability is deteriorating and 
automatic control is becoming increasingly neces- 
sary. We must now learn to pull “‘big tricks’’ with 
automatic control to attain the ability to fly at higher 
speeds and under more adverse conditions than we 
have presently been forced to operate under. 


Harry J. Goett, National Advisory Committee for 
Aeronautics: The Institute is to be commended both 
for its timely choice of subject for the Fourteenth 
Wright Brothers Lecture and for the selection of a 
man uniquely qualified to discuss the subject. Dr. 
Bollay’s coherent synthesis of the concepts and 
methods from at least three fields 
servomechanics, and electrical engineering 
most useful to the aeronautical engineer in his attempt 


aerodynamics, 
will be 


to understand the problems of the expanding field of 
automatic control of aircraft. 

In his introduction, Dr. Bollay remarked that the 
airplane is going through a period of transition be- 
cause of its increase in speed and the resort to powered 
and automatic controls. It also became evident, as 
the talk progressed, that the aeronautical engineer 
engaged in the field of stability and control is going 
through a transition. He is being called upon not 
only to master more complex mathematical tools, but 
also to take for his own the techniques and methods 
of the servo designer and the electrical engineer. In 
this process, there is the danger of losing a grasp of 
the manner in which the various parameters influence 
the physical behavior of the final airplane-autopilot 
combination. This situation promises to get worse 
rather than better as the added complexities of non- 
linearities or added degrees of control are introduced. 
It is interesting to note that one of the factors of 
importance in Dr. Bollay’s evaluation of the various 
methods of analysis and synthesis is the degree of in- 
sight they give as to the importance of the parameters 
involved. Incidentally, in this regard, an aerody- 
namicist finds some consolation in Dr. Bollay’s com- 
ment that there appears little reason to introduce a 
concept of decibel per octave into the aeronautical 
dynamics literature. 

Dr. Bollay’s remark that in the use of Routh’s 
criterion one loses insight into the physical behavior 
as the various parameters are changed deserves com- 
ment. Somewhat of an extension of this method is the 
construction of stability boundary charts that, in the 
past, proved useful in the analysis of unregulated air- 
craft and which are proving equally useful at the 
present time in the analysis of automatically controlled 
aircraft. The advantage of using these charts is that, 
in addition to having the stability boundaries defined 
in two of the important system parameters, such as 
control gearing and static margin, it is possible with 
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little extra computation to superimpose grid lines of 
constant damping and constant period. The pres- 
entation of all this information in the form of one com- 
pact plot is advantageous in many ways, and to pro- 
duce this information using the root-locus or any of 
the frequency-response methods would require much 
more time and effort. The necessity of ‘‘translating”’ 
the meaning of points on the p-plane to points on 
the ¢-plane is 
other methods have their own peculiar advantages, 
as pointed out by Dr. Bollay, and the relative merit 


also eliminated. Of course, the 


of the various methods depends on whether they are 
being used to analyze or synthesize a system. 


Possibly an even more basic question than the fore- 
going is that relating to just what physical behavior is 
desired of the airplane-servo combination beyond the 
obvious requirement of stability. The phrase ‘‘opti- 
mum airplane-autopilot combination” is frequently 
used, but I personally have never encountered any 
quantitative this The 
analysis techniques permit assessment of various per- 


definition of “optimum.” 
formance criteria of which the following four appear 
to be the most significant: (a) speed of response, (b) 
steady-state error, (c) time to damp, and (d) amount 
of initial overshoot. Any numbers that have been 
assigned in regard to these characteristics have come 
from functional requirements in fields other than that 
of aircraft control and, therefore, are deserving of re 
newed scrutiny. The situation we find ourselves in 
reminds me somewhat of that which existed around 
We 
knew then that we wanted a stable slope of Cy vs. 
C,, but, beyond that, the aeronautical engineer had 


1935 in regard to static stability requirements. 


to select his own requirements as to the magnitude of 
such quantities as stick force per g or stick force varia- 
tion with speed. This vagueness was finally resolved 
by the development of “‘flying qualities requirements” 
pioneered by the work of Soulé and Gilruth of the 
N.A.C.A. and eventually adopted by the Services. 


The parallel to our present situation in the airplane 
servo field is obvious. We need a more explicit ap- 
preciation of the relation of the several oscillation 
characteristics previously mentioned to the ultimate 
function of the servo-controlled airplane. In asking 
for this, I realize that the human pilot, whom we pre- 
viously considered the somewhat variable and tem- 
peramental quantity to be satisfied, is relatively unt- 
form when compared to the wide variety of functions 
that the 
perform. 


servo-controlled aircraft will be asked to 
The question is, therefore, how long do we 
have to remain in this stage of somewhat intuitively 
adjusting the rates of response, overshoot, or time to 
damp of our airplane-servo combination, after which 
we wait for the final judgment on the basis of such 
things as miss distance or compatibility with associ- 
ated radar or automatic computer equipment as 
determined by actual flight tests. 
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W. F. Milliken, Jr., Manager, Flight Research De- 
partment, Cornell Aeronautical Inc.: 
The subject of this year’s Wright Brothers Lecture js 
The speaker has 
chosen to focus primarily upon the analytical aspects 


Laboratory, 
obviously one of considerable scope. 


of dynamic analysis or what is commonly referred to 
as servomechanism theory. It is important that this 
material become a part of the tools of the stability 
and control designer, for many of his problems are no 
longer completely solvable by aerodynamic means 
alone. The paper serves to point this up in a timely 
manner, and Dr. Bollay is highly deserving of con- 
gratulations. 


Ic is indeed fortunate, as the speaker points out, 
that these 
able to us from the electrical engineering field and 


mathematical methods have been avail- 
that we have been saved from a large part of their de- 
There are, however, unique areas for 
research the itself 
must pursue if it is to be able to wisely and economi- 


velopment. 


which aeronautical profession 


cally avail itself of what the servomechanist has to 


offer. These areas evolve around two items: 1) 
the rational determination of the extent to which 
servomechanisms are desirable and necessary for a 


given operational situation and (2) the determina- 


tion of the aerodynamic transfer functions. 

Present applications of servomechanisms to aircraft 
divide into those whose primary function is to im- 
prove the handling qualities of the manually flown air- 
craft and those in which fully automatic control is 
provided. While it is easy to envision completely 
automatic control from take-off through landing, it 
should be noted that such ‘“‘push button”’ operations 
are seriously lacking in tactical flexibility. For many 
problems, the partial use of servomechanisms in the 
form of artificial stability and control will be the bet- 
ter solution. The full implications of artificial stabil- 
ity, of which the rate damper on the yaw axis is but 
one example, are yet to be realized. By such means, 
it appears practical to eliminate the Dutch roll, the 
spiral instability, and the phugoid modes of the con- 
ventional airplane and replace them with highly 
damped subsidences. Such an aircraft—in no wise 
programmed or in any way confined in its operational 
use—may greatly influence our philosophy in the 
combining of servomechanisms with air frames. 

With regard to the aerodynamic transfer functions, 
they are vitally important, since our ultimate objec- 
tive is the control or stabilization of aerodynamic de- 
vices through the characteristics and limitations of 
aerodynamic flight controls. The still incomplete 
knowledge of these aerodynamic quantities remains 
as a stumbling block in systems analyses. It may be 
in order to supplement the author's excellent intro- 
duction to these matters by a statement of ‘‘where we 
stand”’ in this country in relation to the determination 


of aerodynamic transfer functions. 
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Response methods have been devised which deter 
mine in flight the aerodynamic transfer functions, and 
appropriate analysis schemes are available for ex- 
tracting the derivatives. Rather complete measure 
ments have been made on several moderately rigid 
aircraft and have shown the practical equivalence of 
transient and steady-state oscillation testing. <A 
variety of types of transient inputs have been tried, 
and the optimum—if such exists for a given situa- 
tion—is yet to be determined. The measurements 
have now extended on manned aircraft to high sub- 
sonic Mach Numbers, with a small amount in the 
trans- and supersonic regime, some information here 
being available from missiles. Some work is under 
way on large flexible aircraft. 

The full-scale flight methods now yield repeatable 
data, and, while a considerable degree of reliance 
may be placed upon them, they are not simple to ob- 
tain. Neither can the methods be said to have reached 
routine use. Theory and experiment by no means 
exactly agree, even in the subsonic range, and the 
evidence is that the theory may not adequately take 
into account all of the physical facts. At any rate, 
much solid research remains to tie up the loose ends. 

The question of the importance of nonstationary 
flow effects is still under much discussion and some 
investigation. The best experimental evidence indi- 
cates that these effects may not be nearly so great as 
the two-dimensional theory earlier indicated. The 
only comprehensive three-dimensional theoretical 
treatment agrees well with experiment in showing the 
minor effects of nonstationary flow to high subsonic 
Mach Numbers. By all means, let us hope that these 
effects are negligible in normal stability and control 
problems. 

Some evidence exists that the dynamic downwash 
effects for both rigid and flexible airplanes need much 
more investigation. Experimental measurements are 
required on the amplitude and phase of the downwash. 

The necessity for more research on the responses of 
aeroelastic aircraft, which the author mentions, needs 
greater emphasis. This experimental work should go 
far in orienting theoretical effort and, at the same time, 
should have immediate usefulness in connection with 
servomechanism applications. 

All of the transfer functions referred to above are for 
control surface inputs, and only recently has any 
appreciable work been undertaken in the response due 
to gust inputs. A method for determining in flight 
the frequency response in phase and amplitude to a 
gust input is required. Certainly continued research 
on the application of statistical methods to this prob 
lem is in order. 

Groundwork, leading to the stabilization of heli 
copters, which the author mentions as a future auto 
matic control application, is in fact well under way. 
An appreciable number of full-scale dynamic responses 


to impulsive inputs has been determined and is now 
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available. These should be usable for some purposes 
immediately, though they may well defy detailed 
analysis to the derivatives. 

For at least one research airplane, a sizable body of 
transfer data is available from stalled flight, and an 
automatic control is being devised to fly in the stalled 
regime. 

There is a great need at the present time for the 
development and general adoption of new  wind- 
tunnel techniques for the determination of transfer 
functions and derivatives using ordinary nondynami- 
cally similar models of good size. The flight methods 
should be used as a basis of development of dynamic 
wind-tunnel balances, so that more complete dynamic 
data can be made available in the design stages. 

Little aerodynamic information is available for 
missiles at large displacements relative to the air 
stream. Here, nonlinearities will doubtless be en- 
countered, and, as the author points out, we may well 
expect that nonlinear mechanics will come in for much 


research in the future. 


Allen E. Puckett, Head, Missile Aerodynamics Sec- 
tion, Hughes Aircraft Company: The paper presented 
by Dr. Bollay has pointed out most clearly that the 
science of controlling and stabilizing aircraft is in a 
transition phase-—the transition from manual control 
using human intelligence and muscles to various de- 
grees of automatic control using electrical and me- 
chanical brains and muscles. This comprehensive 
and detailed survey of the subject will demonstrate 
to the engineer the nature of this transition and will 
form an excellent basic catalog of the types of new 
problems which are now being encountered. More- 


~ most of the new 


over, it provides a description 
tools and techniques that are being used in the solu- 
tion of these problems, with particular emphasis on 
the vocabulary of the communications and servo engi- 
neers. This is a timely and forceful reminder that the 
aeronautical engineer must apply these methods in 
conjunction with his aerodynamic analyses. 

The most important reason for this transition can be 
described in terms of time constants. The important 
time constants in the dynamic performance of an air 
plane are various natural periods and damping times. 
As Dr. Bollay has shown, these vary inversely with 
velocity. The pilot, manually controlling an air- 
craft, can be regarded as a servomechanism with 
rather long time constants and also limited ranges of 
gain factors. At high velocities, the aerodynamic time 
constants may become short compared to the pilot 
time constants: moreover, the required range of 
control system gains may vary by factors larger than 
the pilot can accommodate. An automatic control 
system, for at least some functions, then becomes in 
evitable. 

It appears, incidentally, that an important possibil 
ity in the study of dynamic stability of manually con 








622 JOURNAL OF THE 
trolled aircraft may have been largely neglected. 
lr. Bollay observed that the study of dynamic stabil 
ity in the past has occupied a relatively academic 
position because of the artificial character of the 
fixed, 
Stated an- 
other way, the real artificiality is in the omission of 


boundary conditions imposed —controls con- 


trols free, or various other conditions. 


the pilot as an essential feedback element in the closed- 
loop system; the real “dynamic stability’’ of the air- 
craft depended upon the pilot’s ability, or lack of 
ability, to react. For example, one of the classical 
modes of oscillatory longitudinal motien is the slow 
“phugoid”’ oscillation, which is usually highly un- 
eliminates 


damped. The pilot, as a 


this oscillation so easily that he is hardly aware of it. 


servosystem, 


Several other types of ‘“‘instabilities,’’ such as the 
“Dutch roll,’ or the spiral mode, may exist without 
appreciably affecting the ‘‘flying qualities’ of an air- 
The 


definition of satisfactory dynamic stability, or ‘‘flying 


plane if the time constants are sufficiently long. 


qualities’ of an airplane, has, in the past, depended 
upon a highly empirical correlation of pilot reaction or 
opinion with the ‘open loop’’ dynamics of the air- 
craft. 

There appears to be no reason why a complete 
closed-loop stability analysis of the manually con- 
trolled airplane could not be made for at least some 
of the simple modes of motion. The pilot would be 
represented by a servosystem with particular reac- 
tions and time constants to signals such as changes in 
air speed, normal and lateral acceleration, etc. By 
making a reasonable representation of an ‘‘average”’ 
pilot, the actually dynamic characteristics of a par- 
ticular plane could be investigated, or, for a particular 
airplane, the necessary pilot gains and time constauts 
required to achieve certain specified dynamic be- 
haviors could be calculated and used to describe the 
“flying qualities’ of the airplane. This is only one 
step beyond the specification of some of the steady- 
state pilot responses, such as stick force or stick motion 
per g. 

The introduction of the automatic control system, 
particularly in pilotless aircraft, points up another 
basic change in our philosophy of aircraft control. 
The pilot was not only an essential part of the control 
system but also of the guidance system. Having 
rather clearly understood, but never too formally 
defined, capabilities, he could respond to certain types 
of information. Even modest departures from a con 
ventional presentation and processing of this informa- 
tion, as, for example, in the question of the I.L.S. 
versus the G.C.A. type of landing systems, can cause 
With the auto- 
matic control system, the whole question of what 


badly nonlinear reactions at first. 


guidance information is presented and how it is proc- 
essed becomes completely open and, moreover, be 
comes an integral part of the control and stability 


problem. Therefore, the aerodynamicist must con 
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sider the guidance system together with the contro] 
system in making a complete dynamic analysis. 

We must indeed express our congratulations and 
appreciation to Dr. Bollay for this comprehensive 
presentation of topics and techniques with which the 
aeronautical engineer must now rapidly become famil- 
iar. 


George S. Schairer, Staff Engineer, Aerodynamics 
and Power Plant, Boeing Airplane Company: Dr, 
Bollay, in this Wright Brothers Lecture, has com- 
pletely and concisely summarized the basic facts re- 
This 
is timely and important, and I would like to suggest 


lated to the servo-loop problem in aeronautics. 


that this lecture marks the official point at which the 
fundamental understanding of this knowledge be- 
comes a necessary tool of the aeronautical engineer and 
the aerodynamicist in particular. 

As one who has worked in this field over 8 years, | 
am glad to be able to tell the uninitiated that after 
an initial struggle with the subject, a fairly good under- 
standing of the problem can be reached which requires 
little further reference to the more complicated equa- 
tions. The fundamental conceptions are of vital 
importance and must be learned. 

This subject is still new, and many years will be 
required before a mature view of its relationship to 
aeronautics can be evolved. Since much of the back- 
ground material originated with the electrical engi- 
neer, we find the literature filled with electrical ter- 
minology and mathematics. We also find a strong 
tendency to solve all servo-loop problems with elec- 
trical gadgets. The next step is towards hydraulic 
The step after this appears to be 
With time, we will 
find simple solutions to most of the problems that we 
The under- 


servomechanisms. 
a return to aerodynamic servos. 


now solve with complicated mechanisms. 
standing of these solutions will not be simple, but the 


devices almost certainly will be simple. It is doubt- 





ful whether we can count on any of the complicated 
mechanisms to accomplish any large percentage of our 
first-order vitally important tasks. Complicated de- | 
vices are not that reliable yet. 

After digesting the straightforward linear theory, 
we will most always find that the actual solution of a 
problem is only started when we have solved those 
Nearly 


every problem involves a number of important non- 


problems that the simple linear theory treats. 
linearities. The most common ones are aerodynamic 
dead spots, static friction, and slop in our mechanisms. 
There are ways to get around or minimize many of 
these problems, and theoretical treatments are useful. 
However, we will spend a great deal of time on prac- 
tical cut-and-try solutions of these problems. 

Much effort in automatic stability has been di- 
rected toward the directional stability problem. The 
paper by R. J. White published in the March, 1950, 
issue of the JOURNAL OF THE AERONAUTICAL SCIENCES, 





“Tn 
XB 
tro’ 
alo: 
oth 
yea 
lar$ 
lost 
rud 
ten 
dire 
cast 
der 
dire 
trai 
stat 
dire 
Boe 
trai 
out 
ver 
dan 
tot 
The 
bee 
tior 
cil: 
mo’ 
hov 
pro 
wol 
ap] 
lati 
the 
wo! 

I 
cou 
add 
tas] 
of 1 
dur 
met 
tha 
anc 


nau 
Tec 
Ha: 
his 
pris 
sta 
esti 
tur 
cal 
wor 
Dr. 
cur 


nitro] 


and 
sive 
h the 
amil- 


Mics 

Dr. 
com- 
S re- 
This 
rest 
1 the 

be- 
and 


rs, | 
ifter 
der- 
lires 
yua- 
vital 


| be 
D to 
ack- 
ngi- 
ter- 
‘ong 
lec- 
ulic 
» be 
will 
we 
ler- 
the 
ibt- 
ted 
our 
de- 


di- 
he 
50, 











AERODYNAMIC STABILITY 


“Investigation of Lateral Dynamic Stability in the 
XB-47 Airplane,’ presents the problem and an elec- 
tromechanical solution to the problem. The effort 
along this line has been continued with studies of 
other solutions to the same problem. A number of 
years ago aerodynamicists became aware of the 
large amount of directional stability which was being 
lost with trailing rudders. This resulted in the use of 
rudders with small trailing tendencies or even leading 
tendencies. This was followed by a bad epidemic of 
directional oscillations. 
cases was caused by rudder movement of leading rud- 


A large percentage of these 
ders. A recent attempt has been made to solve the 
directional damping problem by the use of a strongly 
trailing rudder, wherein the rudder moves to reduce the 
static stability but increases the damping of the 
directional oscillation. Such a rudder was tried on a 
Boeing 377 stratocruiser. 
trail strongly through neutral and to thereby wash 
out about half the basic directional stability of the 
This trailing action was sufficiently 


The rudder was made to 


vertical tail. 
damped by aerodynamic forces to give a phase shift 
to the rudder motion and thus an increase in damping. 
The response of the airplane to gusts appears to have 
been materially reduced by the reduction in direc- 
tional stability, and the damping of the ensuing os- 
cillation was greatly increased. The period of the 
motion was also increased. This is an example of 
how the basic servo-loop theory can be used to solve a 
problem by simple means. Before leaving this, it 
would be well to mention that this solution does not 
appear to be a general solution to the directional oscil- 
lation problem and will not work in many cases. The 
theory given by Dr. Bollay can be used to separate the 
workable cases from the unworkable ones. 

I hope that the aerodynamic educators of this 
country will take note that this subject has been 
added to the required list. This will not make their 
task any easier but is necessary just the same. Many 
of us were students of Prof. Stark Draper at M.I.T. 
during the period when he was first adding the funda- 
mentals of this subject to his courses. We all feel 
that his leadership in this field has been outstanding 
and should be recognized. 


R. C. Seamans, Jr., Professor, Department of Aero- 
nautical Engineering, Massachusetts Institute of 
Technology: In the late thirties, Dr. Bollay was at 
Harvard University, and I was fortunate to be one of 
his students. Therefore I am not in the least sur- 
prised to find that his discussion of aerodynamic 
stability and automatic control is up to date, inter- 
I feel that this lec- 
ture will serve as a standard reference for aeronauti- 
cal engineers in the field of stability and control. | 
would like to take this opportunity to congratulate 
Dr. Bollay and to make several observations that oc 


esting, and especially thorough. 


curred to me while listening to his address. 
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In this paper, the transformation into the complex 


Vip) = 


domain is carried out by use of equation 
pf e ” V(t) dt. This transformation is equal to p 
0 


times the conventional Laplace transform used in 
standard textbooks. Dr. Bollay includes the p so 
that the transformed quantity will have the same di- 
mensions as the original quantity in the time domain. 
I feel that this argument hardly justifies the added 
complication of the extra term and the confusion 
that will arise out of the lack of uniformity between 
the transforms listed in this paper and those found in 
most textbooks. 

It is stated in this paper that control engineers pre- 
fer an airplane that is approximately neutrally stable. 
The advantage of such an airplane is that it can be 
made maneuverable even with a low elevator effective- 
ness, since its static stability is low. However, if a 
small elevator displacement gives a large angle of at- 
tack, it is impossible to limit a maneuver with stops 
on the elevator servo drive. Automatic limiting is 
required as a function of either angle of attack or nor- 
mal acceleration. It is felt that this problem is suf- 
ficiently important to warrant thorough study by some 
research agency. 

By means of dimensional analysis, Dr. Bollay de- 
velops a relationship between the natural frequency of 
a missile and an airplane. He concludes that the 
missile has 100 times as large a natural frequency in 
the longitudinal mode. This figure, although pos- 
sible, is not representative. 
longitudinal responses, the important linear dimen- 
sion is the length, not the diameter. In addition, the 
speed ratio is larger than is ordinarily encountered 
For this reason it is felt that a natural frequency ratio 
of between 10 and 20 is more likely to be experienced. 

Dr. Bollay outlined the four steps required when 


When comparing the 


developing a control system for an airplane. These 
four steps are: (1) theoretical design, (2) simulation, 
(3) flight test, and (4) refinements based on funda- 
mental noise studies. 

These four phases overlap one another on a chrono- 
logical basis. For example, asa resuit of flight test, the 
dynamic performance of an airplane may be found 
different from that predicted by theoretical and wind- 
tunnel considerations. When this is the case, the 
theoretical design of the control system must be modi- 
fied, and consequently step (1) may depend on step 
(3}. 

Also, when the automatic control equipment is 
undergoing flight test, it is highly desirable to check 
the flight results against simulator results so that 
valuable flight time is not wasted. One highly satis 
factory procedure involves the use of a high-speed 
electronic simulator. This simulator is a combination 
of the straightforward type and the electronic inte 
grating type mentioned in Dr. Bollay’s paper. It 
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A Limiting Case for Missile Rolling Moments 


ERNEST W. GRAHAM* 
Douglas Aircraft Company, Inc. 


SUMMAR\ 


Missile rolling moments are studied by linearized theory for the 
limiting case in which infinitely many wings are symmetrically 
arranged around a circular cylinder. 

This case is physically interesting when interference effects 
prevent the rolling moments from getting too far out of the range 
of those practically obtainable with a finite number of wings. 
Mathematically, the limiting case is convenient because it can be 
analyzed through knowledge of the rate at which angular momen 
tum is being added to the air, the rolling moment being propor 
tional to this quantity. 

The analysis is independent of Mach Number, and, since there 
are few restrictions on the geometry of the system, it should 
furnish an upper limit for rolling moments for a variety of con- 
figurations and flight speeds. Many cases with a finite number 
of wings are difficult to analyze, and knowledge of a limiting value 
may help in checking intermediate cases. 

For a delta wing system with leading edges on the Mach lines, 
body of zero radius, and each wing acting as an aileron, the roll- 
ing moment coefficient for an infinite number of wings is w times 
that for the planar wing and 1.85 times that for the cruciform 
wing. For narrow delta systems, the above values apply ex 
actly for the planar case and approximately for the cruciform 

Wing-tail interference effects on aileron rolling moments may 
also be studied for the limiting case. Here, all angular momen- 
tum is removed from the disc corresponding to the maximum cross 
section of the tail (unless opposite surfaces of the tail are designed 
to ‘‘float’’ relative to each other and are controlled only through 
the sum of their displacements). If the tail disc is as large as the 
wing disc, then no resultant rolling moment is produced by wing 
ailerons in the limiting case. This assumes the dises to be con 
centric. For a lifting wing, the effective displacement of the disc 
centers might introduce important nonlinear effects. 

Damping in roll can be studied by assuming the wing angle of 
attack proportional to the radius. Equating the damping mo 
ment to the aileron moment gives an equilibrium rate of roll which 
may be expressed as the ratio of tip helix angle to aileron dis 
placement angle. For the delta wing system previously con- 
sidered (with sonic leading edges ), this ratio decreases as the num 
ber of wings increases. For an infinite number of wings, the 
ratio is two-thirds of the planar value or 0.78 times the cruciform 
value. For a given aileron displacement, the equilibrium rate of 
roll is therefore decreased by increasing the number of wings. 
Physically, this corresponds to the fact that interference effects 
on damping in roll are less than on aileron rolling moments 
The above figures are approximately correct for narrow delta sys 


tems 


INTRODUCTION 


Ir MISSILE DESIGN, the wing surfaces do not necessarily 
lie in a single plane but may, for example, be given 
a cruciform arrangement. The analysis of such sys 
tems is more difficult than the analysis of planar sys 


Presented at the Supersonic Aerodynamics Session, Nineteenth 
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1951. 
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tems, and for this reason it is desirable to consider a 
limiting case in which infinitely many wings are sym- 
metrically arranged around the body axis. This case 
is of physical interest because interference effects may 
limit the rolling moments to values not too greatly ex 
ceeding those for the cruciform narrow delta wings. 


DEVELOPMENT* 


Procedure 


The rolling moment acting on a wing is ordinarily de- 
termined by finding the spanwise lift distribution and 
integrating the rolling moment contributions of all the 
differential elements of span. Alternatively, the rolling 
moment acting on the wing may be obtained from a 
knowledge of the rate at which angular momentum is 


being added to the air, as shown in Eq. (1). 
Mp = / | puV gr? dé dr (1 


Mp, = rolling moment 


where 


p = mass density of the air 

u = axial velocity of the air 

V, = tangential velocity of the air about the roll- 
ing axis 

r = radius out from the rolling axis 

8 = angular displacement about the rolling axis 


The integral is taken in a plane behind the wing and nor- 
mal to the flow direction. Fig. 1 shows the coordinate 
system and type of problem being considered. 

Eq. (1) can be obtained by applying a general mo- 
mentum theorem! stating that in a steady flow the 
moment exerted on the surfaces bounding a region fixed 
in space is equal to the net rate of flow of angular mo- 
mentum out of the region. If the region is chosen to 
be that between two planes normal to the flow direc- 
tion, one ahead of and one behind the wing, then Eq. 
(1) is obtained, since no net angular momentum flows 
through the upstream plane and the only moments 
acting on the bounding surfaces are those contributed 
by the wing itself. 

Eq. (1) may be written as 


Mr = J [ (po + Ap) (to 4 
0 0 


* Valuable contributions have been made by Drs. P. A 
Lagerstrom and Z. O. Bleviss, of the California Institute of Tech 


Au) Vor? dé dr (2 


nology, and by Harold Luskin, Aerodynamics Research Engineer, 
Douglas Aircraft Company, Inc., Santa Monica Plant 
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where Pp = po + Ap, “a= & } 


indicates free-stream conditions. 
Effect of Linearization 


For a flow in which the velocity perturbations are 
small, the equation may be linearized to give 


Mr = pollo / / V,r° dé dr (3) 


2 


This can be further simplified, since 7 J V, dé 


Jo 
taken outside the wing tip circle is the circulation 
around a closed curve in the fluid. This circulation is 
zero, since the flow was originally irrotational and no 
If R = 
tip radius, then (as indicated in Fig. 2) Eq. (3) can be 
rewritten as 


R *2a 
Mr = polly [ / Vr? dé dr (4) 


Introduction of an Infinite Number of Wings 


viscous effects have been encountered. wing 


If a planar wing with deflected ailerons is considered, 
then the value of tangential velocity (V4) is known at the 





THE SMOOTHING OUT 
OF THE TANGENTIAL VELOCITY PATTERN 
THROUGH THE ADDITION OF WINGS 
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wing surfaces within linearized theory, since I, = 
au, Where a is the local angle of incidence. However, 
this determines |’, only at 6 = 0, 7, and 27 as shown 
in Fig. 3, and, presumably, I’, is less at intermediate 
angles as shown. For a cruciform wing with four 
ailerons deflected, V,is known at 6 = 0, 7/2, 7, 37/2, and 
27, and, if the same chordwise dimensions are retained as 
on the planar system, the fluctuations in |’, should be 
decreased between the wings. Since this still does not 
provide sufficient information to permit solving Eq. 
(4) for the rolling moment, it is necessary to consider 
the limiting case in which infinitely many wings are 
arranged about the axis so that adjacent wings are 
separated by an infinitesimal distance. Here, the 
fluctuation in tangential velocity between two wings is 
negligible (assuming finite chordwise dimensions), and 
therefore the tangential velocity distribution is known 
Only the tan 
gential velocities at the trailing edges of the wing sys 
Behind the trailing edges 


everywhere inside the wing tip circle. 


tem are important. 
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remains constant for a given radius, since the circulation 
is constant for that radius. Therefore, the rolling mo- 
ment (.W/,) can be evaluated from Eq. (4) using the 
trailing-edge values of V,. The trailing edges are de- 
fined as all edges that are the last encountered by air 


leaving the wing system. 


The above method of analysis is related to methods 


used in propeller theory.* 


DISCUSSION 


Limitations 

The preceding analysis of the limiting case is inde- 
pendent of Mach Number and independent of the longi- 
tudinal shape of the wing system. It is only necessary 
that, for a given finite Mach Number, the ratio of the 
distance between adjacent wing tips to chordwise and 
spanwise dimensions of the wing system approach zero. 
For a system with a finite number of wings, such dimen- 
sionless ratios give a crude qualitative indication of how 
closely the limiting conditions are approached. 
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Modifications and Extensions 


The presence of a circular cylinder in the system 1s 





accounted for by modifying the inner limit of integra- 
tion. Other bodies of revolution produce a radial dis- 
placement of the stream lines, which should be con- 
sidered. | 
Damping moments in roll can be calculated by as- 
suming the wing angle of attack proportional to the } 
radius, and equilibrium rates of roll are then obtained 
by equating damping moments to aileron moments 
When the equilibrium rate of roll is attained, the net 
angular momentum left behind the wing system is zero. 
It is also possible to determine rolling moment inter- 
ference effects between wing and tail if one assumes that 
the tail also has an infinite number of symmetrically 
arranged surfaces. This is actually a special case of a 
wing system having an irregular longitudinal shape 
In general, the tail (assumed to be in the rear) merely 
removes all angular momentum from the disc that tt 
covers, and, if it is as large as the aileron disc and con- 
centric with it, no resultant rolling moment is produced 


by wing ailerons. However, the tail does not affect the 
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A LIMITING CASE FOR 
damping moment in roll* (unless its diameter is greater 
than the wing disc diameter). This follows because the 
air flowing through the tail system has already been 
given an angular velocity by the rolling of the wing 
and, hence, has no angle of attack relative to the tail 


surfaces. 


Nonlinear Effects 

If a lifting wing is considered so that the wing wake 
is displaced relative to the fuselage axis, then, in effect, 
the tail and wing discs are not concentric and some of 
the angular momentum of the wing wake may not be 
removed by the tail even though the discs have the same 
diameter. This is a nonlinear effect that may be im- 
portant* and tends to reduce the unfavorable rolling 
moment produced on the tail for the case where the 
wing and tail discs are of equal diameter. For the 
case where the tail disc is of smaller diameter than the 
wing disc, it is not necessarily true that the nonlinear 
effects on rolling moment are favorable. The effective 
displacement of the dise centers may also introduce 


nonlinear yawing moments. 


* Suggested by Dr. Lagerstrom. 
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RESULTS AND COMPARISON 


Table 1 shows some comparative values of C,, and 
C,, and equilibrium rate of roll for delta wing systems 
having two wings, four wings, and an infinite number 
of wings. The values are obtained from Examples (1) 
and (2) of this paper and from the references listed. 


EXAMPLE (1). THe LimitiInG CASE OF INFINITELY 
MANY AILERONS DEFLECTED IN THE SAME DIRECTION 


Assume that each surface in Fig. 4 is an aileron de- 
flected through an angle 6. Then, the tangential ve 
locity (V4) = 6m and, from Eq. (4), 


"R *2n 
Mr = polly J / bur? dé dr 


Mr = potty?6 (27R*/3) 


or 


letting 
q = (1/2)potto?, Mp = 5g(4/3)rR? 


The rolling moment coefficient (C;,) defined as 
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TABLE |. Rolling moment coefficients and equilibrium rates of rell for delta wing systems 
4 < 
“\ / 
IYPE ws w / > 
OF le Uo * \. Uo “i Uy Pi \ 
ack —_ -—- ‘ —~_ - on od ~ 
W ] N G ¥ me 
SYSTEM : \ 
~ \ 
— f -——3- 0 1.0 co 
tan W 
cc. /ren 2 WINGS 0, 66719) 6. 66707) 0. 667/£?) 
’ ‘ (3) 132!3) 1.016/¢°7 
‘ WINGS 1. .O1€ 
AILERON 1 WINGS 1. 333 | 
EFFECT co WINGS co 2.094 2.094/f 
FC, /ranp| 2 WINGS 0. 333(°) 0. 33309) 0. 393 /¢°) 
4 7 
DAMPING |4 WINGS 0. 667°4) 0. 623! ) 0. 637/£' ) 
INROLL leo WINGS oO 1.571 1. 571/£ 
TIP HELIX] |2 wiINnGs 2. 000 2. 000 1.697 
ANGLE 
AILERON 4 WINGS 2. 000 1.817 1.595 
ANGLE . 
EQUIL. ROLL |eO WINGS 2. 000—1. 333 1. 333 1. 333 
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Fic. 9. 


Ci, = Mr/q(cR)b6 


becomes 


Ci, = (27/3) tan w 


This is plotted in Fig. 5 to compare with values for 
plane delta and cruciform delta ailerons determined by 
Bleviss.* The cruciform configuration considered by 
Bleviss has a fixed vertical fin and horizontal ailerons. 
However, by rotation and superposition, it is easily 
shown that C;, values for horizontal ailerons in the 
presence of a fixed vertical fin are doubled by deflecting 
the vertical surfaces through the same angle. 
EXAMPLE (2). DAMPING IN ROLL FOR THE LIMITING 
CASE OF INFINITELY MANY WINGS 


Assume that the angle of attack of each wing in Fig. 
6 varies linearly from zero at the axis to —pr/u at the 
tip, with p = rolling velocity “ radians per second, 


Then, V, = —pr and, from Eq. ( 
Mp = —potto f te prr? dé dr 

or 

Mr = —(gp Uy) rR! 
and 

C, = —(pb/2u) (1/2) tan w 

or 

Ci, = —(m/2) tan w 


This is the rolling moment coefficient due to damp- 
ing. If the equilibrium rate of roll resulting from a 
given aileron deflection is desired (assuming that the 
complete wings act as ailerons), the sum of the rolling 
moments due to aileron deflection and damping is set 


equal to zero, giving 


(pb 2u) = (4 3)6 
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Figs. 7 and 8 show the effect of number of wings op 
damping in roll and equilibrium rate of roll. In Fig, § 
it is seen that a singular point occurs at f = 0 for an 
infinite number of wings. If f is set equal to zero and 
the singularity is approached by starting with a finite 
number of wings and increasing the number without 
limit, then the value of 2.0 results for the ratio of tip 
alternatively, the 


helix angle to aileron angle. If, 


singularity is approached by starting with an infinite 


number of wings and letting f approach zero, then a 
value of 4/3 results. 
EXAMPLE (3). THE EFFECT OF A CENTRAL CIRCULAR 


CYLINDER 


If in Example (1) the inner limit of integration is 
d/2, corresponding to the presence of a central circular 
cylinder of diameter d (as shown in Fig. 9), then the roll- 
ing moment coefficient produced by an aileron deflection 
6 becomes 


C, = (27/3)6 tan w[1 — (d/b)*] 


Similarly, in Example (2), if the inner limit of integra- 
tion is d/2, then the rolling moment coefficient pro- 
duced by damping in roll becomes 


; (2) 7 tan “| (; 
4 ,; = | —_ ) 
219 2 b 


Plots of the ratio of C, to (C))4g=0 are shown in 
Fig. 9. The coefficients are based on the span and area 


corresponding to (d/b) = 0. 
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Interference Between Wing and Body at 
Supersonic Speeds—Theoretical and 
Experimental Determination of Pressures 
on the Body 


R. H.C 


RAMERt 


Cornell Aeronautical Laboratory, Inc. 


ABSTRACT 

Experimental data are adduced to show the pressure dis- 
tribution on an axially symmetric body, aft of a lifting wing, 
in supersonic flow. The axial distribution of integrated lifts, 
as well as the peripheral pressure plots, is shown to be predicted, 
to engineering accuracy, by use of the linearized theory pro- 
pounded by Ferrari, even when only a few terms are used in the 
Fourier development required by the theory. The theoretical 
curves were obtained by means of a new simplified solution to 
Ferrari's differential equation; it is expressed in terms of the 
hyperbolic functions of successive integral multiples of an argu- 
ment, so that the work of numerical computation progresses 
rapidly with the aid of standard tables. 


INTRODUCTION 


— THEORY OF WING-BODY INTERFERENCE, pro- 
pounded by Professor Ferrari in reference 1, 
presents a method based on the familiar linearization 
of the equation of motion for calculation of the aero- 
dynamic action of the body upon the wing and also 
that of the wing upon the body. Experimental pres- 
sure distributions are now available to compare the 
results of theory and experiment for the induced pres- 
sures on the body (the results for the body induced 
effects on the wing have already been noted in reference 
2). These data were obtained on a wing-body com- 
bination that closely reproduced in the wind tunnel the 
features of the mathematical model treated by Pro- 
fessor Ferrari in the numerical application of his theory. 


DISTRIBUTION OF INCREMENT OF LirT ALONG Bopy 


It was demonstrated in reference | that the distribu- 
tion of the incremental lift along the body which arises 
from the normal component of the induced velocity 
(which must be negated at the body surface by the new 
distribution of doublets along the body axis) has, for 
the configuration shown in Fig. 1, a form and magnitude 
such as given in Fig. 2. To this distribution must be 
added the effects of the longitudinal components of the 
induced flow from the wing (a uniform, two-dimen- 

Received December 29, 1950. This investigation, both ex- 
perimental and theoretical, was carried out under Navy BuOrd 
sponsorship 


t Principal Research Engineer. 


sional velocity distribution over a rectangular wing at 
an angle of attack, ®», is being assumed at this first 
stage of the iteration process; the process is described 
in detail in reference 2). The way these two distribu- 
tions are related and summed is shown in Fig. 3. 

A comparison of this total increment in lift along the 
body with the experimental result is given in Fig. 4 
for 8° angle of attack. Even though the distribution 
of lift along the body is predicted to good accuracy, one 
still has no assurance that the actual pressure distribu- 
tion around the periphery of the body is represented 
as Closely by the theory as is the lift. 


LocAL PRESSURE DISTRIBUTIONS AROUND Bopy 
BEHIND WING 


A representative comparison of theory and experi- 
ment is presented here as Fig. 5 for the pressures that 
exist around the periphery of the body at an axial 
station denoted by &* = 1.15. <A section through the 
body at the leading edge of the wing has an axial loca- 
tion given by &* = 1.0, where, in general, §&* = x/A 
i.e., the nose length depicted in Fig. | is denoted by X. 
Thus this particular section selected for illustration 
of a typical peripheral pressure distribution is located 
slightly behind the trailing edge of the wing whose 
axial coordinate is * = 1.14. The physical conditions 
of the airflow in the wind tunnel have been previously 
listed in the above-cited paper? by Ferrari. The 
reason for illustrating these results by selection of data 
for a cross section lying aft of, but near, the wing is 
merely in order to obtain large values of the increments 

















= 
Fic. 1. Configuration used in the numerical application and 
upon which the wind-tunnel tests were performed 
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in pressure due to interference; in this way the con- must be made to the formulas in order to obtain th 


tribution to the pressures that arise from the flow correct values for the individual pressures. It was 
about the body alone are rendered relatively inconse- noted over a year ago by Lagerstrom and Van Dyke 
quential (although they are actually accounted for in in a footnote to reference 3, that an operator had beep 


the theoretical distribution plotted in Fig. 5). applied incorrectly in obtaining the formulas for the 

; 7 » individual pressures according to reference 1. Re. 
ORRESPONDING THEORY FOR THE PRESSURES ; , 

C . putes cently, Bolz! has worked out the details of the firs 


The determination of the theoretical results requires two terms; it will require several more terms, however. 
some note of explanation. Although the formulas’ in order to attain a practical description of the pres. 
given by Ferrari in reference | are valid for calculation sures around the body. Therefore, the continuatioy 
of the lift distribution as given by Fig. 2, a modification of the suggested method is apt to become unwieldly, 


SIMPLIFICATION OF FERRARI'S FORMULAS FOR BODY PRESSURES 
Fortunately, the whole problem can be obviated by noting that the operator procedure is not actually neces. 
sary; an identical solution can be obtained by use of hyperbolic functions as follows: 
The differential equation to be solved {see Eq. (4) of Professor Ferrari's original paper] is 
0*¢,,‘” 4 om + | O¢,,‘*® -0°o,** 


= @& 


or? r Or Ox? 


The required solution of this equation can be shown by direct substitution to be 


m m 0 
. (—1)"a : 
Pn” = fm (x — ar cosh yu) cosh mu du (2 


m 
r arcosh ¥ /ar 


This basic result appears to have escaped the attention of Lamb® and later researchers. 
Ward,* however, has presented another simple solution, but the subsequent evaluation of the arbitrary function 
in terms of the body coordinates would involve the Laplace transform, which may be inconvenient to evaluate. 


It follows immediately from Eq. (2) that 


O(r" bn) a ; = (= | 
= (—1)"+lqm+l fm(x — a@R cosh u) cosh mu cosh u du = —F.,, (= 3 


or r=R arcosh x/aR 


Then, by employing the same technique of division of the interval of integration in small subintervals, through 





out which the individual /,, are assumed constant (as was illustrated in Ferrari's original paper), it follows that 


ear ow 4 x x mx, . xy . JB ' 
r Rig \ — 1 cosh m arcosh — sinh m arcosh = F,, t) | 
m* — | a’R,? ak, aR, ak, R, 
for m # 1, while the formula previously given by Ferrari in the case of m = | is still valid. The A,,") employed 


in Eq. (4) are, of course, the constant values of the /,, in the first small interval; while, in general, if x;—;' and 
x,’ are the abscissas of the left and right ends of the ith interval, respectively, and K,,” is the constant value o!} 
/, in that interval, then one has, at the point P’”’, with coordinates (x,, R,), that 


, 


( ede ] \*a™ 4-1 n x; — x 1 


a? S. = ! 1 \2 Xn a | 
. > Ky \ —1 cosh m arcosh 
mem L oss aR, ak, 





od te Xn — Xi M(X_, — X-1') Xe — Mir 
— lcosh m arcosh - — sinh m arcosh + 
aR, aR, aR, aR, 
/ / 

M(Xn — Xi) . Xn — Xj <i : 
sinh m arcosh =F, d 

aR, aR, Re 

for m # 1, while the formula previously given by Ferrari in the case of m = | ts still valid. 
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INTERFERENCE 


CONFRONTATION OF EXPERIMENTAL AND THEORETICAL 
RESULTS 


Calculations based on Eq. (5) have been carried out 
in order to obtain the theoretical values for the pressure 
coefficients pertaining to the wing-body configuration 
depicted in Fig. 1 and more fully described in reference 
Rs 

The correspondence between theory and experiment 
shown in Fig. 5 is considered to be good, especially in 
view of the fact that only an m = 7 representation 
(actually only four nonzero terms computed) is em- 
ployed in the determination of the effect of the normal 
components of the induced velocities in producing the 
graphed pressures; in addition, these four terms are 
the opening sequence of only the first step in the 
iteration process mentioned previously.” 

Further steps in the iteration process would be 
necessary to give a more precise determination of the 
flow along the side of the body behind the leading edge 
of the wing, since the secorid stage of the iteration 
process for the distribution? of velocities over the wing 
shows marked changes from the two-dimensional dis- 
tribution being assumed to exist over the wing when 
carrying out this first stage of the iteration process for 
the body. Nevertheless, it is highly significant that 
close correspondence is obtained with the experimental 
results even without recourse to such further mathe- 
matical refinements provided by Ferrari’s theory; the 
more complicated formulas required for carrying out 
the second and higher numbered steps of the iteration 
process are available in reference 7. Note, however, 
that the sequence of steps in the iteration process ad- 
vocated in reference 2 is reversed in reference 7; the 
sanction for this change is found in the smallness of the 
subsequently computed second-order effects. The ex- 
planation of why, in the first step of treating the body 
upwash effect on the wing, a nil amount of upwash is 
assumed to exist in the region of the wing actually 
covered by the body acquires an acceptable logical 
foundation on the basis of the more extensive argu- 
ments advanced in this latter reference. 
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of attack (Mach Number = 2.0). 


PRESSURE “‘JUMPS’? DEMARCATING SHOCK SURFACES 
FROM WING LEADING AND TRAILING EDGES 


There will, of course, always tend to be some dis- 
crepancy between the experimental distribution and 
the theory in the neighborhood of the jumps in the 
distribution which occur at the edges of the region 
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Fic. 5. Final determination of pressure coefficient obtained 
by adding to distribution resulting from the normal components 
of the induced velocities that part of the local pressure contributed 
by the body alone, plus increment to pressures arising from 
longitudinal components of induced velocity at Station P® where 


£* = 1.15 (2.25 in behind wing leading edge on model). Mach 


Number = 2.0; 8° angle of attack. 


influenced by the longitudinal components of the 
induced velocities. Because of the ameliorating in- 
fluence of the boundary layer it can be expected that 
such sharp discontinuities will not, however, be ex- 
perienced in physical reality. An experimental value 
of the pressure at such a jump point which turns out 
to be approximately the average of the values for the 
superior and inferior limits of the theoretical curve on 
either side of the jump can be considered as proper 
confirmation of the predicted pressure in such a 


location. 
AVAILABILITY OF ADDITIONAL DATA, INCLUDING WING 
PRESSURES 


These data are but a small part of the results ob- 
tained in the wind-tunnel tests on this configuration. 


It is believed, however, that they are suflicient to 
show the essential confirmation of Ferrari's theory ag 
it relates to pressures on the body. Confirmation of 
the theory with regard to pressures on the wing was 
presented in reference 2. Additional wing data over 
and above those presented there are also available 
although there is little need to further substantiat, 
the evident relationships by means of similar plots 
These data include pressure measurements on the aff 
region of the wings, close to the body, where the lack 
of information has been deplored by Morikawa.’ Thy 
complete set® of experimental data and _ correlations 
with theory may be obtained on request from The 
Johns Hopkins University, Applied Physics Laboratory 
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RIEF REPORTS of investigations in the aeronautical sciences 
and discussions of papers published in the JOURNAL are 
presented in this special department. The publication will be 
comple ted approximately 8 fo 10 weeks after receipt of the material 
The Editorial Committee does not hold itself responsible for the 
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Lateral Elastic Instability of Hat-Section 
Stringers 


C. M. Tyler, Jr., and T. J. Walker 

Department of Mechanical Engineering, Carnegie Institute of 
Technology, Pittsburgh 

June 11, 1951 


— AND SHEET COMBINATIONS under longitudinal com 
pression loading may fail in several ways. 


fail as columns by buckling perpendicular to the plane of the 


The stringers can 
sheet; they can fail due to local crushing or crippling; or they 
can fail by the edge of stringer, which is not attached to sheet 
buckling parallel to the sheet. This last mode of failure is known 
as lateral instability and commonly occurs with open section 
stringers. 

It is generally believed that closed stringers such as hat sections 
are unlikely to fail due to lateral instability. As a result of some 
tests in which this type of failure apparently occurred, Goodman! 
has developed two energy methods for calculating the lowest critical 
strain corresponding to this mode of failure. The simpler method 
involves two arbitrary constants and assumes that the legs of the 
hat are fixed at the sheet. The second inherently more accurate 
and more complicated method involves five arbitrary constants 
and assumes that the legs are hinged at the sheet. The second 
method gives appreciably lower numerical values 

The writers felt that the difference in the numerical results of the 
two methods was caused primarily by changing from fixed-leg to 
hinged-leg conditions. Consequently, they adapted Goodman’s 
two-constant method to hinged legs and have found that numeri- 
cally it checks closely the five-constant method. Thus, it may be 
concluded that the accuracy of this simpler method is essentially 


the same as the five-constant method. 
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This two-constant hinged-leg method has been put into dimen 
sionless form and used to calculate Fig. 1. These curves indicate 
that lateral elastic instability is not likely to occur unless the 
heigiit-width ratio is above 5; they also indicate that the critical 
stress can be rather low for large values of this ratio 

For height-width ratios above 5, Fig. 1 shows that the critical 
stress is rather insensitive to thickness. Neglecting thickness, a 


simplified conservative critical strain formula can then be used. 


a l \ l l | 
E 2D + 1<2.6 n? l l 
} 6.76 (1 + 2R) + + 2 \ 
12S? 5.2\R 
(1) 
where 
8 (C2 = 
D = = + — CR? +—-R 
105 R38 30 ‘ 
yy (1/R) + [1/(11 + 2R)] 
R h/w, height-width ratio 
S = L/w, length-width ratio 
Eq. (1) is applicable to hat sections with parallel legs and is 


based on the two-constant hinged-leg method. A lower bound 
expression can be found by taking the length as infinite 

The possibility of lateral elastic instability being the indicated 
mode of failure can be checked with the aid of Fig. 1, 
be, the critical stress can be conservatively calculated by Eq. (1) 


and, if need 
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Reflection of Strong Shock at Nearly Glancing 
Incidence* 


Donald R. White 

Charles A. Coffin Fellow, Palmer Physical Laboratory, Princeton 
University 

June 7, 1951 


ING AND LUDLOFF! HAVE OBTAINED a solution for the density 
yp pressure fields in the reflection of shock waves from a 
wedge of small angle (nearly glancing incidence) and have com- 
puted explicit results for a shock of strength § = P;/P: = 0.137 
which sets up a flow of Mach Number 1.21. During the past year 
the shock reflection problem? has been the subject of an experi 
mental survey at this laboratory utilizing a Mach-Zehnder inter- 
ferometer in conjunction with the shock tube* to determine the 
entire density field for a range of shock strengths and wedge 
angles. In the course of this program results have been obtained 
which may be compared with their predictions as presented in 
Fig. 1 of reference 1. Features of their results to be noted are the 
region of uniform flow near the leading edge of the wedge, the 
singular point on the reflected shock where the contours of con- 
stant density converge, and the separation of the density and 
pressure contours in the region where a slipstreamT would appear 
if the reflected shock were not of zero strength at the point where 
it meets the incident shock. This separation is also implied by 
the previous result of Lighthill' that the Mach stem for strong 
shocks has a sharp maximum in curvature near the triple point, re 
sulting in different entropies for adjacent streamlines 


* This research was supported by the Office of Naval Research of the 


Navy 
t The term 


Department 
as used here refers to the density and velocity 
three 


“‘slipstream”’ 


discontinuity surface necessarily resulting from the intersection of 


shocks, all of nonzero strength 
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Our experimental results for a shock of this strength and a 
wedge half-angle « = 5.4° are shown in Fig. 1. This angle is such 
that the shock is not quite attached, so that one does not observe a 
region of uniform flow near the leading edge. The density con- 
tours do seem to converge toward a point just beyond the re- 
flected shock, and their behavior near the expected position of the 
Fig. 2 represents an attempt to obtain a 
To do this it 


slipstream is similar. 
picture qualitatively similar near the leading edge. 
was necessary to use a stronger shock and a smaller angle, but the 
same general features again appear. 

This convergence of the density contours appears to be charac- 
teristic of the pattern found when a weak shock is intercepted by a 
rarefaction and has been observed in several otherwise unrelated 
experiments, such as the regular reflection of a weak shock and 
the diffraction of a weak shock over a corner. In the present case 
the shock is the attached bow wave and the rarefaction is the 
signal sent out from that portion of the fluid which was at the 
corner when the shock first hit the wedge. It is observed, in 
general, that for a shock of small but finite strength the apparent 
point of convergence lies just beyond the shock, approaching it as 
the shock strength decreases. Thus, the singularity appears only 
in the limit of an infinitely weak shock, where the magnitude of 
the disturbance vanishes. 

The similarity of the density contours near the slipstream is 
striking, but too much significance should not be attached to this 
point, since the duration of the spark used in taking the interfero- 
grams tends to smooth out fringe discontinuities so that the 
physically sharp slipstream that we would expect for a finite angle 
wedge would appear as a region of continuous but rapid density 
change. The maximum shift of the fringes in the interferogram 
from the undisturbed flow is about one fringe, so the error in the 
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numbers attached to the contours is of the order of the chang 
from one contour to the next, but the shape of the contours i 
relatively accurate. 
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Ring Shape Source Distribution in Axially 
Symmetrical Supersonic Flow 


S. |. Pai 
Institute for Fluid Dynamics and Applied Mathematics, University 
of Maryland, College Park, Md 


June 14, 1951 


oo LINEARIZED THEORY Of axially symmetric supersonic floy 
can be used to study the supersonic flow past a tube of nearly 
constant radius for both the external and the internal flow con 
ditions. For simplicity we consider here the case without angle of 
attack. The fundamental equation is 

' 0*o O° 1 Op 


(Af? ) } 
“one or? r Or 


where J/ is the Mach Number of the basic supersonic flow, ¢ is the 
velocity potential of the disturbance, x is the coordinate along the 
axis of the tube, and 7 is the radial distance from the axis of th 
tube. 

The most general solution of Eq. (1) is! 


= i f(x — Brcosh u) du 4+ J F(x + Brecosh u)du (2 
0 


where 8? = 1/* — 1 and both f and Fare arbitrary functions to be 
determined for a given problem so that certain boundary values 
are satisfied. 

In case of external flow past a tube, the condition of no incom 


ing disturbance gives 
eo 
o= J f(x — Br cosh u) du (3 


The solution of the problem is to find the source distributio! 


f(x) on the x-axis (ry = 0) such that the boundary condition on the 


average radius of the tube is satisfied. 
In case of internal flow, it is often erroneously assumed that the 
solution of the problem is 


(® «x 
= i F(x + Br cosh u) du 4 
( 


Because the Mach lines of Eq. (4) are x + Br = constant, it looks 
like an incoming wave. However, if one investigates Eq. (4 

little in detail, it is easily seen that Eq. (4) does not represent the 
incoming wave at all. The proof is as follows: 


If one uses the transformation 
x + Brcoshu =é 
Eq. (4) becomes 
F() dé 
eo = f 
xt+er V (é — x)? — 8*r? 
where F(~) represents a source distribution on the x-axis. Fof 
simplicity of illustration, we consider only a simple source situated 
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at the origin —i.« 0, F(é) #0. It is easily seen that 
the integral (6) is different from zero for all the points on the left 


of the surface x + Br = 0 and equal to zero for all the points on 


, only at x = 


The region of influence is upstream! 
flow. 


the right of the surface. 
(See Fig. 1. It 
Hence, Eq. (4) or (6) cannot be used to represent the flow inside 


contradicts the behavior of supersonic 


Furthermore, it artificially introduces source distribu- 
y-axis inside the tube (in the region investigated ); 
authors found that the solution breaks down 


the tube. 
tion on the 
that is why 
behind the leading Mach line from the edge of the tube. 

The proper physical picture of this problem is shown in Fig 


some 
9 
The present author would suggest that a cylindrical source dis- 
tribution would correctly represent the actual flow. The cylindri 
cal source distribution is made up of ring shape source distribu- 
tion. The source strength is constant on each ring whose radius is 
roand whose center is on the x-axis. Each source on the ring will 
have the velocity potential, 


f® 
o= | fix BC ro 
J0 


— 7) cosh ul | du = 


Px 
fix + (Br - 
0 


It is evident that the Mach line is x + 8r = constant and the re- 
Furthermore, if one 


Bro) cosh u|] du (7) 


gion of influence is downstream (Fig. 3). 
chooses 79 > R, where R is the average radius of the tube, no arti- 
ficial sources are introduced in the region where the flow is investi- 
gated. 
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The application of this method of ring shape source distribution 
to axially symmetrical supersonic flow will be reported later. 

A similar attempt of this method was given in reference 2, but 
our jurmulation of the formula is different. We use a method that 
is a modification of Weinstein’s method of ring shape source in 


potential flow.* 
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The Effect of Surface Cooling on Laminar 
Boundary-Layer Stability 


Martin Bloom 

Research Associate, Polytechnic Institute of Brooklyn Aeronautical 
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May 10, 1951 


i be STATEMENT IS MADE in reference 1 that the minimum 
critical Reynolds Number for laminar boundary-layer sta- 
bility is increased in a pronounced manner by surface cooling, 
especially at supersonic speeds. This is based on results of 
reference 2, which state that, at a free-stream Mach Number (1/7) 
of 4, complete stability in an isobaric flow can be achieved by 
cooling the surface to 0.86 of its equilibrium temperature (7,) 
1, that the 
derived therein for adiabatic flow are conservative 


Polytechnic 


On this basis, it is estimated in reference results 


An investigation, carried out at the Institute of 
Brooklyn under the sponsorship of the N.A.C.A., 
wise. Results of this study (Figs. 


ing has little stabilizing effect on an isobaric flow in the supersonic 


indicates other- 
1 and 2) show that surface cool- 
range, so that the calculations of reference 1 are less conservative 
than it is believed. 
plete stability can be achieved by surface cooling for 1 < 


The accompanying figures show that com 
M < 5.7 
Moreover, whereas surface cooling is effective for stabiliza- 
tion at WV < 1 
the amount of required cooling increases sharply with increasing 
M = 4, the surface 
temperature required for any practical amount of stabilization is 
namely, 0.257, 


only. 
5, its effect at higher Mach Numbers is small, since 
M in the supersonic range. In particular, at 
close to that required for complete stabilization 
These results seem reasonable in view of the known destabilizing 
effect of high Mach Numbers on adiabatic flows. Comparable re- 
sults of reference 2 appear to be in error in this respect; this is 


probably due to the assumption of a small value for a parameter 
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Fic. 2. Difference between equilibrium and surface tempera 
ture required to stabilize the laminar boundary layer by cooling 
at various free-stream Mach Numbers and Npr = 1. 


[Eq. (15), reference 2], which in these cases may be large, as well 
as to numerical errors. 

The stability calculations represented in Figs. 1 and 2 were 
carried out according to the theory developed in reference 2, 
modified somewhat for these applications. For purposes of com 
parison, mean-flow properties were calculated by two alternative 
methods—the sixth degree polynomial integral method (also used 
in reference 1) and the mathematically exact type of analysis of 
reference 3. The surface temperatures required to stabilize the 
flow completely and to stabilize the flow only a small amount 
(Re,er /C = 500) were plotted. The symbol @ denotes momentum 
thickness, and C denotes a constant, of the order of one, which is 
used to approximate linearly the Sutherland temperature-vis- 


cosity relation.* 
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Astroballistic Heat Transfer* 


R. N. Thomas and F. L. Whipple 
University of Utah and Harvard College Observatory 
June 1, 1951 


ABSTRACT 


Some current work on heat transfer to rapidly moving objects is summa 


rized The work embodies free-flight study of ablation 


INTRODUCTION 


W: PRESENT HERE A BRIEF RESUME of some recent work on 
heat transfer in the astroballistic régime. We use heat trans- 
fer in the conventional aerodynamic sense of energy received by 
the body from the original translational plus internal energy of the 

* An investigation supported in part by Naval Ordnance Contract No 


10449 Task I, of the U.S. Naval Bureau of Ordnance and in part by the Arma 
ment Branch, Office of Naval Research 


gas. The measures reported, however, relate to an inference of 
heat transfer from ablation measures. We have introduced the 
term “‘astroballistics”’ to define that region where the heat transfer 
to a rapidly moving object is so great as to cause appreciable 
ablation from the surface of the body. The results on the heat 
transfer may probably also be taken over directly to some cases of 
nonablation; at the least, the results suggest the desirability of re. 
assessing current methods used in calculating such heat transfer 
The following is a summary of results presented from the engi- 
neering-aerodynamic-ballistic aspect. The complete details of th, 
work, presented from the standpoint of the application to the 
theory of meteors, will appear partly in a forthcoming article jy 
the Astrophysical Journal and partly in a later paper. 


SUMMARY OF RESULTS 


The customary aerodynamic-ballistic representation of heat. 


transfer calculations (per unit area ) is based upon the expressions 


heat transfer = /i(t tn) 
hx 8 
= constant P,*R 9 
k 


where / is the heat-transfer coefficient, Rk is the thermal con- 
ductivity, x is the distance of the surface element from the tip of 
the body, ¢, is the air temperature outside the thermal boundary 
layer, t) is the body temperature, and ?, and R, are the Prandt| 
and Reynolds Numbers. In the present note, the significant 
quantity is 8, which is usually taken to have the value (0.5 for 
laminar, and 0.8 for the turbulent flow.' On the other hand, the 
heat transfer in the case of meteors and other phenomena occur 
ring in rareified gases is generally taken as? 
heat transfer = A p(v?/2) 3 
where A is a dimensionless heat-transfer efficiency coefficient, p is 
the air density, and v is the velocity. If A is a sufficiently slowly 
varying function, Eq. (3) is obviously greatly to be preferred to 
Eq. (1), since no knowledge of body or gas temperature is re- 
quired. Both of these temperatures are usually difficult to de- 
termine and/or predict, particularly in free-flight experiments 
It is, of course, clear that the Eqs. (1) and (3) are similar under 
certain conditions at very high velocity, as may be seen by r 


writing Eq. (8) as 


nA Roe 
heat transfer = ~ 
2Cp x k M p,a* p 
(: 7 + zt P-R.) a 
Cy } M xv 


where we have used the standard formula for the pressure rise 
across a shock wave normal to the flow (for high velocity 
bi/p ~ [2¥/(v + 1)](v?/a? (5 


where kp is the Boltzman constant and V/ is the molecular weight 
of the gas. Thus, in the case that fy > t and 


cpaP,)-*R,-? 


Eqs. (1) and (3) are seen to be identical. There is, however, no 
reason to assume such a law of variation forcy. The fundamental 
experimental difference between Eqs. (2) and (3) lies then in the 
variation of heat transfer with air density 

We have performed experiments on free-flight heat transfer, 
using the controlled-pressure free-flight range of the Naval 
Ordnance Laboratory through the kind cooperation of the Aero- 
ballistics Department, under R. J. Seeger, with the Harvard Ob- 
servatory Meteor program. We denote these results by the 
letters HMP-NOL. We have also a few results on some CO: 
model wind-tunnel experiments, made available to the Meteor 
Program through the N.O.L. Aeroballistics Department through 
H. H. Kurzweg and J. M. Kendall. Some experiments on very 
high velocity pellets have been performed by J]. S. Rinehart and 
his group at the Naval Ordnance Test Station at Inyokern to pro- 
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TABLE | 


Experimental Results on the Heat-Transfer Efficiency Factor 


Logw ¢ 
(Minimum) Range of 2 Range of logio p 
Ergs Gm. > A (Km. Sec Gm, Cm 
Free-flight 3.92 to 
HMP-NOI 0012 0.76 to 1.8 2.62 
Free-flight 0.012 « 
N.O.T.S A < 0.024 2.0 to 2.5 2.92 
Free-flight 
N.O.T.S A < 0.04 v< 5.0 2.92 
Wind-tunnel 0.014 0.48 to 0.67 —3.48 to —4.57 
12.2 A > 0.05 8 to 65 -5.5to —9.3 


Meteor 


vide a third direct measure of A. Rinehart and the Mechanics 
Branch of N.O.T.S. have kindly permitted us to publish an ab- 
stract of these results here to complement our own measures and 
as part of some cooperative research between N.O.T.S. and the 
University of Utah currently under way. We denote these re- 
sults by N.O.T.S. All these exveriments refer to direct observa 
tions of mass loss by heat transfer resulting in ablation. Finally, 
we have analyzed the considerable material now available from 
the Meteor Program being conducted by the Harvard Observatory 
under a contract with the Naval Bureau of Ordnance. From the 
free-flight and wind-tunnel experiments we obtain Ccirectly the 


quantity A. The meteor observations provide only values of 


A/2T¢t (6) 


c= 
where ¢ is the heat of vaporization or fusion according to which 
process gives the mass loss and where T is defined by 


drag force = I'pv?A (7) 


with A the projected area perpendicular to the velocity. (IT and 
A are thus directly efficiency factors expressed relative to the 
Newtonian putty ball. We regret the introduction of two new 
coefficients into an already confused ballistic-aerodynamic litera- 
ture, but these coefficients have been in standard meteoric use.) 
We conclude that the best value of ¢ and IT for the meteor case are 
7.10" ergs Gm. ~! (representing vaporization from iron meteoroids ) 
and 0.5, respectively. There is considerable dispersion in the o 
value obtained from meteors, but no obvious correlation with 
velocity or air density. 

We tabulate the mean values of A for the free-flight and wind- 
tunnel experiments and o and A from the meteor experiments, in 
Table 1. We give also the range of v and p covered by the experi- 
ments. 

A combination of HMP-NOL and N.O.T.S. values may provide 
one more bit of evidence on the internal consistency of the results 
In the HMP-NOL experiments 
melting was observed at 


under certain assumptions. 
carried out on woodsmetal projectiles 
velocities around 3,500 ft. per sec. and not observed at velocities 
around 2,500 ft. per sec. It can be shown theoretically, and has 
been demonstrated experimentally in lined shaped charge experi 
ments, that luminosity produced by rapidly moving bodies comes 
mainly from the interaction of the ablated material with the at- 
mosphere and not from the solid body itself. At velocities high 
enough to melt iron, the iron oxide formation should, at least, be 
visible. If then we use Eq. (1), we have 
mass loss/unit area = (A/f)(pv*/2) 

Thus, assuming that iron should become luminous when it begins 
to melt and taking ¢ = 0.5 * 10° ergs per Gm. for woodsmetal 
and 1,10" for iron, we predict that an iron missile should become 
which woodsmetal 
(We have 
Lumi- 


luminous at velocities about (20)'/? those at 
began to melt or between 6,800 and 9,500 ft. per sec. 
assumed A constant so the estimate is a lower bound.) 


uosity is observed to begin at about 8,000 ft. per sec. 
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If we attempt to represent the free-flight data by formulas of 
the type of Eqs. (1) and (2), using computed temperature behind 
the shock wave for fj and room temperature for f, we obtain 
(HMP-NOL free-flight results alone) 


B = 1.06 + 0.15 (std. dev.) 


From the wind-tunnel results (only two measures were available), 
0.97. From a combination of wind-tunnel and free-flight re- 
sults, B = 0.95. 

We conclude, therefore, that the data strongly suggest that Eq. 
(3) provides a satisfactory representation of the heat transfer over 


B= 


a sufficiently wide range of velocity and air density to be generally 
adopted for current work in high-speed aerodynamics and ballis- 
tics, as well as meteorics and astroballistics. We wish, however, 
to emphasize the desirability of further free-flight experimentation 
on this heat transfer, for our data are far from complete. With 95 
pe® cent of current research in aerodynamics and ballistics con 
fined to wind tunnels, we cannot overemphasize this need for free- 
flight experimentation, particularly in the astroballistic régime, 
since wind-tunnel measures cannot induce mass loss from ma- 


terials of engineering interest. 
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On an Extension of Glauert’s Theory of Thin 
Airfoils to Profiles in Cascade 


N. Scholz 
Dr.-Ing., Institute of Fluid Mechanics, Technical University, Braun- 

schweig, Germany 
May 12, 1951 

N A FORMER NOTE! we have shown how the two-dimensional 
I potential flow through a cascade of airfoils can be calculated 
by replacing the airfoils by a distribution of vorticity and sources 
and sinks (method of singularities). We dealt mainly with the 
“Indirect Problem,’’ where the distribution of vorticity (and 
sources and sinks) is given and the shape of the profile in cascade 
is asked for. More important but also mathematically more dif 
ficult is the ‘‘Direct Problem,’’ where the shape of the airfoil in 
cascade is given and the distribution of velocity over the surface 
of the profile is asked for. In this note we want to contribute to 
this ‘‘Direct Problem. 

For the very thin single airfoil there exists the well-known solu- 
due to Glauert,? who replaces the 


” 


tion of the “Direct Problem, 
airfoil by a vorticity distribution that is expanded into a Fourier 
This method has been further extended by Allen* and 
We want to show that the method 


series 
Riegels‘ for the single airfoil. 
of Glauert can also be applied to a cascade of airfoils 


Each airfoil of the cascade in the complex plane = = x + iyis 
replaced by a vorticity distribution along its chord (see Fig. 1 
Then the induced complex velocity distribution w = u iv on the 


chord of the profile (x-axis) is, according to reference 1, 


tr es 9 » at 
r c/2 , 4 v P 

f ¥(x’) coth {| we” ) di (1 
2 Jx’=—c/2 f 


chord c, t/¢ 


wx) = 


Here, y(x) is the vorticity distribution along the 
means the solidity, and \ is the angle of stagger 

t—> o, one obtains from Eq. (1) the well-known equation for the 
The shape of the profile 


In the limit, 


induced downwash of a single airfoil. 
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y(x) being known, the induced velocity w(x) is also known, so 
that Eq. (1) is an integral equation for the unknown vorticity dis- 
tribution y(x). In the same manner as in Glauert’s theory of the 
single airfoil, the solution of the integral Eq. (1) can be given in 
the form of a Fourier series for y(x). Replacing x, which is run- 
ning from —c/2 <x < +c/2,byacomplex angular coordinate ¢, 

—cos ¢ = tanh [(rx/t)e™]/ tanh [(rc/2t)e] (2) 
yields a real part of g running from 0 to z, whereas the imaginary 
part of ¢ is zero for x = +c/2. With this coor- 
dinate, ¢, the Fourier expansion of +(x) can be given by just the 


—c/2and x = 


same series as for a single airfoil, due to Pistolesi,® 


¢ 
ao cot = 


v(¢) = Ua + » 3 a, sin (ng) (3) 
n=1 


Here, the Fourier coefficients do, a), . . ., and also y(¢) are complex 
quantities, but only the real part of y(¢) has a physical meaning. 
Introducing Eq. (3) into Eq. (1) and integrating by terms yields 


for the complex velocity, 


do me ~ 
‘ Xr a 
re jg cosh ae" —_ > “cos (ng) + 
2 2t 2 


n=1 < 


—~ + a, Lat 
p Ds tanh” e (4) 
n=2 2 4t 


(Here >* means that only the even terms have to be summed up.) 
By the aid of the series expansion of the hyperbolic functions, it 
can easily be shown that, in the limit t ~ © (single airfoil), Eq. 
(4) becomes purely imaginary and results in the well-known 


wx) = 


Glauert downwash equation 


ao ay, 
W(x) = Uo ‘OS 5) 
= + >» , cos (ng) (5 
2 n=1 2 
As in Glauert’s equation [Eq. (5)], as well as in our Eq. (4), the 


first term yields an induced velocity that is constant along the 
chord. It means the induced velocity due to the vorticity dis- 
tribution of a flat plate in cascade; but, furthermore, in Eq. (4) 
the even sin terms of Eq. (3) by the terms tanh” (e*2c/4t) also 


give a contribution to the constant part of the downwash. The 
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physical meaning of this is that, for a cambered airfoil in cascade 
the “ideal angle of attack’’t depends on the solidity. 

In order to calculate the unknown coefficient do, a, . . ., we tise 
the kinematic flow condition that the direction of the resulting 
velocity is tangential to the profile. Since in our linear theory the 
induced velocity in the x-direction can be neglected, the kine. 
matic flow condition in complex velocities is 

Us — 1U2(dy/dx) = Us + w(x) (6 


This is equivalent to (dy/dx) = o(x)/U.. 
Eq. (4) into Eq. (6) yields the coefficients a, in the same manner 


as in Glauert’s theory by simple quadrature 


4 4 dy ( j 1.2 
1, = cos(n Ge, n=I1,Z,... 7 
: rw? dx re ? 


Because of the last sum in Eq. (4) the coefficient a9 depends on the 


Introducing w(x) from 


coefficients a,, » > 2. Therefore it is necessary to calculate a, 


last: 
I 


~ cosh [(rc/2t) e* |] 


2 (x (dy i alt » 
Ba (2) ae 4 > a, tanh (= e )] (8 


Eqs. (7) and (8) result in the well-known 


ag = 


In the limit t ~ o~, 
formulas of a single airfoil due to Glauert. As in the case of a 
cascade, Eqs. (7) and (8) are complex; they have to be split up 
into two real equations for numerical evaluation. To each value 
of the coordinate x there correspond a real and an imaginary part 
of the complex angle ¢, which can be calculated universally from 
Eq. (2), not depending on the shape of the airfoil but only on the 
solidity ¢/c and the stagger A. After having evaluated the real 
and imaginary parts of ao, . . ., a, from Eqs. (7) and (8), the vor 
ticity distribution can be calculated as the real part of Eq. (3 

The final result is the velocity distribution on the chord of the 


cascade airfoil 


W(x) = Uo + u(x) + [y(x)/2] (9 


with + for the upper and — for the lower side of the airfoil. For 
the cascade airfoil there is, in the general case, an induced velocity 
in the x-direction on the profile, u(x), which has to be determined 
as the real part from Eq. (4). 

The total circulation of one cascade airfoil is obtained from 
Eq. (3) on integrating in the complex form 


TC . 
Uet | aoe~*® sinh (= >) + 


Dk 
° ugé 
e~™ > a,, tanh” (« ) 
2@ i 


n= 


(10 


where, again, only the real part has a physical meaning. In the 
limit > ©, Eq.(10)also results in the well-known Glauert equa 


tion 


TY = Uec(r/2) [ao + (1/2)a] (11 
For staggered cascades with highly cambered profiles it may be 
necessary that, in the kinematic flow condition, Eq. (7), the in- 
This leads toa 


nonlinear relation between the profile shape y(x) and the vor- 


duced velocity in the x-direction, u(x), is retained. 


ticity distribution y(x). 


1 This being defined as the angle of attack at which the flow enters the 


leading edge smoothly 
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Past an Infinite Row of Thin 
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The Flow Around a Rotating Disc in a Uniform 
Stream 


H. Schlichting and E. Truckenbrodt 
Prof. Dr. and Dr.-Ing., Respectively, Institute of Fluid Mechanics, 
Technical University, Braunschweig, Germany 


May 12, 1951 


OR LAMINAR FLOW the solutions of the following two problems 
pee known: 

(1) The flow in the vicinity of an infinite plane rotating in a 
fluid that is at rest far from the plane (von Karman! and Coch- 
ran’) 

(2) The flow with axial symmetry normal to an infinite plane 
flow at a stagnation point—(Homann‘). 

The solutions given in references 2 and 3 are exact solutions of 
the Navier-Stokes equations. 

We have solved the more general problem of the laminar flow 
with axial symmetry normal to an infinite plane that is rotating 
with constant angular velocity about an axis normal to the plane. 
This case contains both problems (1) and (2) in the limit. If 
r, y, and z are cylindrical polar coordinates, the plane is taken to 
be z = 0; it is rotating with constant angular velocity w about the 
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there is a boundary layer in the vicinity of the plane. If u, v, and 
w are the components of the velocity in the boundary layer 
and U, V, 


the potential flow outside the boundary layer in the directions of 


and W are the corresponding velocity components of 


r, ¢, and s, respectively, the potential flow is given by 


—2az (1) 


Uz=ar, V=0, We= 


with a as a constant, which gives the intensity of the potential 


flow. 
According to the solution of (1), the thickness of the boundary 
layer is 65 ~ V v/w; whereas, according to (2), 6.5~ VY v/a. In 


our general case the thickness of the boundary layer 6 depends on 
v/aand a/w. We have calculated the velocity distribution in the 
boundary layer from the momentum theorem, which for the 7 and 


¢ direction yields the following two equations: 


d ie : dU ; 
r u(U — u)dz|+r (Ll u) dz 
dr 0 dr Jo 
Me r, 
vds=r (2) 
0 P 
d ; » Te 
rv? (uv) dz | = r? — (3) 
dr 0 p 


Here, 7, and rg are the components of the shearing stress at the 
disc, z = 0, in the direction of r and ¢, respectively. Eqs. (2) 
and (3) have been solved by taking, for u and v, polynomials of 
6, which fulfill the 
most important boundary conditions at the wall and at the edge 
We have taken a polynomial of fourth 
The polynomial for u con- 


the dimensionless distance from the wall, 


of the boundary layer. 
degree for v and of fifth degree for uw. 
tains an arbitrary constant that, together with the boundary-layer 

































axis r = 0. From the solutions of (1) and (2), it is known that thickness 6, is left to be calculated from Eqs. (2) and (3), From 
80 7 
40 CulRe 
60 ‘ 
20 40 
<p exact, Cochran [. 2) 
10 | 
. na ew oom wee le 
z]83 ~) 
i . Fig.1: Moment Coefficient of a rotating aisc 
" as a function of a/a 
44 
laminar. + 
exact, Cochran[2} | 
| 
2 | 























Moment Coefficient 100 Cn 


~— 
| 
—$——}—}—__ 
























































ge} - 
a6;— ar ; | 11 | turbulent, 
= —— ee t TT - Karman [1]- 
nn Sa SS = | eared 
| | 
| | ges i | || 
02 ; + | + ; + 1 } | ; i LJ r = = 
70 2 4 6 810 2 4 6 8 2 4 6 8 10 2 
Fig.2: Moment Coefficient of a rotating disc Reynolds Number Re = Ree 


for «various We/Re . 





640 JOURNAL OF THE 


these calculations one obtains the boundary-layer thickness 6, 
the velocity distribution in the boundary layer, and the com- 
ponents of the shearing stress 7, and ry, all as functions of the 
dimensionless quantity a/w, reducing to the known solutions of 
(1) for a/w — O, and of (2) fora/w— ~. 

From 7, the total moment J/ for one side of a disc of radius R 
can be calculated. Defining the moment coefficient by 


Cu = M/(p/2) w?R® (4) 


yields Cu'V Re as a function of a/w as given in Fig. 1, with 
Re = R*%w/v as the Reynolds Number. The dimensionless quan- 
tity a/w of the infinite rotating plane can be related to Wa2/Vo = 
W./Rw of a dise with radius R, where W is the velocity normal 
to the disc and V) = Rw the velocity of the edge of the disc due to 
rotation. One gets, from the condition of similarity of both flows 
in the vicinity of the stagnation point, 

a/w = (2/7) (Wao/Rw) (5) 
Therefore, the moment coefficient can also be given as a function 
of Re and W../Rw (see Fig. 2). For Wao/V» > 1, the results can 


be interpolated by 
/ / 
Cu = (2.53/V Re) V Wa/Vo (6) 


For constant angular velocity the moment increases considerably 
with increasing Wo. 
More details of our calculations will be published shortly.# 
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Erratum 


T. H. H. Pian 
Aero-Elastic and Structures Research Laboratory, Massachusetts 
Institute of Technology, Cambridge, Mass. 


May 28, 1951 


HE FOLLOWING ERRORS appeared in the paper ‘‘Analytical 
Bsa Experimental Studies on Dynamic Loads in Airplane 
Structures During Landing’ (JOURNAL OF THE AERONAUTICAL 
ScreNcES, Vol. 17, No. 12, pp. 765-775, December, 1950) by 
H. I. Flomenhoft and myself. 

The expressions for Bo’s in Eqs. (25), (26), and (27) should be 
as follows: 
Equation (25): 


*b/2 
J eM i. P 
g i" S$ — ma ~ e)l> de, x8 
a 4 ym 1, | ma €) ( dx, xX > xp 
) 


m(a — e)] bax ie u) 


xX <x 


Equation (26): 


a ie eM. 
Bo = Mmaz. 2 <m + IS 
6/2 b 2\ . 
x > 


nd eM. | 

Bo = Mmaz- & ym + [S — m(a — e)] ¢ dx dx - 
b/2J »/2\ Fr. f 

»), x <M 


|+ sf dx, x > xr 


- e)] t dx dx, 


l 
M(xr — 
9 
Equation 


n 


max 


20. ; 
= Nmar. 8 ' [I — Sla -e)] + Stas — 2 aa) 
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Discussion of Wright Brothers Lecture 


(Continued from page 623) 


represents a small enough investment so that it can be 
assigned to a single program; yet it is sufficiently 
flexible so that it can accommodate changes in aero- 
dynamic and control parameters and can include cer- 
tain types of nonlinearities. 

Frequently the noise problem is so severe that it 
must be properly taken into account even before the 
control system can be made to operate. In such cases, 
a noise study does not lead to a refinement; rather, 
it dictates the control system that will give satis- 
factory performance. Hence, the noise study be- 
comes part of the theoretical design 


In conclusion, I want to commend Dr. Bollay for 
his statement that much closer teamwork is required 
the 
signer. The 
much of the 
control system complexity will involve increased costs 
and decreased reliability. Going back to Wilbur 
Wright’s address in 1901, we find that he was con- 
“steering.” 


between aerodynamicist and the autopilot de 
aerodynamicist should not expect too 


autopilot designer, otherwise resulting 


cerned with the problems of power and 
Enough power is currently available to fly a barn door. 
However, let us not expect the autopilot designer to 
make it maneuverable. 
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